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Preface 



Very recently, the PVLAS collaboration has reported the observation of two unexpected effects. 
Studying the propagation of linearly polarized laser light through a strong transverse magnetic field in 
vacuum they find an anomalous rotation of the polarization plane as well as an induced ellipticity of the 
outgoing beam. None of these two effects has a standard explanation within conventional physics at this 
moment, but they converge into a frequent prediction of physics beyond the standard model of particle 
physics: the existence of neutral and nearly massless bosons coupled to light. 

The key word in the above paragraph is "unexpected" . Indeed, the small collaboration was actively 
looking for these so called axionlike particles (ALPs) leading to the two mentioned effects. The important 
point is that these particles seem to have couplings to the electromagnetic field which are far to be allowed 
by other experiments and, in particular by astrophysical arguments. 

A conservative look would then discard immediately the ALP interpretation of the PVLAS signal but 
I shall show in these few pages that this decision would be a bit premature. In these years under the 
supervision of Dr. Masso we have showed that there are particle physics models in which the discrepancies 
coming from astrophysics and other laboratory experiments can be circumvented. Notably, these models 
involve new physics at low energy scales. In this PhD thesis I address the issue of presenting all our 
conclusions. 

The realization of this work has been a great pleasure for me, in part because of the confluence of 
the many fields of physics involved, optics, astrophysics, cosmology and, of course, particle physics, but 
mostly because it has been carried out in friendly touch with many experts on all these fields. It has 
been then twice an enriching experience also, but moreover, I think that has given our work a soundness 
which otherwise would be lacking. The multidisciplinarity of this thesis has a natural drawback, however, 
one would be tempted to review all the fields involved and produce a too long writing. I was not really 
seduced with this idea and I have tried only to touch the points I really consider important to understand 
my original work. Wherever I feel there is more to be said I give the needed references. 

Since the publication of the PVLAS result we have been witnesses of a real "boom" in the field 
with both, theoretical efforts and the advent of many experimental proposals to test the PVLAS ALP 
interpretation. This reflects, among other things, the great importance of having new tests of physics 
beyond the standard model. In the years before the Large Hadron Collider, which will scan TeV energies 
with an unprecedent budget in the history of physics, we could be finding a new frontier of knowledge in 
a much more modest experiment, and at very low energies. 

In Chapter Q] I briefly motivate low mass particles coupled to light and introduce the required ingre- 
dients that will show up in my models. I am concerned with scalar particles such as pseudoGoldstone 
bosons, candidates for the PVLAS ALP, paraphotons and millicharged particles. Chapter [2] is devoted 
to introduce the PVLAS experiment as well as its ALP interpretation in terms of a new scalar particle 
coupled to two photons. Further comments will be made about other possible interpretations. Finally, in 
Chapter [3] I review the experimental knowledge about the low mass particles coupled to light that could 



be involved in the PVLAS signal and state the inconsistency of the bare PVLAS ALP interpretation with 
astrophysics and two very sensitive experiments: ALP helioscopes and 5 th force searches. 

Chapter 4 constitutes an introduction for the original work presented in this thesis. My contributions 
are organized as a compendium of articles published between the years 2005 and 2007, each in a separate 
Chapter. The chronological order in which they appear reflects somehow the evolution of the theoretical 
efforts required during these last three years. Chapter 5 (reference pQ) gives two general ideas for 
reconciling the ALP interpretation of PVLAS with the physics of stars introducing new particles and 
interactions at low, accessible, energy scales ~ eV. The concrete models proposed seem to be now a bit 
obsolete but the message imprinted remains valid, there is need for additional low energy physics in order 
to evade astrophysical bounds on the PVLAS particle. 

A further refined model is presented in Chapter 6 (reference [2]). There we include two paraphotons 
and a millicharged particle, together with the PVLAS ALP, in order to circumvent the astrophysical 
constraints. The model is simple and, although it has some fine-tuned quantities, turns out to be one of 
the few models that could to the job if the PVLAS ALP is confirmed. 

Even before the time this paper was published, the community was starting to be active proposing 
other models and we had the feeling that some model-independent study of evading astrophysical bounds 
could be of interest. Therefore we prepared [3], which is presented in Chapter 7, as a tool of evaluating 
more precisely the models in [IJ [2] (which lack precise calculations) as well as any other of similar 
characteristic^. The general conclusion is again that the new physics should appear at scales accessible 
for Earth-ground laboratory experiments (even lower than naively expected) and therefore experimental 
efforts were strongly encouraged. 

Last summer, the PVLAS collaboration released new data suggesting the parity of the ALP. Mea- 
surements performed in gas indicated that the ALP should be a parity-even particle and therefore it 
would mediate long range forces between macroscopic bodies. In the last article presented in this thesis 
(Chapter 8 from reference [1]) we analyzed such a possibility. The very precise measurements of the 
Newton's law and Casimir effect rule out contributions of the PVLAS particle and the impact in the 
parameter space of ALPs is really impressive. Moreover, we realized than in the model presented in 
Chapter 6, the force is suppressed with respect to the naive expectations and the ALP interpretation of 
PVLAS is possible. 

Finally, In Chapter 9 a short summary and a final discussion are presented. 



1 Interestingly enough, we noticed a new way of evading the astrophysical bounds if the PVLAS particle is a Chameleon 
field. 



Chapter 1 

General concerns about light weakly 
coupled particles 

1.1 Introduction 

Despite the considerable success of explaining almost every observation within particle physics, the stan- 
dard model (SM), based on the gauge symmetry SU(3)c <8> SU(2)i <8> U(l)y, is generally thought not to 
be the final theory explaining the interactions between elementary particles. There are several reasons 
for that, and we can arrange them into three big sets of problems and objections. 

The first one is concerned with very fundamental objections to the quantum theory: the postulate of 
measurement, the lack of a direction of time, etc. The second is less fundamental and simply concerns the 
elusion of the gravitational interactions. This problem is addressed at present within quantum gravity 
and string theory, being the last the most promising but still far from providing mature answers. 

Finally, the third set are problems regarding the concrete form of the standard model as a quantum 
field theory. Most of them can be regarded as purely "aesthetical" features. As examples we find the 
large number of parameters (coming from our ignorance of a possible flavor structure at higher energies?), 
the flavor problem (why three generations?), the hierarchy problem (why such a difference between the 
electroweak (EW) and the Planck scales?), the missing Higgs boson (where is it?), the strong CP problem 
(why does QCD respect CP?), etc. 

While the first two sets of "fundamental problems" of the SM seem difficult to solve in a quantum 
field theory framework, those from the third are the subject of many speculative but ingenious proposals 
beyond the SM. For instance, the existence of a symmetry between bosons and fermions, the so called 
Supersymmetry, could help to stabilize the EW-Planck hierarchy and theories of dynamical breaking of 
the EW symmetry do not need a Higgs boson. These models generally introduce new fields and gauge 
symmetries beyond the SM. Most of the new particles and gauge bosons are extremely massive and 
we only have chances to discover them by building huge accelerators or by precision measurements of 
carefully chosen observables. 

However there are several of such theories in which also low mass particles are predicted. If these 
particles exist they must necessarily be very weakly interacting, otherwise they would have been already 
discovered.]. On the other hand they can be active and stable at the low energies of our universe and 
give rise to a completely different phenomenology from their massive companions. 

This Chapter is devoted to review some relevant scenarios where those particles arise while I leave 
the study of their peculiar phenomenology to Chapter [3l Particular emphasis will be made on their 



1 This last assertion acts as a definition of "very weakly interacting" and "low mass" in a precise sense. 
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interactions with the electromagnetic field because these turn out to be specially interesting from the 
experimental point of view. Indeed many of the experimental efforts looking for these particles have been 
focused on these couplings, and one of those experiments, performed by the PVLAS collaboration, has 
recently reported a positive signal. 

From a theoretical viewpoint, the existence of low mass particles is closely related to symmetries. 
Regarding the renormalization group approach, masses are dimensionful parameters of the theory that 
tend to develop values near the ultraviolet cut-off unless some symmetry forbids them. Put it another 
way: if we intend to force some dimensionful parameter to be much smaller than the higher energy 
scale in the theory A, radiative corrections involving particles related to A will induce contributions 
proportional to powers of A which will have to be canceled at a very precise level with counterterms of 
the same order. Such a tuned cancelation is generally though to be unnatural and unaesthetical, and 
it engenders the so called hierarchy problem. The only way to control these radiative corrections is to 
protect the small scales with symmetries. Indeed, supersymmetric theories were motivated in part to 
avoid the electroweak-planck hierarchy. As notorious examples we find also the masses of fermions in 
chiral theories and Goldstone and gauge bosons. Let me say a few words about them. 

The lowest-dimension (non trivial) irreducible representations of the Lorentz group in 3+1 space-time 
dimensions are spinors of left and right handed chirality. Parity transformations convert one set into the 
other but, although parity is approximately conserved, it seems not to be a good quantum number at 
the most fundamental level. It is allowed then for LH and RH fermions to be in different representations 
of the gauge group, whatever they are. Dirac mass terms involve just one LH and other RH spinor so 
they have to combine forming a gauge invariant term. From this viewpoint, a pair of LH and RH fermion 
fields living in representations that allow a scalar of the gauge group will acquire a mass of the order 
of the cut-off, A, but LH and RH fermion fields which cannot be paired for any reasoro would lead to 
massless modes. On the other hand, Majorana mass terms are allowed only for fields completely neutral 
under the gauge group. This is the case, for instance, of the RH neutrinos of the standard model which 
for the above reasoning are supposed to be very massive. Interestingly, a huge Majorana mass for the RH 
neutrino field produces a suppression of the Dirac neutrino masses by means of the see-saw mechanism. 

This picture is valid as long as the gauge symmetry is unbroken. In the standard model, the Higgs 
field acquires a vacuum expectation value that breaks the SU(2)j U(l)y sector to U{V)qed- Li this 
way, interactions of Yukawa type of the Higgs with a pair (LH-RH) of fermion fields lead to effective 
mass terms that are naturally much smaller than the cut-off scale of the theory as long as the vacuum 
expectation value of the Higgs, or the Yukawa couplings, are small. Explaining why a scalar field like 
the Higgs could develop a small expectation value ~ 200 GeV (small with respect to the Planck scale 
mpi ~ 10 19 GeV) is again the hierarchy problem. 

Furthermore, the Nambu-Goldstone theorem tells us that theories with a global symmetry sponta- 
neously broken in the vacuum contain massless particles, known as Goldstone bosons. This result goes 
beyond perturbation theory so radiative corrections do not spoil this conclusion. Nevertheless most of 
the known examples of Goldstone bosons come from slightly explicitly broken symmetries which allow 
for small masses that can be again sensitive to the hierarchy. 

Finally, the gauge bosons of a local unbroken symmetry are naturally massless at all orders being 
their masses protected by the gauge invariance of the dynamics. There is a nice exception in pure U(l) 
symmetries because a mass term is also allowed [5j El [TJ [8] . From this point of view it seems natural 
to think that U(1)qed or U(l)y are not fundamental U(l) symmetries but the low energy residual of a 
higher embedding non-abelian group. 

In the remaining of this Chapter I will go deeper into the important theoretical frameworks that 



J For instance, there could be more LH than RH fields in nature, they could belong to representations of different range, etc. 
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provide small mass particles that are crucial for this thesis. I begin with an exposition of the Nambu- 
Goldstone theorem followed by the related history of two of the most important examples, pions coming 
from the breaking of the axial part of the isospin symmetry of QCD and axions, Goldstone bosons of a 
hypothetical broken axial U(l) symmetry proposed to solve the strong CP problem. The discussion is 
focused on this last case, since it has become a kind of paradigm in low mass bosons searches. 

Next I consider the phenomenology of new U{\) symmetries in a hidden sector, giving rise to parapho- 
tons and millicharged particles. These provide the framework for our solution to the PVL AS- Astrophysics 
puzzle so I describe them with greater detail. 

1.2 The Nambu-Goldstone theorem 

The Nambu-Goldstone theorem states that: 

whenever a global continuous symmetry of the action is not respected by the ground state, there appear 
massless particles, one for each broken generator of the symmetry group. 

The demonstration of this result inhabits central pages of many text books, so I believe it is not 
worth to be reviewed here. However, it will prove convenient to express it in mathematical form. To do 
so, consider the global infinitesimal transformation 

fo(x) -> <k(x) + Oatijfoix) (1.1) 

where <f>i(x) are fields (of any type, although only scalar fields are allowed to have nonzero vacuum expec- 
tation values), a a is an array of arbitrary infinitesimals, and t"- are the generators of the transformations 
with the index a running over such transformations. If these transformations are symmetries of the 
lagrangian there are associated currents J a ^ which are conserved in the sense that d^J^ = 0. We can 
write the Goldstone theorem^ as 

([J°»(x), <MO)]>o = / (03 [P»P?(P 2 V P - X + h.c] (1.2) 

i 

We see that the Kallen-Lehmann spectral density p?(^ 2 ) contains a pole in fi 2 = only if the associated 
symmetry generator is broken in the vacuum, ^2it%(<f>(0)j)o ^ 0. The corresponding massless particles, 
known as Goldstone bosons (GB), are created by combinations of the fields that will develop a v.e.v. 

7rf(s)cxt^(x). (1.4) 

It is clear that there are as many GBs as broken generators. The contribution of the massless modes 
to the two-point function f|1.2[) is proportional to (0| J a/i (x)|7r")(-7r"|0j(O)|O). A nonzero value requires 
that nf must be spin-zero fields (because 0j(O)|O) is rotationally invariant) and have the same conserved 
quantum numberg^l than J a0 (otherwise (0| J a0 (x)\irf) vanishes^). 

1.2.1 Approximate symmetries and pseudo Goldstone Bosons. 

It is usually the case that global continuous symmetries present in our theories are just approximate. 
The potential V(0) is then divided into a "symmetric" part Vq(4>) and a breaking term V± (</>). If the 



I display the result for hermitian J£ and 0iS only for simplicity. 



4 The quantum numbers that are conserved also by the vacuum. 

""Provided -nl has spin-zero, fi/0 can not contribute because it transforms as a spin-1 field. 
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breaking term V\ is small we also find GB but the new terms can include non-zero (but small) masses. 
We name such particles pseudo-Goldstone bosons (pGBs). Also, typically, a reduction of the degeneracy 
of the vacuum configurations is found (the so called vacuum alignment). 

1.2.2 Axions and the strong CP problem 

Consider quantum chromo dynamics with only two flavors of massless quarks, u and d. The lagrangian of 
this theory has a global symmetry U(2)l <8> U{2)r that allows to perform rotations of LH and RH fields 
independently, namely 

• ^ Ta (7) ; (7)^ e lX * Ta ( 7 



dh J \ d L J \ dn J \ dn : . 

where T a = {1, t} (a = 0,1,2,3) with r the Pauli matrices and A a ' free parameters. This symmetry 
can be also decomposed into vector and axial parts SU(2)y (g> SU{2)a <8> U(l)v ® U(1)a- We can easily 
associate the SU{2)y part with the isospin symmetrjQ of the strong interactions. Moreover, U(l)y leads 
to conservation of the baryon number, which is an observed property of the strong and electroweak 

interactional!. 

At first sight there is no symmetry (not even approximate) in nature related to the axial subgroups. At 
low energies, a quark condensate forms making (uu)q = (dd)o ^ and thus breaking spontaneously this 
part of the symmetry. The axial symmetries must then show up in the appearance of the corresponding 
Goldstone bosons. This seems to be the case for the SU(2)a part because we observe in nature three 
light bosons, the pions, one for each broken generator of SU{2)a- But there is no such Goldstone boson 
for the U(1)a symmetrjo. Weinberg stated this puzzle as the U(l) problem and suggested that there 
could be no U(\)a symmetry in the strong interactions. 

The solution of this problem was found by 't Hooft who realized that the U{1)a subgroup was indeed 
badly broken in QCD due to its non-trivial vacuum structure |10| . To understand the main points, first 
realize that U(1)a suffers an anomaly due to a triangle loop diagram (Fig. II. ip with quarks circulating 
on it. This anomaly contributes to the divergence of its related current J^ with 

%Ja = ^G^G a , u (1.6) 

where Ga V is the gluon field strength, and Ga U = ^e^apGa is its dual with e^ va p the total antisymmetric 
tensor in 4-dimensions. {N is the number of quarks circulating in the loop, 2 in our case.) This itself does 
not violate U(1)a because (|1.6p amounts a total divergency which in principle could be naively removed 
(as it happens in quantum electrodynamics) 

^-G^G a , u = d,JC» (1.7) 

But the surface integral of the Bardeen current K^ can not be set to zero. This is understood when 
we realize that the QCD vacuum is degenerated being the multiple states gauge transformations of the 
Ga = configuration. They can be classified by the so called topological winding number n which 
characterizes the way these configurations behave at spatial infinity. Notably, the QCD lagrangian allows 



6 In the standard model this symmetry is only approximate because u and d quarks have different electric charge and a 
"small" mass difference compared with Aqcd , both due to the electroweak interactions. 
This is valid up to the electroweak scale where instantons alter the picture. 

s The lightest candidate was the rj, but it turns out to be too big. Weinberg showed [9] that the squared mass of the 
U(1)a Goldstone should be smaller than three times the square of the pion mass. 
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Figure 1.1: Fcynman graph for the color anomaly. The bosonic line labeled J^ stands for a fictitious gauge field 
coupled to J^. The gluon field is labeled g and it is assumed a sum over all fcrmions q carrying color charge. There 
is an additional graph with the fermions circulating in opposite sense. 

for non-perturvative solutions of the equations of motion, the so-called instantons that mediate transitions 
between these vacuum states. But these instantons vanish very slowly at spatial infinity (~ 1/r) so they 
give finite contributions to the integral of the divergence of Bardeen's current. To be concrete, the 
contribution of Ga U (n + ,n_,x), an instanton taking a vacuum state with winding number n_ to one with 
n + is 

otxd^ = ^-j d A xG^G a ^ = v (1.8) 

where v = n + — n_, the diference in winding numbers. Clearly, U{\)a is violated by these instantons. 
This solves the Weinberg U(l) puzzle but unfortunately (or not) it poses a new and very interesting 
discussion remaining nowadays yet unsolved: the strong CP problem. 

Note that any physical amplitude must take these instanton solutions into account in the path integral. 
For reasons of gauge invariance and cluster decomposition one choses a weight factor e~ t9u for any of 
those configurations. This is formally equivalent to include a new term in the effective lagrangian 

C e = e-^G^G a ^ . (1.9) 

Such a term violates parity and time reversal invariance but conserves charge conjugation invariance (and 
thus violates CP). The picture looks a bit more complicated if we take into account that quarks gain a 
mass after the spontaneous symmetry breaking of the SU(2)j ® U(1)y gauge symmetry of the standard 
model. In general this induces a complex mass matrix 

C-mass = -q iL Mijq jR - q iR Mlflj L (1.10) 

(q holds for up and down quarks here) that can be diagonalized by a biunitary transformation that 
involves precisely a U(1)a rotation (jl.5p with 

\ L =i,2,3 = ° : X o = "Ao = ^ArgDetM (1.11) 

This rotation modifies the measure of the quark fields in the path integral of the theory [11] introducing 
another term like (jl.9p in the lagrangiarjj, shifting the value of 9 and allowing us to define an effective 9 

9 = 6 + AigDetM (1.12) 

so actually the CP violating term in the QCD sector of the standard model is the sum of two very different 
contributions, 9 coming from QCD vacuunr dynanrics and ArgDetM from the electroweak breaking sector. 
There is no convincing reason why any of these quantities should be zero, and it is even more unlikely 
that being both non-zero they cancel exactly 10 ! |13| . 



See the exposition in |12j . 
10 If U(1)a would have been a true symmetry then this transformation would have had no consequences and the phases of 
quark masses would be unobservable as well as the term (|1.9[) which could be rotated away. 
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However, as mentioned before, a nonzero value of 9 implies CP violating observables related to the 
strong interactions and, remarkably, it induces an electric dipole moment for the neutron which up to 
now has eluded detection. The latest measurements [2] set the impressive bound 

#<1(T 10 , (1.13) 

which finally poses the strong CP problem: 

why is 9 so small being the sum of two, in principle unrelated, quantities? 

Trying to solve the strong CP problem has motivated a lot of work and attractive speculations 
|13| . Here I focus on the proposal of Peccei and Quinn [15^ [16] which has a relevant phenomenological 
consequence for this work: the introduction of a pseudo Goldstone boson, the axion [17U18| . 

Let us first assume a new axial symmetry U(1)pq having a color anomaly N c . If the symmetry is 
spontaneously broken at a very high energy /pq (it has to be like this since it is not apparent at low 
energies) then the low energy effective lagrangian would include a massless Goldstone boson, £. Because 
of the anomaly, £ develops an interaction terrri 11 ! 

-^J^G%"G a , v . (1.14) 

In the rest of the low energy effective theory £ enters only through its derivatives as a Goldstone 
boson should. Now, provided that the rest of the theory preserves CPi 12 l except the #-term (II. 9|) then the 



effective potential is even in 9 + AA c £//pq, having a minimum at 9 + N c £,/fpQ = which conserves CP. 
The consequence is that the £ field develops an expectation value (£) such that 9 is canceled. We can say 
that the Peccei and Quinn mechanism relaxes 9 dynamically to zero solving the strong CP problem. 

Note that the PQ symmetry is almost completely similar to the QCD U(1)a- The subtle point that 
makes the difference is that now /pq is a scale much higher than Aqcd and therefore the instanton effect 
can be considered now a small explicit breaking of Upq(1) while it was a 0(1) breaking of Ua(1)- It is 
because of this reason that the axion exist and the GB of the Ua(X) does not. Due to this small explict 
breaking, however, the axion acquires a small mass, and hence it is rather a pseudo-Goldstone boson. 

Interestingly enough, many of the relevant properties of the axion are independent of the exact 

realization of the PQ symmetry. For example the £ field mixes with the neutral pion 7r° resulting in two 

massive eigenstatea^j, one very close to ir° and the other called axion whose mass can be easily obtained 

to be 

. , c y/m u m d F-k W 9 GeV 

m a = m^oN — — ~ 1.3 meV N — . (1.15) 

m u + m d /pq /pq 

Note also that being J7(1)pq an axial symmetry, the axion must be a pseudoscalar spin-zero boson. 

Also, the axion field always appears in the lagrangian suppressed by the Peccei-Quinn scale, /pq. The 
original Peccei and Quinn model set /pq close to the electroweak scale but such a low energy scale was 
soon excluded. The so-called invisible models [191 HO] propose a much higher value of /pq which makes 
the axion quite difficult (though not impossible) to detect. Remarkably, the Peccei-Quinn solution to the 
strong CP problem works independently of the size of this parameter. 



11 We can read it in eq. (|1.2p where we see how Goldstone bosons couple their relative currents. This mechanism is the 
same that provides ty° an interaction j^j-F^F^,; responsible for the anomalous ty° — + 77 decay. For further details see 

CP violation coming from the electroweak sector apart from the quark mass matrix is just a small perturbation to this 
situation. 

13 Corrections due to the other 0~ mesons are small and irrelevant for this discussion. 
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Figure 1.2: Feynman diagrams for the axion coupling to two photons through its components along the original 
7T° and £. Here, / stands for a sum over fermions carrying Peccei-Quinn charge. 

Furthermore, it is quite generic that the axion develops a two photon coupling, denoted as 

- 9 -fa F^F,,,, (1.16) 

with F^ v the field strength tensor of the electromagnetic field. The coupling g ai is the sum of two 
contributions, one coming from a plausible C/(1)pq electromagnetic anomaly and the other coming from 
a — 7T° mixing (See Fig. II. 2[) 

9a 7 = , A/ ■ (1.17) 

27t/pq V 3 m u + m d J 

Where a is the fine structure constant and M e is the electromagnetic anomaly of l7(1)pq. A cancelation 
of this coupling is rather unlikely due again to the different nature of the two contributiona 14 !. The axion 
has also couplings to quarks and leptons which have been extensively considered in the literature but, as 
it will become clear, it is the electromagnetic coupling which has the richest implications nowadays. The 
phenomenology of this coupling will be the main subject of Chapter [3] 

There are beautiful examples of Goldstone bosons other than pions and axions. For instance the so- 
called Majorons |22l [23] come from the hypothetical spontaneous breaking of the "Baryon minus Lepton" 
number U(1)b-l which was invoked to provide a Majorana mass for neutrinos in a Higgs mechanism 
fashion. Other celebrated examples are Familons [24], arising from the spontaneous breaking of family 
symmetries, which are very much related to axions and the strong CP problem. 

The general lesson to be learned from this section is that Goldstone bosons arise whenever there is 
a global continuous symmetry which is spontaneously broken by the vacuum state. If the symmetry 
is slightly explicitly broken the corresponding Goldstone excitations become massive. Interactions of 
Goldstone bosons are always suppressed by the energy scale of the spontaneous breaking and involve 
derivatives except in the case of the anomalous interactions (|1.14p coming from the axial anomalies of 
their corresponding symmetries. The parity and internal quantum numbers of the Goldstone are those 
of the charge density J° associated to the broken symmetry. 

Quite naturally these Goldstone Bosons can acquire a coupling to two photons through quantum 
corrections like those giving the anomalies (See Fig. II. 2\) . But there is one important point to remark 
here. In what follows I will be considering the phenomenological implications of vertices like (|1.16j) not 
restricting myself to the axion case. It happens that they are essentially similar (See Chapter [3]) to those 
coming from the corresponding parity-even gauge field structure 
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F^F^ , (1.18) 



which I will consider altogether. There are however important theoretical differences between them. 
F^F^ is CP and P conserving (See section I2.3.ip and thus (|1.18|) respects these symmetries for a CP 



An accidental tuning however is possible in some concrete models. See [21] as a notorious example. 
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and P even field <f>. This Goldstone mode is to be originated from the spontaneous breaking of a vector-like 
symmetry (rather than axial). We know that this is not possible in theories with only vector- like fermions 
(with non-chiral interactions) like QED or QCD because of the Vafa-Witten theorem [25]. However, this 
is certainly not a strong restriction; the standard model itself is not subject to it, indeed. 

However f 1 1 . 1 8 [) cannot have the same origin than (J1.16P since there is nothing like vector anomalies. 
Put it another way, F^ U F^ V is not a total derivative. As Goldstone bosons have derivative interactions, 
(|1.18p cannot arise, even for a parity-even Goldstone boson. 

These theoretical considerations disfavor somehow parity-even massless scalars coupled to two pho- 
tons. However (as we saw in the case of the axion) typically global symmetries are not exact and the 
full effective description of particle physics at low energies will contain terms that explicitly break the 
symmetry. We know that these terms would provide a Goldstone boson mass so in principle they can also 
naturally lead to interactions like (|1.18|) . 
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1.3 Paraphotons and millicharged particles 

Extensions of the standard model quite generally contain extra U{\) gauge factors. These gauge factors 
can be fundamental or sit within a non-abelian gauge group. Such gauge symmetries are necessarily 
accompanied by their corresponding spin-1 massless gauge boson which is automatically a new low energy 
candidate for new physics. 

Gauge bosons can obtain mass from the spontaneous breaking of its related symmetry but purely U(l) 
symmetries allow also for a Stuckelberg mass-term [51 (6] [7] [8] . At the same time there is at first glance 
no reason for their masses to be small, on the contrary, it is more natural if they are of the order of the 
highest mass scale in the theory, which it is usually assumed to be the Planck mass m-p\ ~ 10 19 GeV. Of 
course, these arguments do not rule out a small mass-scale, we know at least two completely "unnatural" 
mass parameters in nature, the electroweak scale (~ 100 GeV), and the cosmological constant (~ 0.1 
meV). 

There are many examples of such particles in the literature: mirror photons, Z's, shadow photons, 
U-bosons, etc. They differ in the interpretation that is given to the new U(l) symmetry or the way it 
links with the standard model sector. Here I will refer mainly to "paraphotons" , a term first introduced 
by Okun [26]. To my knowledge he was the first one to consider these particles, at least in the framework 
I will discuss here. 

On a basic level, the impact of paraphoton models is to proportionate deviations from the predictions 
of QED, at least at low energies. But also, paraphoton models include naturally millicharged particles, 
i.e. particles whose electrical charge is not a multiple of the electron charge but a small numbed 15 ! 
They are, together with paraphotons, necessary ingredients of our solution to the PVLAS-astrophysics 
inconsistency. Interestingly enough, if the millicharged particles are fermions, they can have naturally 
small masses in chiral theories like the standard model. 

Quantum electrodynamics is usually said to be the most precise theory of all sciences so deviations 
from their predictions are extremely restricted. This implies that paraphotons and particles coupling to 
them have to be somehow isolated from our standard particle physics. It has become a fixed tendency 
to consider that paraphotons and other light paracharged particles coupled to them inhabit a kind of 
hidden sector. The connection between this and our sector would be then driven by radiative corrections 
with very massive particles acting as mediators between the two worlds. 

Consider a toy model in which we have quantum electrodynamics^! in our sector (labeled 0) and a 
similar theory for the hidden sector (labeled 1). The complete gauge symmetry is then U(l)o x U(l)i 
and the gauge piece in the lagrangian is the sum of two similar parts 

£qed = £o = --F Qllv Fg v + e QoV'o4oV ; o (1-19) 

A = -jF llu ,F^ + eiQi^A^i 
plus a term that mixes kinetic terms of the gauge bosons 

£mix = —-^Fq^Fj 1 (1-20) 

where Fq" = d^A^ 1 — 8 u Aq ± are the field strength tensors of the gauge bosons Aq 1 , eo,i are the respective 
gauge couplings and Qq i are the 0, 1-charges of the different fermions ijjq \. We are interested in models 



See the discussion in [27] about charge quantization. 
16 For my purposes it is sufficient to consider how paraphoton models alter the standard QED scenario. A wider scope 



would consider the whole standard model 
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Figure 1.3: One loop contribution to e. Only multichargcd fcrmions ipo\, acting as mediators, run in the loop. 

in which there is isolation of the and 1 sectors so I put the tree level term (I1.20P to zero at zeroth 
ordeio and consider that the masses of mediators, fermions charged under both £/(l)'s (V'oi's) are very 
high so they can be removed from the low energy physics. 

Although decoupled, these fermions (of mass rrij) lead to a small correction 18 ! to e through the diagram 
depicted in Fig. 11.31 

^S^E^o 01 ^ 01 ^^! • (1-21) 

The dimensionful logarithm is understood to be regularized in any of the standard ways but a simple 
choice of the charges of the heavy fermions can evade this complication. For instance, I can choose two 
fermions, t/jgi ^ (1, 1) and V'oi £ ( — 1,1), to get e = ^^Log ^-. If one does not assume a fine cancelation 
due to mass degeneracy and use e\ ~ eo, the obtained value is of order 10~ 3 . However one expects 
indeed higher degrees of symmetry and degeneracy as the energy scales involved (here the f/'oi masses) 
are increased and thus it is "natural" to obtain much smaller values for e for example in supersymmetric 
models or in string theory [29j . 

We call {^cb^i} the "interaction basis" of the gauge bosons because the interactions with their 
relative fermions are diagonal by definition. In opposition to this we can find a "propagation basis" 
{^o^i} in which (|1.2U|) is rotated away and the gauge fields are eigenstates of propagation. In the 
course of diagonalizing the kinetic gauge sector we are led with mixed charge assignments. Then, in our 
calculations, mixing terms can be traded for off- diagonal charges, a fact that I reviewed in \30\ I3TT] in the 
context of the model proposed in [2] solving the PVLAS- Astrophysics puzzle (Chapter [6|). This can be 
shown in two ways, either by direct diagonalization of the kinetic part of the lagrangian or by studying 
the ipo^i elastic scattering amplitude. Let me show both of them. 

1.3.1 Diagonalization of the kinetic matrix 

First write the kinetic part of the lagrangian in a matrix form, adding a general mass term that will 
prove to be very convenient later on 

Ckinetic = --F T M F F + -A T M A A . (1.22) 

Where F, A are vectors with 0, 1 components and Lorentz indices have been suppressed. I choose the 
mass matrix to be diagonal in the interaction basis since it is more likely that this mass is gained through 
interactions, then 

1 e \ , , ( m\ 



M "=[ e i ; M * = ( o m ? ■ < 123 > 



17 There are two natural ways to obtain this: let U(l)i sit within a non-abelian gauge group and the tracelessness of its 
generators will ensure e = 0, alternatively consider a discrete symmetry Aq «-> —A\ at high energy broken for instance by 
fermion masses at lower scales. Both cases have been discussed in [29] . 

18 Other radiative corrections produce the usual charge, mass and field strength renormalization which do not interfere 
with my conclusions here. 
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To diagonalize Mf and express the result in canonical form we perform a rotation R of 45° followed by 
a contraction C ([CR]Mf[CR] t = 1)- In the resulting basis, Ma is not diagonal but it is symmetric so 
a further rotation R' is able to diagonalize it also ([R'CR]Ma[R'CR] = DiagjmQ, mf}). The matrix 
transforming from the "interaction" to the "propagation" basis is therefore U = R'CR and takes a very 
simple form when we consider the interesting asymptotic case e <C 5 = |mg — m\\/{m^ + m\) 

■1 -2 



U(e«5)= ml « 

mf—uiQ 

The fermion interaction in this basis acquires the promised non-diagonal charges. For instance, we get a 
e-sized electric charge for the paracharged fermion 

2 

eiQiffi^iVi-freiQi^Odfi + e ™° 2 io)^i , (1-25) 

m| — TTIq 

and a small paracharge to standard model charged particles, for instance electrons t/Jq, with Qq = — 1, 
implying 

2 

e Q r4o^ e -> e Qo^(io + € Z" 1 2 ^)^ e • (1-26) 

ttiq — mf 

The massless case mo = TOi = is special. The freedom to choose the matrix R' becomes an arbi- 
trariness of the definition of the propagation basis and thus of the charge assignments. It is conventional 
then to define the propagation photon Aq along the direction of the interacting photon so standard model 
charged particles, V'o'S) do not interact with the resulting paraphoton A\ but parafermions ipi acquire an 
e-sized electric charge. In this case we have 

f/(m = mi = 0) = ( ^ °A . (1.27) 

However, this is nothing but a definition without physical meaning. We could also have chosen a situation 
in which standard model particles have paracharges, or even a mixed situation. I have found that one 
can get the most general result from (jl.24p performing the limit mo = rm\ — > ending up with 

C7(m ,i =0)=(\_ "Y 1 ) ( L28 ) 

for arbitrary r as long as r ^ 1 to ensure that the off-diagonal elements of (|1.24p are always small. 

It is interesting to note that one cannot get this result by putting tuq = mi into the corresponding 
U(e ^> 5) and letting iuq — ► . Although mo = mi seem to imply also a rotational freedom, a small 
e-sized mass difference arises in the propagation basis mo 7^ mi breaking the rotational invariance (This 
mass difference is proportional to mo = mi so this does not spoil the freedom in the massless case.). 

In general So, ei, mo, mi, etc. will differ from the interaction values in small e-sized quantities. 

1.3.2 /0/1 scattering amplitude 

The charge assignments of the last subsection can also be obtained in a very simple calculation, a 
procedure that turns out to be much simpler in models involving several paraphoton fields (I will need 
two of them for my purposes). Consider the elastic scattering amplitude of a 0-fermion with a 1— fermion, 
by using the mixing term in (|1.20|) as an interaction (It can be done provided e is small) 

A = (eoQoJ°)^^(^ 2 )^^(eiQij^) (1.29) 

p q z — wiq q z — mf 
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Figure 1.4: Mixing between massive photons is equivalent to non-diagonal e-charges for the fermions. 



(q is the 4- momentum carried by the bosons, jn,jn the currents of 0,1 particles and I have used 
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(1.30) 



This amplitude is equivalent to the sum of two single boson exchange diagrams as shown in Fig |1.4l 
which require that we assign an e-charge to ip\ and a e-paracharge to ipo: 



QT = eQ\ 
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Q-mfei 



(1.31) 



We naturally associate these intermediate bosons with the propagation states we calculated before, and 
we find (always at first order in the small quantity e) from (11.19j) . (|1.24j) and (|1.3ip that, as expected, 
both approaches give identical charge assignments. 

The diagram studied, however, does not give us information about the values of the order— e cor- 
rections to mo, mi, eo, e\. To calculate them we shall study e~e~ or f~f~ elastic scattering. In that 
case, beyond the th order diagram in which the two electrons (/s) interchange a 70(71) there is another 
contribution in which the 70(71) intermediate state turns for a while into a 71(70) to reconvert again into 
the original "flavor". These diagrams provide the o(e 2 ) corrections to the masses and gauge couplings. 

If mi ~ mo the expansion (11.3QH would break down and this method probes to be not very useful. 
However, again the massless case is special. Use again the limit mo = rm\ — ► with \r — 1| > e in (|1.30p 
or (jl.3ip to get the same assignments (jl.28p we found within the previous diagonalization procedure. 



1.3.3 Some remarks 

First note that only particles charged under U(l)o will receive a e-sized charge of type 1 and viceversa. 
Neutral particles as the neutrino will continue being completely neutral if photons mix with a new gauge 
boson from a given sector. Notably this holds for composite "states" as for instance atoms. Also, notice 
from (jl.26p that for instance, electrons will gain an opposite paracharge than protons because their charge 
is opposite. 

Second, note that the matrix (|1.24p is not unitary, reflecting the fact that "interaction" states are 
not orthogonal, since they have a kinetic mixing term. Such a complication does not make difficult the 
calculation of the 70 — 71 oscillation probability. If we prepare (70,1) states using (jl.24p their temporal 



1.3 Paraphotons and millicharged particles 



13 




Figure 1.5: The two basis of the /7(0)o X E/(l)i theory. The interaction basis is labeled 70,71 (and displayed in 
black) and the propagation basis 70,71 (displayed in red). The first one is non orthogonal because of the small 
mixing of size — e. The two angles are assumed to be small (but magnified in the picture) and they also differ in 
the small quantity e. The sterile state j s with respect to the U(l)o interactions of the standard model is displayed 
in blue. 



evolution will be easily expressed in the propagation basis 
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(1.32) 
(1.33) 



where E{ is the energy of the state |7j) and I have not taken into account the normalization, irrelevant 
for the following computation since it just provides 0(e 3 ) terms. We find 



-P(7o -> 7i) 



l(7i|7o(t))P 



1 + ; o o,„ sin 2 
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AE 



t))+ 0(e' 



(1.34) 

(1.35) 



Where AE = Eq — E\ = (ttiq — m\)/2E for the relativistic limit E » mo,i- We see clearly that in 
addition to the typical oscillation sinus, a t = contribution has arisen, due to the fact that 70 has a 
small component along 71 (another way of looking at the mixing term) of size — e. This is made clear in 
Fig. O 

Another interesting difference with neutrino oscillations is that in the case mo = 0, there are no 
oscillations, even if ttoq = 7^ m\. This situation corresponds to one where Wgoo = ' and then 71 = 71 so 
clearly no oscillation could happen since one of the "interaction" states is a "propagating" one, namely 
<7i|7o)W = -e(7il7i>e i£;ii . 

It is common to consider experiments of disappearance of photons (70), as for instance light-shinning- 
through- walls experiments (See subsection 13.3.31) . It will prove useful then to define the sterile states 
7s0, Tsi, as the states orthogonal to the interaction ones 70, 7i- The probability of 70 — 7^0 oscillations is 
therefore derived in complete analogy to the neutrino oscillations case to be 



-P(7o^7so) = |<7sol7o)(*)| : 
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(1.36) 
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These sterile states were called " paraphotons" by Okun [26] although here we have used the name for 
a general abelian gauge boson. It is mandatory to compare this construction with the original work of 
Okun [26J. There are two main differences. First, he did not consider a current of type 1. Second, he did 
not wonder where do the mass terms come from. In my work I need particles charged under the hidden 
sector gauge group U(l)i, but most importantly, I assume that the mass terms come from interactions 
within the two gauge groups U(l)o,U(l)i separately, i.e. the interaction states have a diagonal mass 
matrix. 



Chapter 2 

The PVLAS experiment 



2.1 Introduction 

A photon can convert into a neutral boson eft in an electromagnetic field in the presence of a ^77 interac- 
tion. This is the so-called Primakoff conversion first proposed by Henry Primakoff to measure precisel}o 
the 7T° — ► 77 decay constant [32] . 

Moreover, if the <j> particle is light enough, this vertex would lead to 7 — eft "mixing" where a coherent 
superposition of the two particles arises, in complete analogy to the iiT°-meson or neutrino systems. 
Notably, in this case and due to the presence of the magnetic field, the angular momentum needs not to 
be conserved, allowing the mixing of photons with particles of higher (but integer) spin. 

This immediately suggests a possibility for looking for novel <fi bosons coupled to light in a classical 
optical setup. Consider Fig. I2.1j q. a photon from a laser beam interacts with an external magnetic field 
and it is converted into a <f> boson. If this absorption is selective, i.e. is different for photon polarizations 
parallel (||) and perpendicular (_L) to the magnetic field direction, the vacuum permeated by magnetic 
fields would act exactly as a dichroic medium since both polarizations will be depleted at different rates. 
The net effect on a polarized beam will be a rotation of the polarization plane which can be measured. 
This can be seen in l2.11 6. where I depicted a view of the plane perpendicular to the laser direction showing 
in and out polarizations of the incoming laser. It is straightforward to see that the sign of the rotation 
will be positive if the coupling favors 711 — (j) transitions and negative in the opposite case. As we will see, 
the sign of the rotation will be used to determine the parit>o of (/). 

Maiani, Petronzio and Zavattini were the first to realize that the two photon coupling of low mass 
bosons (in fact they considered only spin zero bosons) also leads to birefringence for light propagating 
in strong magnetic fields [33J . The <fi boson can be produced and, after some distance, reconverted_| as 
shown in Fig 12.11 c resulting in a phase delay of the photon. Again, if 711 — <p — 7| transitions are more 
likely than 7_i_ — <^> — 7_i_ or viceversa there would appear a phase difference between the two polarizations 
and then the outgoing beam would show a counter-clockwise or clockwise ellipticity (See Fig. I2.1l rf). We 
can conclude then that the vacuum permeated by magnetic fields would behave as if it were a dichroic 
and birefringent medium. 

Indeed it was known long ago (See [33] and references therein) that pure QED effects already produce 
birefringence in a transverse magnetic field in vacuum (a vacuum Cotton- Mouton effect). This is due to 
virtual electron loops like the one depicted in Fig. 12.21 which lead to the following indices of refraction 



In this case they use the Coulomb potential generated by an atom as the electromagnetic field. 
2 More precisely the parity structure of the two-photon part of the ^77 interaction. 
3 This can happen even off-shell, i.e. with 4-momentum squared different to m|. 
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Figure 2.1: A laser beam traveling through a transverse magnetic field develops a small rotation of the polarization 
plane and a small ellipticity in the presence of a new boson cf> coupled to two photons if this coupling is different 
for both photon polarizations. The real production of (f> depicted in (a) accounts for the dichroism, as shown in 
(b) where we see the plane transverse to the propagation and the two laser polarizations are labeled Eu and E± 
(Electric field of the laser beam parallel and perpendicular to the external magnetic field). The </>'s reconverted 
into photons (c) are phase-shifted and account for the ellipticity shown in (d). 
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Figure 2.2: QED Feynman diagram for the vacuum birefringence of light in a magnetic field. There are other 
diagrams differing from this in permutations of all photon lines. 
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It was T. Erber [35] the first to propose the use of a polarized laser beam propagating in a magnetic field 
to measure the vacuum QED effects. However, he was wrong thinking that it leads to a rotation of the 
polarization plane, not to an ellipticity. Some years after Iacopini and Zavattini [36] wrote an interesting 
experimental proposals to detect this ellipticity at CERN. Time passed by and the paper of Maiani, 
Petronzio and Zavattini [33] came out supporting the idea of such an experimental with the possibility of 
detecting low mass bosons coupled to light (at a time in which axions and other light bosons like arions 
[38], majorons [22] or familons [23] were very fashionable). 

Finally, a pioneer experiment was performed in the Brookhaven laboratory by the so-calleqj BRFT 
collaboration. They found no signal of this kind of 0^77 interactions^ [391 I40j so they were able to establish 



Apparently unaware of Erber's idea. 

5 Little time after, Raffelt and Stodolsky [37] showed that this paper contained some overestimated predictions. They had 
claimed that they would be able to measure the QED birefringence effect and values of g al ~ 1CP 10 GeV~\ competitive 
with the strong astrophysical bounds. 

6 From Brookhaven, Rochester, Fermilab and Trieste, homeplaces of the members of the collaboration. 

7 This experiment, as we will see, was also sensitive to the existence of paraphotons. 
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Figure 2.3: Brief scheme of the PVLAS experiment essentials. A magnetic field between two crossed polarizers. 

exclusion bounds on the ALP parameter space. Most of the collaboration moved then to the National 
Laboratory of Legnaro, near Padova (Italy), where, under the leadership of E. Zavattini, they started to 
build the PVLAS experiment [H] as a natural continuation of the BRFT apparatus and technique. Next 
I will comment on its set up to end up with their latest results. 



2.2 PVLAS experimental setup and recent results 

The PVLAS experiment is as simple as beautiful. The basic idea consists in sending a laser beam along 
a strong transverse magnetic -B e xt field located between two crossed polarizers, all in a cavity at high 
vacuum, see Fig. 12.31 It is clear that if light does not interact with B ext , then no light would came 
out from the second polarizer. Conversely if the magnetic field in vacuum has optical properties such as 
dichroism or birefringence then the situation would change and some emerging light would "enlighten" 
the phenomenon. 

As I mentioned this should be the case if one takes into account quantum fluctuations of the QED 
vacuum structure, producing a net birefringence but no dichroism. Also, in the case that a light boson 
with a two-photon coupling exists, the real and virtual production shown in Fig. 12.11 will lead to both a 
rotation and an ellipticity of the outgoing beam. 

In the PVLAS experiment the magnetic field intensity reaches 5.5 T and the length of the magnet is 
1 meter. With these parameters, the two perpendicular components of light emerging from the magnetic 
field would have accumulated a QED-induced phase shift ~ 10~ 16 from (|2.ip implying an ellipticit)o of 
~ 10~ 16 which means that only a ridiculous amount of light will traverse the two polarizers. Even using 
a MegaWatt laser (which probably would destroy the polarizers), the emerging light would be far beyond 
the reach of any present (and maybe future) photodetector. 

For this reason, the PVLAS signal is enhanced in two different ways. First of all, the magnetic field is 
located inside a pair of highly reflective mirrors (actually a Fabry-Perot cavity) which increase the path 
of light through the magnet by a huge factor N ~ 44000. This might look quite a cheap way of building 
a 44 km length magnet, but it is not exactly the case as I will show below. Second, the outgoing signal 
is enhanced by heterodyne combination with an ellipticity modulator (called SOM for Stress Optical 
Modulator) before entering the detector. The principle of heterodyne detection is very simple, two waves 
at frequency lo\ and 002 traveling through a non linear material will produce outgoing waves at frequencies 
uj\ ±u>2- The amplitude of these waves depends linearly on both amplitudes of the incoming waves. This 
trick can be used to amplify a very small amplitude by combining it with a larger one and look at any 



* Ratio of the minor to the major axis of the polarization ellipse shown in l2.1l d. See section [2. 3. 3 1 
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Figure 2.4: A closer look to the PVLAS experiment. The incoming light is linearly polarized and travels several 
times the interior of the Fabry-Perot cavity where a strong magnetic field is placed. The magnetic field is transverse 
to the propagation of the light and slowly rotates its direction at an angular frequency loq. At the exit of the Fabry- 
Perot the light could have developed a small rotation and a small ellipticity due to mixing with ALPs or pure QED 
effects. A quarter wave plate can be used to convert rotations into ellipticities and viceversa. The SOM superposes 
a time dependent ellipticity to the incoming light. The small ellipticity (beating at 2o>o) acquired at the magnet 
will beat with the SOM time-modulated ellipticity allowing heterodyne detection after the second crossed polarizer. 
The rotation entering the SOM will be left unchanged and because its smallness will pass undetected. 



of the interference waves. This procedure, however, demands that the small incoming wave to detect 
(with our desired excess ellipticity) is beating at some frequencjo and in order to do so, the PVLAS 
collaboration managed to make the magnet rotating at a small frequencM 10 ! ~ 0.5 Hz. 

The whole resulting setup is depicted in Fig. 12.41 and is able to measure ellipticities as small as 
10~ 7 rad Hz - 2. This setup can be adapted very easily to measure rotations because a properly oriented 
quarter wave plate transforms small rotations into ellipticities and viceversa [42] . so if placed just before 
the SOM it will convert a possible rotation generated by the magnetic field into an ellipticity which will 
be again enhanced by heterodyne combination in the SOM and thus detected. Also, interchanging the 
slow and fast axis of the quarter wave plate produces a sign flip in the outgoing ellipticity, providing a 
check of the nature of the signal detected. 

With this setup the PVLAS collaboration has recently reported [53] an excess rotation of 3.9 ± 



0.5 10 12 rad/pass for laser light of 1064 nm traveling through aim length, 5 T, magnetic field. This 
signal has been reported after careful search of possible systematic effecta 11 ! over several years. Also they 



have reported in several conferences |45[ |46| [44] positive results for the search of birefringence, again 
at the 10~ 12 rad/pass level. This amazing result exceeds by three/four orders of magnitude the QED 
expectations and nowadays remains without a conclusive explanation within conventional physics. 

However, there are two non standard candidates that in principle could account for the PVLAS 
results within particle physics. The first one is a neutral boson weakly coupled to two photons, one of 
the original motivations of PVLAS. Such a particle will be called from here on axion-like particle (ALP). 
I devote next section to examine this possibility. 

The second possibility has arisen very recently [27] and assumes the existence of very light millicharged 
particles (MCP). Again, the reactions distinguish both photon polarizations and will lead to dichroism 
and birefringence. 



9 The original setup in [33] used a modulated Faraday cage. PVLAS is nowadays involved in the modifications required 
to implement this other technique as a crosscheck of their results. 

10 The angle twisted in the photons-time-of-flight is a ridiculous 10~ 3 rad. 
11 For a review see [44] 
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A very interesting paper has recently appeared comparing both ALP and MCP scenarios and it seems 
that we will have to wait until more experimental data is released to decide which of them, if any, is more 
likely to be the responsible for the PVLAS signal |47j . 

Both explanations have many things in common however, the first one being that they are not 
satisfactory at all because they are excluded by astrophysical arguments, as I will show in Chapter [3j 

But moreover, they share a much more interesting feature, both scenarios can be made compatible with 
astrophysical observations within a concrete model I presented in [2] and constitutes the central work of 
this thesis. 
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2.3 Axionlike particle interpretation 

The study of coherent 7 — 4> oscillations in presence of a long, static homogeneous magnetic field has been 
addressed by many authors [33l E3 HS1 H9] . Here I follow the exposition of Raffelt and Stodolsky [50] , 

I do not feel the necessity of including the effects of the QED vacuum since they are very small 
compared with the PVLAS signal. But it is worth to recall that any birefringence will show up into 
a rotation when the magnetic field rotates (as in the PVLAS setup) [5U E2]. Nevertheless, it is quite 
intuitive (and has also been proved rigorously in [52]) that these effects are negligible for the PVLAS 
setup 12 l. In what follows I neglect the rotation of the magnet. 

Consider the following situation. A beam of laser light is shone against a region where there is a 
strong magnetic field. I choose the beam direction as the z-axis and the magnetic field polarization to 
lie in the xz-plane. It is convenient also to make the problem one-dimensional, in the sense that all fields 
we want to determine are functions of time and z only. 

A»{x u )=A»{t,z) ; <t>(x") = <p(t, z) (2.3) 

This is formally the case if the magnetic field would be infinite in the x — y directions and the laser is 
described as a plane wave. Effects from deviations of this simplistic picture are small. 

Next, I consider that there is an interaction between two-photons and a new neutral boson, called <f>, 
which is very light (we will see soon how light this "light" means). The most discussed case is that <j> 
represents a spin-zero particle, but also one might consider that 4> is a spin-1 particle like a paraphoton 
or a spin-2 particle like the graviton. 

In Chapter [1] we have found motivations for two different couplings associated with opposite parity 
structure (II. 16ft and (j!.18[) which I rewrite here for completeness 

\g + F^F, u <p ; \g-F^F^<l> (2.4) 

We know that parity is not exactly conserved by the weak interactions so there is no strong reason a 
priori to think that it should be conserved by new additional interactions. This suggests not to either 
assign a definite parity to <f) nor to remove one of these two interactions. However, if ^ is a true Goldstone 
boson then it could easily happen that it has a definite parity. Explicit breaking terms can modify this 
statement, but corrections should be small and then one of the g's will be much smaller than the other. 
Moreover, unless g + and g~ are fine-tuned 13 ! these considerations will not change our main conclusions. 
This holds also for the case in which there is more than one neutral boson <j>. 

It is then enough to consider only one of the two interactions. Following criteria of theoretical and 
historical importance I choose the parity-odd interaction which arises for axions and other celebrated 
pseudo Goldstone bosons and drop the — sign of the coupling for simplicity. 

The lagrangian for the <fi — 7 system i a 14 l 

C = —F^F^ + J^A* + \d^d^ - imV + \gF^F^ (2.5) 



whose equations of motion are 



dpF"" = r + gF^drf (2.6) 

(d 2 + mm = -\gF, v F^ (2.7) 



j_> 



The attempt made in [51] to explain the PVLAS signal in these terms is mistaken. 



Maximal parity violation in a hidden sector like this might not propagate to the SM sector and thus it should not be 
excluded at first glance. 

14 In all this thesis I am using the metric r]^ — Diag{l, — 1, — 1, —1}. 
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where I have included the electromagnetic current J u . We can divide the field strength into two parts 
corresponding to the external magnetic field (with the definition 15 ! d^F^ t = J u ) and the laser beam, 



F^ = F e Z + Fj^ . 



(2.1 



If we consider that the external magnetic field is much more intense that the laser beam we can ne- 
glect terms involving 
differential equations 



gleet terms involving two fields, c/>F~/ and FJ^ v F-y^ u to find a linear system of two coupled second order 



vr 



(d 2 + m$) 






(2.9) 
(2.10) 



Note that this is the case in the PVLAS experiment, since | -B cx t| ~ 5T ~ 1000 eV and the laser power is 
~ 1 mWatt in a 2 mm diameter gaussian beam which means \E~\ = \B~\ ~ 0.1 eV (taking into account 
the Fabry-Perot cavity which enhances \B\ by a factor ~ 10 4 because the coherence of the reflections.) 
Then the wave equation for the time-varying potential A (in Coulomb gauge V • A = 0) is 



d 2 A = gB T d t 



(2.11) 



(d 2 + m 2 



-gB T -d t A 



Where Bt is the transverse part of the external magnetic field (Bt • z = 0). As the coefficients of the 
unknown fields A, <f> are time independent it is natural to guess a solution of the form 



which implies 





A(t,z) = 




<t>(t,z) = 


(lu 2 + d 2 z )A ± (z) = 


^ 2 +(. 


uj 2 ig 
-igBxto to 2 



e iujt A(z) 
e iwt 6(z) 



igB^u) 



m~ 



A\ 



(z)=0 



(2.12) 
(2.13) 



(2.14) 

(2.15) 



Recall that the electric field will be E = dtA = iuoA. For simplicity in further expressions I will evade 
imaginary quantities redefining A — > i A . 

The solution of the first equation is trivially a plane wave with momentum k± = uj. In the coupled 
set we can made a further approximation if the photons and (j) are very relativistic, in this case uj 2 + d 2 = 
(uj + id z )(u) — id z ) ~ 2u)(ui — id z ). The resulting equation is therefore of first order and easier to solve 



The solution reads 



UJ 



-4, 



id z + 





gB T 

2 



gB T 

2 

9 

Tlf. 


2uj 



A 



(z) = Exp (—iuiz) Exp (— iMz) 



Ai 



(0) 



(2.16) 



(2.17) 



15 This is not a "local" statement. A far away J v generates the magnetic field at the region where the laser is shinning by 
contour conditions. 
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being M the 2x2 mixing matrix in (|2.16p . To get a useful expression for An (z) we must diagonalize M 
with a rotation R of an angle 9 defined by 

/eos* -sin0\ 1 tan2d = -^ (2.18) 

y suit' cos# J 2 m. 

The resulting states (A±, A2) = (j4ii, (/>)-R t are eigenstates of propagation because for them the evolution 
is decoupled 

^ ) (z) = Exp (-iuz) Exp (-iDz) (^) (0) (2.19) 



-^0 + V / ^+(23^tc) 2 n ^ 

with V = I 4w 



(2.20) 



n ^_v £ / 

U iui ) 



These solutions carry momentum p\ = w+Poooo andp2 = oj +£>£<= and invariant masses m\ = to 2 —p 2 and 
m\ = to 2 — p\. When we specify the initial conditions for 

A\\(z = 0) = l,4>(z = 0) = we find 

A\\(z)= [R T Exp(iDz)R] u = e~ lulz (e~ lVllZ cos 2 9 + e - iV22Z sm 2 9) (2.21) 



z 



[i? T Exp (Wz) R] 2 i = e~ luJZ (e- Wllz - e'™™*) sin9 cos 9 (2.22) 



In order to get formulas for the probability of appearance of <j) particles and disappearance of photons, 
we interpret these Au(t, z),(j>(t, z) fields as probability wave functions. Then we find the probability of 
7| 1 — (j) oscillations at a distance z and time t 

P(7||-$= |<K*)M N (0)| 2 = sin 2 2fl sin 2 Qa osc zJ (2.23) 

P( 7 ||-7||)= |A N (z)M N (0)| 2 = l-Pfy|-0) (2.24) 

with 

\H + (2ffi?Tw) 2 
A osc = Pu - P 22 = -* — (2.25) 

The typical oscillation pattern arising from these formulae is shown in Fig. 12.51 In order to maximize 
the signal there is an optimum distance at which to place the detector at z = mr/A osc , with n an integer. 

There are two very interesting limits of these expressions, the weak and the maximal mixing limits. 

Weak mixing limit, 9 ~ .- 
In this case tan 29 ~ 29 ~ —2gBT0j/'m\ ^C 1 and the corresponding eigenstates of propagation are very 
close to the original A\\,<j> 

A 1 = A\\-?-^-4>~A\\ (2.26) 

11 m z II 

A 2 = <j> + ?^A\\~<I> (2.27) 
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Figure 2.5: Oscillation probability of photons into <j> bosons as function of the distance z. 

1 

We also can expand the square root of the diagonal elements of T> to find that at first order 

g 2 B*u 



pi=u + 



P2=UJ 



m' 



2u> 



2mi 



+ 0(. 



and m 



-5— + 0(. 



m- 



g z B^oj 
2ml 



-.2 d2 , ,2 



+ o(.) 



„2 u2 , ,-i 

A 2 2,9 B T UJ 

and m 2 = m 6 -\ ^ \- Of.) 

mi 



(2.28) 
(2.29) 



also very close to the original photon and 0. To get the results usually quoted as the "weak mixing" 
limit, one has to assume also thall 16 ! ujz tan 29 <C 1, then we obtain 



Au(z) 



,2d2,.,2 



l - t^L | W 



m\z 



m. 



4lo' 



+ i 



mlz 



2^! 



mhz 



sin- 



2a; 



+ ... 



P (7|| 



4g 2 ^u; 2 . 2 
4 sin 

ml 



m^z 



4w 



(2.30) 
(2.31) 



Interestingly enough this weak mixing limit applies in the case of the PVLAS ALP parameters since 
for B r ~5T,w = 1.2 eV, cT 1 ~ 10 6 GeV, m^ ~ 1 meV we find 



tan 26* ~ 10"" 



(2.32) 



Maximal mixing limit 9 ~ 7r/4 .- 
In this case tan~ x 29 ~ cos 29 = sin(vr/2 - 20) ~ vr/2 - 20 ~ -m\j2gB T uj ~ 0, implying 20 = vr/2 + 
m 2 ,/2gBTOJ. The corresponding eigenstates and eigenvalues are 



A, = -= 



Ao 



V2 



V2 



A\\ + 



mz 



AgBxto 

'2 



(4 



Ai 



mt 



AgB-ruj 



(il||+0) 



(2.33) 
(2.34) 



Pi 

P2 = W 



W H h0(.) 



2 
5^r 



+ o(0 



and 



and 



m-i 



-£lB t w + 0( 



m 2 = gB^io + 0(.) 



(2.35) 
(2.36) 
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A\\(z) 



2 gBz _. mj _ : _ 5 5z 



1 — 2sin i „ sin 

4 gBw 



(2.37) 



P ( 7 ||-0) = sm 2 (^) (2.38) 

Note that in the m^ — ► limit there is no imaginary part, as expected from the fact that it comes 
from the difference in the propagation speed of photons and 0's of the same energy, which vanishes for 
m^ — > 0. Still there is a real part that provides 711 — <j> transitions. This is an opposite situation to 
neutrino oscillations because here the oscillations are driven by the interaction ^77 and not by mass 
mixing. 

Regarding the formulae for P (711 — 4>) we see that the limits z — ► in maximal and weak mixing 
coincide. However, for long distances the maximal mixing formula (|2.38p can grow until P = 1 while the 
weak mixing result (|2.3ip has a small maximum value. These results will show their interest later on, 
when I consider the possibility of measuring 7 — <j) oscillations in a low pressure gas. There we will find 
that controlling the gas pressure we can effectively set m^ — ► and get this maximal mixing situation. 

2.3.1 Parity arguments and parity-even ALPs 

It is quite remarkable that the set of equations (|2.1ip decouples into two sub-sets, (|2.14p and (|2.15p . This 
can be understood with parity arguments as shown in [37] . Instead of using P x = — x let us define the 
parity transformation 17 ! P* as 

P*X||=X|| ; P*x ± = -x x (2.39) 

Where now I denote as xn the coordinates projected onto the plane containing the magnetic field and the 
direction of the laser beam (the xz plane) and x± the orthogonal direction (the y axis). Then, electric 
and magnetic fields transform according to 

P*E\\(xx,x\\)=E\\(-x±,x\\) ; P*E ± (x ± ,x\\) = -E ± (-x±,x\\) 

P*B ]] {x ± ,x ll ) = -B ll (-x ± ,x ll ) ; P*fli(sL,X||)=BL(-x_L,Z||) (2.40) 

where again 1 1 , _L refer to the vector components in the xz plane and along the perpendicular direction 
y. This can be easily seen if one considers how E and B change if we apply this P* transformation to a 
charge or a current distribution (See Fig. I2.6[) . Moreover the charge conjugation operation changes sign 
of both electric and magnetic fields CE = —E, CB = — B . The external magnetic field changes sign 
under both C and P* so is left invariant under the combined transformation CP*. Finally, also the ALP 
interactions (|2.4p can respect this symmetry. The pseudoscalar interaction can be written in terms of the 
electric and magnetic fields 

^gF^F^j = -gE ■ B<j> (2.41) 



which at the view of (|2.4U|) is CP* invariant when 4>(x) is CP*-odd, i.e. if (j)(x) is P-odd 18 !. Conversely, 
the scalar interaction is written 



lg F » v F»^=-\ 



gF^F^cf) = --g(E 2 - B 2 )<p (2.42) 



and will be CP* invariant for a P-even <j)(x) field. 

17 This is nothing but the conventional parity P, with a further rotation of 180° in the xz plane. 
18 The 4> field is neutral so Ccj> = 4> and scalar so invariant under the rotation added to P* 
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Figure 2.6: P* transformation properties for electric and magnetic fields contained in the x — z plane (||) and 
orthogonal to it (_L). The transformation acts on their relative sources (a positive charge for the electric field and 
a circular current for the magnetic field) which then determines the transformed fields. 



In any of these cases, the whole dynamics is invariant under CP* and thus the eigenstates of the 
system can be classified by the eigenvalues of CP*. Therefore, at the view of Fig. 12.71 we can conclude 
that: 

\\-plane waves are CP*-odd and .L-plane waves are CP*-even 19 \. In the presence of Parity conserving 
interactions they can only "mix" with particles of the same parity. 

Splitting (|2.42p and (|2.4ip into the external magnetic field and the laser contribution we get : 



I 



IgF^F^t 



gF^F^cf) ~ 9 B-B 



-gE ■ B ext (, 



-gEu ■ B ext c 



CXt'i 



gBn ■ B ext i, 



-gE ± ■ B ext c 



(2.43) 



(2.44) 



where for the last equality I have used B = z x E for a plane wave in vacuum. At this point we 
understand that starting with the parity-even structure (|1.18j) we would have ended up with exactly the 
same system of equations (|2.1ip interchanging the role of Au and A± states. Finally, (|2.2ip and (|2.22p 
hold as the solutions for A± and a parity-even-^) mixing with the coupling (|2.42[) . 



2.3.2 Angular momentum conservation 

It is also is interesting to note that only transverse magnetic fields are capable of inducing transitions 
between photons and spin particles. In order for these transitions to occur, the "mixing agent" (the 
external magnetic field) has to match the quantum numbers, in this case the spin component s z . 



This is true even if some small longitudinal component arises in a medium. 
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Figure 2.7: CP* transformation properties for plane waves. ||-waves are CP*-odd and _L-waves are CP*-even. 

Formally we can write the external magnetic field as proportional to the angular momentum operator 
B ex t oc J and consider the 7 — (ft transition amplitude. Using y± for the photon eigenstates of J z with 
eigenvalue ±1 we find that 

(<ft\ J t \ 7 ±) = ±(<ft\y±) = (2.45) 

by conservation of the total angular momentum. So a longitudinal magnetic field (B z ) cannot mediate 
such transitions, although allows transitions between 711 and "y± states, as occurs in the Faraday effect to 
be discussed later on. 

Finally, I have also considered together with E. Masso and C. Biggio the case in which (ft is an spin- 
2 particle [53] . In this case we have explicitly checked that a transverse field can mediate transitions 
between photons and s z = ±2, states while a longitudinal one will be responsible for the transitions to 

4 = ±1. 

2.3.3 Rotation and ellipticity 

A general light wave propagating along the z-direction can be parameterized by means of the Stokes 
parameters 50,1,2,3 defined by 



SO — \E X \ + \Ey\ 

S 2 = 2Re (E*E y ) 



Si — \E X \ 



\E,,, 



S 3 = 2Im(E*E y ) 



(2.46) 
(2.47) 



Where the electric field is E = —ujA. 

From them it is easy to derive the formulae for the angle of polarization £, the ellipticity angle ift, 
and the ellipticity e, of a generically polarized plane wave 



tan2C = |^ 
.Si 



■ 1 ; ^3 
smzyj = — 

^0 



e = arctan ift 



(2.48) 
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Figure 2.8: Generic elliptical polarization. £ is the angle between the mayor semiaxis and the x direction (Parallel 
to -Bext)- The cllipticity e is the ratio of B and A, i.e. the minor over the major semiaxis. Understanding e as a 
tangent, the corresponding angle is called the cllipticity angle ip. 



The physical interpretation of these quantities is made clear in Fig. 
We can parameterize a small change in amplitude and in phase as 



En ± (z)=E lu (z = 0)(l-r ] 



'l-L 



llf: 



1,-L 



(2.49) 



-up 



and let rj, (p be small real quantities, for which we can use 1 — n — i<p ~ (1 — r\)e~ 

Then, if the initial beam is linearly polarized at angle Co (See Fig. I2.1i a ), after passing the magnetic 
field it shows a rotation of the polarization plane 



20 



<5C 



sin2£o 



and a small ellipticity or ellipticity angle 



ip ~ 



sin2£ 



Vn _ Vj 



(2.50) 



(2.51) 



The sign of the ellipticity (or the ellipticity angle) tells us if we have clock-wise (sign e > 0) or anti- 
clockwise (sign e < 0) elliptical polarization, and it is also measurable. 

It is easy to get 77,, and tp., from equations (|2.30p and (|2.37p . If we had selected a parity-even ALP 
we would have got the same solutions but replacing || for _L. In this case, we see that the rotation (|2.50p 
and ellipticity (|2.5ip will change sign, as I have argued from Fig. 12.11 



2.3.4 Further comments on the ALP interpretation 

We have calculated the resulting rotation and ellipticity on an initially polarized laser beam that traverses 
a strong transverse magnetic field once. However, in the PVLAS setup, the magnet is located inside a 
Fabry-Perot cavity so that the laser beam performs N reflections before entering the detector. If the <fi 
field does not interact with the cavity mirrors, then at every reflection, the (f> particles produced scape 
the cavity and the <f> field is reset to zero. The final effect is analogous to shinning a laser beam through 



20 We see here how if the magnet is rotating "adiabatically" 
twice the frequency of the magnet's revolution 2u>o- 



(£0 = tuot), as in the PVLAS setup, the rotation beats with 



28 The PVLAS experiment 



N copies of the magnet separated by barriers that absorb the ^-particles. The resulting amplitude for 
the photon field is then 

A\\(z) = 4||(0) [1 - ?](z) - icp{z)] N ~ ^||(0) [1 - N?]{z) - iN(f)(z)} (2.52) 

In [33] this point was unnoticed and they got a much more optimistic 

A\\(0)[l-ri(Nz)-i<j>(Nz)] (2.53) 

which will hold in models in which the (j> field is reflected from the cavity mirrors. Incidentally I found 
such a model in [53], studying general frameworks in which the astrophysical bounds on these particles 
could be evaded. 

Note that we have started with an interaction g^F^F^ which will be most probably the first order of 
a complete effective lagrangian. With the mixing approach in [37] we have managed to obtain a solution 
that involves further powers of g 2 . As long as additional effective vertices are expected in the realistic 
lagrangian, we cannot trust in principle our results further than the g 2 level. However, as long as we 
can neglect the terms involving two or more fields in the equations of motion the further interactions to 
consider, will give the same oscillation pattern, i.e. the same functions (|2.23j) . (l2,24|) but with corrections 
oc g 4 to the mixing angle and A osc . The same will hold for corrections to the relativistic assumption 
m^ <C oj. In this case the expansion in g is to be understood in 6 and A osc and not on Aii(z), (j)(z). 

2.3.5 Gas measurements 

Up to know we have considered 7 — <f> oscillations in vacuum permeated by a strong magnetic field. 
The PVLAS collaboration has also performed measurements filling the magnetic field region with a low 
pressure gas that may complete and give more soundness to the ALP interpretation. It is mandatory 
here to include the theoretical tools to interpret these results. 

It is well known that a gas in a transverse magnetic field becomes birefringent. This is the so called 
Cotton-Mouton effect I mentioned before. The physical idea is simple to explain. The forward scattering 
of light in neutral matter does produce a phase delay in its propagation. This is accounted for with 
an index of refraction which in isotropic materials cannot depend on the light polarization. A gas at a 
low temperature can be considered such an isotropic medium because their spins are randomly oriented. 
However, when you provide an external magnetic field these neutral atoms tend to align their spins with 
the magnetic field and produce an anisotropy of the medium. The forward scattering is thus different for 
the two linear polarizations (|| and _l_) and then we expect a different index of refraction for En and E±_ 
waves, i.e. the medium becomes birefringent. 

The resulting indices of refraction are quadratic in the external transverse magnetic field, Bt, and 
are given by 

'"II 

where A is the light wavelength and C is the so called Cotton-Mouton constant, that depends linearly on 
the pressure for low pressures. 

Remarkably, if there is a longitudinal component of the magnetic field also 711 — 7_l transitions are ex- 
pected, a phenomenon called Faraday effect. Hence, this effect also produces a rotation of the polarization 
plane for initially polarized light. In this case, however, the effect is linear in B^. 

Both magneto-optical effects can be included in our mixing formalism by including the corresponding 
forward-scattering terms to (J2.15P getting 



(nil - n ± ) c .M. = CXB 2 (2.54) 



uj 2 + d 2 + 2w 2 





(2.55) 
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When only a transverse magnetic field is present our results get modified in a very interesting and simple 
way. Again the mixing matrix splits in two parts (||,_L). The rotation that diagonalizes the resulting 
mixing matrix in the (_L) sector becomes 

~ tan 29 = 2 , ^ , S -^ (2.56) 

2 2w 2 (n|| -1) -m 2 Am 2 v ; 

and the nontrivial gas results follow from the substitution rtv 2 , — > Am 2 = m? — 2u; 2 (ni — 1) in (|2.23|) . 



(I2.24|) and in fj2.25|) which gives an oscillation length 



(m 2 -2a; 2 (n||-l)) 2 + (2< 7 i? r u;) 2 
A osc = J! (2.57) 

In particular, maximal mixing can be achieved easily by adjusting the gas pressure (and thus nn) to 
get Am 2 = 0. This is very interesting because one can benefit of using not only strong, but also long 
magnets, like those used for bending particle beams at HERA or LHC, without fear that the length could 
produce de-coherence of the beam. 

Notably PVLAS can measure also ellipticity. This depends on the difference of the indices of refraction 

of the two perpendicular components of light. Remarkably, it turns out that the sign of nil — n± is different 

for different gasses. For instance, it is positive for nobel gases like Neon or Helium while for Nitrogen 

and Oxygen is negative. This provides an independent check of the sign of the ALP coupling. Notice for 

instance the ellipticity (I2.51J) in near maximal mixing regime!! (I2.37J) for a parity- odd ALP and in the 

limit z — > 

iVsin2Co (m 2 \ 

e= I -^ - w(n\\ - n±) I z . (2.58) 

We find that if we scan the values of nn — n± by varying the pressure we will find a zero of the ellipticity 
only if riii — n± is positive, i.e. like Neon. It is straightforward to see that in the parity-even ALP case, 
the sign of the mass in (|2.58p is negative and thus the zero of the ellipticity requires a negative value of 
nil — n±, i.e. like Nitrogen. This last possibility is the one that PVLAS measurements in gas support 
empirically. 



21 The conclusions hold for any situation but the formula is much simpler in this case. 
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Figure 2.9: Degeneracy in g — m^ plane from the PVLAS rotation signal (green) 
excluded regions of the BRFT experiment [39[ [4*0] . 

2.3.6 Results and discussion 



together with the 



As was already advanced, the PVLAS collaboration has recently published [43] an excess rotation signal 
of 

|5C| = 3.9±0.5 10~ 12 rad/pass . (2.59) 

Note that, as the rotation depends on the two parameters of the ALP model, g and ms, (|2.59p cannot 
be used alone to determine univocally their values. Using values of B = 5 T, L = 1 m, and u = 1.18 eV 
we find that the degeneracy fills the space between the green curves plotted in Fig. 12.91 We see that in 
any case 

\g\ > 1.16 10" 6 x GeV" 1 (2.60) 



To provide finer values, the PVLAS collaboration combined (|2.59j) with the exclusion limits of the 
preceding BRFT experiment [391 HD] getting 



1 meV < mA, < 1.5 meV 



2 x 10 5 GeV < \M\ < 6 x 10 5 GeV 



(2.61) 



We can see this preferred region and some other smaller islands in Fig. 12.91 There we see that the BRFT 
sensibility was higher at low masses. Still the BRFT detector seems to have been placed near a minimum 
of the photon- ALP oscillation pattern. This can be considered a very unlikely casuality. Indeed, it has 
generated several complaints in the community (See for instance [55]). 

Other interesting possibility to extract these values is to combine results from rotation and ellipticity 
measurements. The PVLAS collaboration has performed several measurements of ellipticity during the 
last years and some of their results were available at the moment of the publication of [43] . However these 
measurements are far more sensitive to systematics than those of rotations and the final publication is 
suffering a considerable delay. However, measurements of ellipticity have been released in many confer- 
ences (See for instance [55j HU [561 HI] ) ■ A very interesting and useful paper [57] has recently appeared 
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providing fits of the combined results of the PVLAS rotation and ellipticity measurements together with 
the exclusion limita 22 ! from the BRFT experiments. Their best fit to the ALP hypothesis narrows the 
PVLAS values to 

m ~ 1 meV ; \M\ ~ 5 x 10 5 GeV (2.62) 

As mentioned in section 12.3.51 the ellipticity measurements point to a parity-even ALP. 

Notice the absolute values of (J2.59P and (|2.60p . I explained before that the sign of the rotation directly 
indicates the paritjo of the ALP. Unfortunately the PVLAS results [43] did not include such a sign. Only 
very recently PVLAS has made public new results that include the sign of the rotation, pointing to an 
odd-parity ALP [58]. Notice that if this claim and the ellipticity measurements are confirmed, the bare 
ALP interpretation I have presented would be ruled out, since explaining the two effects would require 
opposite ALP parities 24 !. 

The real problem with this novel extremely weak interacting particle is that it is indeed too strongly 
interacting. Although it might seem that g ~ 10 -6 GeV -1 is a rather small coupling it has very im- 
portant consequences, at least within two mature enough areas of physics: stellar evolution and tests of 
gravitational interactions. The next Chapter is devoted to show these and other constraints on ALPs. 

2.4 A further possibility: production of millicharged particles. 

There are other possibilities that have been proposed to account for the PVLAS signal. From my point 
of view, the one deserving special attention here is presented in [27]. There, the authors notice that a 
millicharged particle (MCP) of very small mass would be pair-produced from laser light propagating in 
a strong transverse magnetic field. 

This is a well-known non-perturbative phenomenon that depletes 1 1 , _L photon polarizations at different 
rates producing dichroism and hence a rotation of the polarization plane of initially polarized laser light. 
Further recombination of the pair will lead also to a birefringence effect so both PVLAS measurements 
can be justified. In this case, however, the signs of the rotation and the ellipticity do depend on the 
concrete mass and the millicharge of the MCP. Notice from Table 12.11 that if the rotation and ellipticity 
measured by PVLAS are confirmed to be of different sign the ALP interpretation would be ruled out but 
the MCP interpretation would have a chance to fit the data. 

The MCP hypothesis has been numerically analyzed in [UJ and it requires a millicharged particle 
whose electric charge and mass are very roughljo 

Q f ~ 10~ 6 ; m f < 1 eV . (2.63) 

Such couplings and mass are again strongly disfavored by astrophysical arguments, indeed the same 
energy loss arguments leading to constraints on the ALP hypothesis. I will present these and other 
constraints in Chapter [3] but I can advance that observation on HB stars in globular clusters still give 
the more demanding bound 

Q f < 1(T 14 for m f < 10 keV (2.64) 

As I will show, light millicharged particles are at the heart of our proposals [U [2]. However, the fact 
that their direct production leads to dichroism and birefringence was unnoticed by us when we wrote [2] . 



Actually they also include data from another ongoing experiment of the kind, Q&A [57], which has already taken first 
results. However, at the moment their setup is not sensitive to the PVLAS signal. 

23 In the case of an ALP without definite parity and both type of couplings (schizon), the sign indicates the relative 
importance of the two couplings g + ,g~ . See Appendix lAl 

24 This would be true even in the case of a parity non-conserving interaction. 

25 The concrete values depend on the spin of the particle. 
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The PVLAS experiment 





il> >o 


V> < o 


<5C>o 


MCPO(x < 1) or ALP+ 


MCP0(x > 1) 


<5C < 


MCPi( X » 1) 


MCPi(x< 1) or ALP" 



Table 2.1: Signs of the rotation <5£ and birefringence angle ip for the ALP and the MCP interpretations. Here 
X = | ^-^r is the relevant parameter for the MCP interpretation, u) is the energy of the photons of the laser, 

nif, e the mass and electric charge of the MCP, e is the electron charge and B the magnetic field, assumed to be 
transverse. MCPO and MCP^ stand for the case in which the MCP is a spin zero boson or a spin 1/2 fermion. 

Therefore, in spite of using this source of rotation and ellipticity, we further required an ALP coupling 
to those particlea 26 !. From this point of view the model in [2] is not minimal. 



26 This opens the possibility of the PVLAS signal to be a combination of ALP and MCP real and virtual production, 
which is nowadays being investigated. 



Chapter 3 

Constraints on novel low mass particles 
coupled to light 



This Chapter is devoted to a brief review of arguments that have been cast against the existence of low 
mass ALPs, millicharged particles and paraphotons. I divide them into three broad classes, Astrophysical, 
Cosmological and those based on laboratory direct or indirect searches. In some cases, however, some of 
them could involve two of these environments. 

Bounds on novel particles coupled to light have been studied by many authors and extensive reviews 
can be found on MCPs [59], axions or general ALPs [60|. I6H [62] and paraphotons [63]. As we will see, the 
astrophysical arguments provide the strongest bounds and typically laboratory bounds are only quoted 
for particles with masses above the stellar temperatures, harmless in stellar dynamics. However, our 
articles [H [2] have pointed out the possibility that the existence of new particles, of low mass for stellar 
standards, could help avoiding the astrophysical bounds allowing a particle interpretation of the PVLAS 
experiment. Therefore, it is important to provide laboratory alternatives to gain more information about 
these low mass particles. Recently, many articles have been concerned about these revisited laboratory 
bounds, but also new ideas have come into the field. In this section I try to review the most significant, 
a hard task nowadays since the topic is already very hot and new ideas arise quite often. 



3.1 Astrophysical bounds 

Novel particles coupling to light modify the properties of stellar plasmas, specially if their mass is light 
enough to allow kinematically their thermal production. These low mass particles, if existing, must inter- 
act very weakly with normal matter (otherwise we would have discovered them) so after their production 
they are likely to scape from the star, as neutrinos produced in the fusion nuclear reactions do. In this 
way, novel particles provide a new, invisible, contribution to the total stellar luminosity, a crucial param- 
eter of stellar evolution. Another, less comfortable, possibility is that the new particles interact strongly 
enough with the stellar medium that they are reabsorbed soon after its production. In this case they 
do not contribute directly to the total luminosity but they accelerate the mechanisms of energy transfer 
inside the star affecting notably the stellar structure. 

In this section I briefly review how we can use the established knowledge about stellar evolution to 
look for novel, low mass particles coupled to light. 
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3.1.1 Introduction: stelar evolution 

Stars form by gravitational collapse of clouds of gas, mostly Hydrogen and Helium from the early Nucle- 
osynthesis. The gravitationally bounded system losses energy by emission of electromagnetic radiation. 
Both the resulting contraction and the rise of temperature, are explained by the Virial theorem that re- 
lates the average gravitational and kinetic energy of a self-gravitating gas supported by thermal pressure 



(Ekin) — —T:\Egrav) (3-1) 



We see that a decrease of the total energy implies that the potential energy becomes more negative (the 
star contracts) and therefore the kinetic energy increases (the temperature increases). 

The configuration reaches a state of "pseudo-equilibrium" when the contraction has increased the 
temperature and density so much that allows the fusion of Hydrogen nuclei into Helium (typically T ~ 1.3 
keV and p ~ 150 gem -3 in the center core of a star with the mass of the Sun). The rates of nuclear 
reactions, depending steeply with the temperature, will produce a big injection of total energy in response 
to an increase of temperature, therefore opposing to the further contraction of the star. 

At a smaller temperature the atoms have ionized and the electromagnetic radiation generated is not 
free to escape from the star, but rather diffuses out following a random walk path. In the absence of new 
physics, the total energy loss per unit time (total luminosity) of the star is due to thermal electromagnetic 
radiation from the stellar surface and neutrino emission, mainly from the center where the star is hotter. 

A young star like our Sun burns Hydrogen until the temperature is so high at the center that also 
Helium starts to ignite. Such resulting stars burn Helium in the central core (at a typical T ~ 8.6 keV and 
p ~ 10 6 gem -3 ) and Hydrogen in an outer shell and they are called horizontal branch (HB) stars. If the 
HB star is very massive (more than 7 times the Sun, typically) it will reach an even higher temperature 
at which the Carbon and Oxygen, produced by the Helium burning, will start to fuse producing heavier 
elements. Further increasing of the temperature will cause new nuclear reactions leading to the production 
of heavier nuclei. The final configuration of this kind has a degenerate Iron corqj, (a core supported by 
the degeneracy pressure of electrons) and subsequent shells where the other light elements burn. At this 
point, the further heating of the star, due to the outer shells burning, makes the Iron core growing until 
it exceeds its Chandrasekhar limit (M ~ 1.4JAq), producing the collapse of the core. The resulting 
implosion turns into the explosion that we see as Type-II Supernovae, leaving a black hole or a neutron 
star surrounded by an expanding cloud of gas as the final stage of stellar evolution. 

However, lighter stars like our Sun will end up as degenerate Carbon-Oxygen stars, cooling by thermal 
radiation from the surface and neutrino emission from the interior (they can not release more nuclear 
power) until they finally fade out. These stars are called "white dwarfs" because being degenerate are 
very compact and thus the temperature is relatively high (a typical surface temperature is 1 keV and 
density p ~ 2 10 6 gem -3 ). 

Between every new burning phase we find an interesting intermediate state. The fusion products of 
the heaviest element store in the center of the star displacing the burning to a surrounding shell. The 
mass of the core grows until thermal pressure is not sufficient to support the gravitational pull and the 
electrons inside become degenerate. As more products are added to the core it tends to shrink, and the 
external layers of the star expand to a configuration in which the energy transport is mainly convective. 
This expansion makes the external temperature smaller and thus the color of the star appears to be 
redder, hence their name: red giants (RGjj. 



1 Selfgravitating systems have a negative specific heat. 
2 This is because the fusion of Iron nuclei requires energy. 

3 Actually the name red giant is usually restricted to stars where the inner core is made of Helium. Further developed 
stars are called "Asymptotic Red Giants" (AGB). 
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In this phase the temperature of the heavier element burning shell is determined mainly by the 
gravitational potential arising just from the mass of the inner core by means of the Virial dynamics. 
As this shell burns, the core gains mass and shrinks more and thus by the Virial theorem it rises the 
temperature of the burning shell. This feedback accelerates the shell consumption so much that these 
RG stages last much less than the first hydrogen burning phase (which is called the main sequence (MS) 
because stars spend most of their lives on it). Even when the temperature is such that the core itself 
starts to ignite, the degeneracy pressure still dominates and this energy gain does not lead to structural 
changes. Therefore, the rise in temperature is unchecked and keeps on feeding positively on the energy 
generation rate. A RG star becomes brighter than ever, reaching the so called Helium-flash. The core 
then expands nearly explosively from p ~ 10 6 gem -3 until the HB configuration is reached at p ~ 10 4 
gcm~ 3 , where thermal pressure ensures the self-regulated burning provided by the dynamics of the Virial 
theorem. 

3.1.2 The energy loss argument 

A new energy loss channel will lead to opposite effects in a star whose burning is regulated by the Virial 
negative specific heat (MS and HB) and in a star dominated by a degenerate core mass (RG). In an RG, 
the exotic energy loss opposes the positive feedback between the core and the external burning shell. As 
the energy is more easily evacuated, the temperature of the core raises slower thus allowing the core to 
be heavier at the He-flash. The main effect of an exotic energy-loss in a RG is a delay of the Helium-flash 
and thus a larger He-core mass at the following HB phase. 

In a self-regulated burning star the effects are more subtle because due to the logic of the Virial 
regulation, any additional energy loss will readjust the star into a hotter and more compact configuration 
in which the nuclear reactions can provide this new energy. In practice, because of the very steep 
dependence of the nuclear reaction rates with the temperature, only a minor adjustment is required and 
the most important implication is that the nuclear fuel is exhausted faster. Thus in the presence of a 
novel, unaccounted energy loss channel 

a star lives longer in a RG phase and less in a MS or HB phase. 

For a HB star the luminosity is almost constant for its entire life so we can argue than a small exotic 
contribution will lead to a HB phase duration of 

j std 
+ + std ^HB (o n\ 

tHB = t H B Tstd JX (3-2) 

'-'HB ~r ^HB 

Where Lhb and tHB are the luminosity and lifetime of a HB star and std holds for models without novel 
particle losses. In a white dwarf a novel energy loss will simply accelerate its thermal cooling competing 
with neutrinos. As in the case of Supernovae, this competition with neutrinos can lead to very interesting 
bounds. However, in practice these are evolved stars are more complicated to modelate. Being hotter and 
more dense than younger stars they could be thought to be more sensitive to exotic energy loss channels, 
but the models I present in this thesis suppresses the impact of such channels more efficiently at higher 
densities and temperatures. Therefore, for the sake of simplicity I will discuss none of these cases and 
focus only in the MS, RG and HB constraints, which allow for cleaner predictions, more statistics and 
consequently stronger bounds. 
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Figure 3.1: Color magnitude diagram for the globular cluster M15. V and B are the luminosities of the stars with 
filters in the visible and in blue. V is the apparent luminosity and V — B is related to the temperature. Stars on 
the left side are bluer and thus hotter. The main stages of stellar evolution are labeled: main sequence (MS), red 
giant (RG) and horizontal branch (HB). The ratio of the stars in the HB and RG stages and AV^g , the difference 
in magnitude of the brightest RG and HB stars are sensitive to new energy loss channels. See the text for details. 
Figure taken from [M] . 



3.1.3 Observations from globular clusters: HB and RG evolution. 

A globular cluster (GC) is a bound system of around 10 6 first generation stars, i.e. composed mainly 
by Hydrogen and Helium. All the stars have nearly the same small metallicit}0, because of their old 
age and only differ by their initial mass. In the Milky Way there are around 150 GCs distributed in a 
spherical halo which extends beyond the luminous disk. The color-magnitude diagram of such clusters 
reflects very nicely the star evolution outlined before. In Fig. 13.11 we can see such a diagram for the GC 
M15. Every point represents a star whose coordinates are color and magnitude, tightly related to the 
external temperature and the apparent luminosity. 

The branch called MS englobes stars, like our Sun, burning Hydrogen and increasing the temperature 
and luminosity. A star in the MS goes up and left in the color magnitude diagram. The red giant phase 
is labeled RG and it is characterized by an increasing luminosity and a progressive reddening due to 
expansion. At the tip of the RG branch, Helium starts to ignite and a self-regulated Helium burning star 
(HB) arises. They are less brighter but much smaller than late RGs so they appear bluer. The luminosity 
of HB stars is independent of their total mass, and a sub group called RR-Lyrae is used as a standard 



4 Mass fraction of elements heavier than Helium. This is a key point to produce precise statistics, because typically Z 
influences notably the star luminosity. 
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candle to determine absolute luminosities. After the Helium burning phase stars became white dwarfs 
(below the displayed region in Fig. 13. 1 j) or proceed to heavier element burning (and eventually explode 
as Supernovae) and disappear from the diagram as black holes or neutron stars. 

More massive stars evolve faster so they enter sooner the RG and HB phases. Interestingly enough, 
the RG and HB phases are much faster than the MS, so we can consider that stars in the RG and HB 
branches have approximately the same mass, i.e. RG and HB branches are an isochrone of the evolution 
of a unique star. 

To study the impact of novel particles in stars of GCs is not an easy task. One identifies observables 
in the color- magnitude diagram that can be related to the duration of the MS, RG and HB phases and 
compares these observables with numerical simulations. The difficult task is to relate these mensurable 
quantities with parameters like masses, metallicities, etc. 

The most important observables for our purposes are /S.V^ P B and the ratio R of stars in the HB and 
RG branches. Clearly R is related to the relative velocities of the RG and HB phases of a single star. As 
a new energy loss channel enlarges one and shortens the other, a negative deviation from the standard 
value of R could be the smoking gun of a new particle. On the other hand AV BB measures the maximum 
brightness of a RG (the brightness at the He-flash), which is expected to exceed a standard value if the 
degenerated core has grown larger than its standard value. 

High statistics have been collected showing the compatibility of these measurements with their stan- 
dard values from numerical simulations up to a typical 10%. Raffelt has converted this agreement into 
a limit on the novel energy loss channels in a HB, and in the RG core just before Helium-Ignition [65j . 
Curiously, both observables, when confronted with simulations, lead to the same exclusion bound for the 
averagq£| energy production rate in exotica per unit mass, 

(H) < 10 erg g-V 1 (3.3) 

for an HB or RG core. These small losses can be considered a small perturbation of the standard picture 
and therefore E can be calculated from the standard values of temperature and matter density given 
above. However, to evaluate the contribution of ALPs, paraphotons and millicharged particles to (3) 
one has to take into account the properties of the stellar plasmas. Next I address this task to end up 
deriving the corresponding bounds on the novel particles we are interested in. 



The spatial average over the HB interior or the RG core. 



38 Constraints on novel low mass particles coupled to light 



3.1.4 Dispersion relations and screening in a plasma 

Dispersion relations.- In a plasma, the presence of charged particles affects the propagation of light and 
hence the electromagnetic interactions. The calculation of photon-initiated processes in a plasma has to 
take into account these modifications. Concerning our models two aspects will be crucial, the non-trivial 
dispersion relation and the presence of screening. 

In a linear mediums the Maxwell equations in Fourier space (in Lorentz gauge d^A^ = 0) are 

(-K 2 if v + W V )A V {K) = (3.4) 

where K^ is the wave 4- vector, rf" the Minkovski metric and Ii^ v {K^) is the polarization tensor that 
accounts for the interactions of light with the medium. It has to be transverse in the sense that K^H^ = 
0, in order to respect gauge invariance. In a homogeneous medium there is a preferred frame in which 
the average constituents are at rest. In this frame I will study the propagation of electromagnetic waves 
along the z direction and use K^ = (u, 0, 0, k), sow and k are Lorentz scalars given by uj = U^K^ and 
k 2 = k 2 = (U^Kfj.) 2 — K 2 , where U^ is the four- velocity of the medium. 

The restriction to the Lorentz gauge leaves three physical degrees of freedom of the vector potential 
A^, two transverse and one longitudinal, whose polarization vectors are 

e£ = (0,l,±»\0)/>/2 ; e £ = __ . (3.5) 



kVK 2 

satisfying^ |ej| 2 = — 1. The longitudinal excitations correspond to waves with no magnetic field and 
electric field equal to El = d° Ai — d A L oc (0,0, 1), i.e. pointing in the direction of propagation, hence 
the name "longitudinal" . In the static case in vacuum they account for the static Coulomb potential of 
a distribution of electric charge, but in a plasma they generally describe "sound" waves in the electric 
field. 

Because of its transverse nature (in a 4-D sense), the polarization vector will be 

n^= J2 mfaQPT ( 3 - 6 ) 

i=±,L 

with PJ J,U = —e^e\ v the projectors into each polarization and 7Tj scalar functions encoding the response 
of the medium to the i polarization. The corresponding dispersion relations will be 

-uj 2 + k 2 + TTi(uj,k) = (3.7) 

If the medium is invariant under parity we will find that it cannot differentiate between transverse 
excitations so tt + = 7r_ = tvt- 

The polarization tensor is, from a quantum field theory viewpoint, nothing but the photon-self energy 
in the medium. In a plasma composed by electrons and other ions IP" can be computed at first order 
from the Compton scattering only from electrons [66] since the result is inversely proportional to the 
scatterer mass. This calculation is complicated because of the average over the thermal distribution of 
electrons which makes difficult to provide expressions for a general plasma. Here I am interested only in 



<; 



The response of the medium is the generation of an induced current, that in linear media is a linear function of the 
electromagnetic fields, i.e. J? nd — W^ v A v . 

7 There are some subtleties here; when K 2 — one has some freedom to define e£. Also |ez, | 2 = —1 only for space-like 
A' M , while it changes sign for K 2 > 0. 
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a classical plasma, i.e. whose electrons are not relativistic and hence T <C m e . In this case the calculation 
is very much simplified and one gets 



k 2 _ T 

.2 



u> = k +uj p 1 H „- — Transverse Plasmons, (3-8) 

oj 2 m e 



( k 2 T \ 

uj = ojp I 1 + 3—^ — ) Longitudinal Plasmons. (3-9) 

V LO z m e J 

where T is the temperature and m e the mass of the electrons. The plasma frequency wp for a classical 
plasma (in which the charged particles move slowly i.e. the temperature is much smaller than particle 
masses) is given by 

2 4-7TQn e / 5 T \ 

where n e is the electron density. We recognize immediately that transverse plasmons have a particle-like 
dispersion relation with an "effective mass" my — up while the longitudinal excitations have an almost 
fixed frequency uj ~ top. 

Screening effects.- In a plasma, the presence of an electron at a given point enhances the probability 
of finding ions near this point and, inversely, it diminishes the possibility of finding electrons nearby. 
Such polarization effect makes the medium to "screen" the electric charges of electrons and ions. We can 
compute screening effects from the Maxwell equations of a plasma in the static limit uj = 

{k 2 + 7T L (0,k))A (k)=p(k) (3.11) 

(k 2 + 7T T (0,k))A{k) =~J(k) (3.12) 

One easily finds that ttt(0, k) = 0, so stationary currents (and thus stationary magnetic fields) are not 
screened. Conversely ttl(0, k) = k 2 ^ and, notably, it does not depend on k so the potential generated 
by an localized electron (p(k) ~ 5 s (k)) is 

(k 2 + k 2 ) A (k) =p(k) -► A (r) oc -e~ kar . (3.13) 

We find a potential of the Yukawa type which "screens" the charge at distances higher than A;" 1 . A 

direct evaluation of 717,(0, /c) shows that in the classical limit the screening scale gives the well-known 

Debye-Hiickel scale 

9 , 9 4ira.n P m P 
k 2 s =k 2 D = ^^ = ^u; 2 P . (3.14) 

In contrast to their humble contribution to the dispersion relations, heavy ions contribute as much as 
electrons to the screening since k s is independent of the mass. Indeed, the screening scale is 

i 

where the sum runs over all charged species in the plasma. 

If, as an example, we calculate Coulomb scattering of electrons in a plasma, we can take into account 
screening effect by using the modified potential (|3.13p with only a minor subtlety. This potential assumes 
that the particles entering the Coulomb potential feel an average charge cloud around the target. This 
could be the situation if the probes were slow particles because the screening particles have time to move 
around in the probe's time of flight, but certainly this is not the case when the incoming particles are 
photons. To compute such a cross section one must think that the photon feels one charge configuration 
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at a time and then the average must be made after the amplitude of the process is squared. This is 
equivalent to modifying the static Coulomb propagator in the cross section 



. (k 2 + k 2 ) 4 Slow probes 

k^^l y sJ (3.16) 

k~ 2 (k 2 + k 2 )~ 2 Fast probes 
where k =\k\ is the momentum transferred. 

3.1.5 Selected processes 

Armed with the tools for calculating reactions in a plasma I shall use them to discuss the main energy 
loss channels in MS, RG and HB stars due to ALPs, paraphotons and millicharged particles. 

Primakoff effect for ALPs .- The main channel for low mass ALPs is their Primakoff production 
[67] from transverse plasmons (See Fig. 13.21 a). This calculation is addressed in detail for the case of 
the Sun in Chapter [7] to be integrated over a detailed solar model but it is valid also for HB and RG 
stars. However, the latter have a much more dense core where electron degeneracy effects turn out to be 
important, suppressing the production. On the other hand, the Sun's interior is typically much colder. 
Therefore the strongest bound for Primakoff emission of ALPs comes from HB stars. One finds from 
(|7.20p that the energy loss rate per unit mass is 

~, = ^ F{x 2 ) (3.17) 

where x = k s /2T and F(x 2 ) is a smooth function of x 2 in the cases of interest. For a HB core one can 
use F ~ 1 and (T 7 /p) ~ 0.3(8.6keV) 7 /(10 4 gcm" 3 ) leading to (E^) ~ 30 (g 10 10 GeV) 2 erg g" 1 s" 1 and 
thus to the bound [65] 

5hb < 0.6 10" 10 Ge^ 1 (3.18) 

Emission of Paraphotons: Compton scattering and Bremsstrahlung .- Massive paraphotons like the 
ones considered in Chapter Q] couple to standard model charged particles like electrons through a small 
paracharge (11.260 that depends on a possible photon mass up. In the presence of a plasma, and for small 
values of the paraphoton mass fh one has 

~2 
Tft 

Q\ * Q 1 ^— (3-19) 

Ulp 

and therefore there is a strong suppression of paraphoton interactions in a plasma if typically fh <C ujp ~ 
keV. The most important energy loss channels for a non-degenerate plasma (like MS and HB stars) are 
Compton scattering off electrons and bremsstrahlung in electron-ion collisions because they involve the 
minimum number of gauge boson vertices (and thus of powers of a). In spite of the relative suppression 
of a — 1/137, bremsstrahlung is competitive because it is not suppressed at high dense plasmas while 
Compton scattering is. 

The energy loss rate due to Compton emission can be extracted from the results in [68] for emission 
of scalar particles coupled to electrons since the total cross section (after average on polarizations) for 
vector and scalar particles only differs in a factor of two, accounting for the two photon polarizations. 
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(3.20) 
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(a) (b) (c) (d) 

Figure 3.2: Feynman diagrams for the main channels for novel particle emission from stars, (a) Primakoff emission 
of ALPs, (b) Compton production of paraphotons, (c) Brcmsstrahlung of paraphotons in electron- He collisions 
and (d) plasmon decay (7*) into millicharged particles. 

where Y e is the number of electrons per baryon (~ 0.5 in a Helium core) and m u , m e the atomic mass unit 
and the electron mass. The derived bound is in dependent. Using wp ~ 2 keV for transverse plasmons 
in typical HB conditions of density p ~ 10 4 gem -3 and chemical density I finco 



. me V \ ' 
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(3.21) 



The emission by bremsstrahlung in HB stars can be calculated again from the bremsstrahlung emission 
of novel scalar particles derived in [69l [70] and gives a slightly stronger bound 



/ meV \ ' 
e<0.1 ^^ 



\ m J 



(3.22) 



Clearly none of these bounds impose a serious restriction on the existence of paraphotonqj. Note that if 
the "paraphoton" had an intrinsic coupling to electrons (and not driven by mixing effects) the derived 
bounds would not be valid, since no suppression would be expected in this case. To recover the bounds 
in this case use m = 2 keV in 13.221 and the restriction reads typically e < 10~ 14 , as it is found in [59j . 

Millicharged particles and plasmon decay .- Contrarily to what happens with the coupling of charged 
particles with paraphotons, the small electric charge of a novel "para-particle" / would not change in a 
stellar environment and their production will be not suppressed. 

The lowest order (in a) production channel is plasmon decay into a pair of millicharged particles (See 
Fig. I3.21 d). This is possible because plasmons have a non trivial dispersion relation. Transverse plasmons 
in a "classical" plasma have an "effective mass" given by the plasma frequency (|3.10p and thus can decay 
into a pair // as long as uop > 2m f . The dispersion relation of longitudinal plasmons in a similar plasma 
fixes the energy ui ~ wp independently of the value of k, so the decay is kinematically allowed only for 
uip — k 2 > (2m f) 2 . Because of this "momentum cut-off" the longitudinal modes contribute much less 
than transverse ones to the energy loss and can be neglected here. 

The decay width of a plasmon into a pair of e-charged particles is therefore T = e 2 auip/3uj (in the 
star reference frame) . After a thermal average of initial energies one finds that the energy loss in MCPs 
per unit mass is eMCP — 5e 2 erg s _1 g _1 in a HB core. In a RG core before the He-flash the plasma 
frequency is four times larger ujp ~ 8.6 keV but the limit does not improve substantially. Accordingly, 
one must impose 



e < 2 10 



-14 



(3.23) 



Note that a more careful bound should consider the dependence of m 7 with the density. 
9 At least for those in the framework outlined in Chapter [T] 
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3.1.6 The case of the Sun 

Main sequence stars are cooler and less dense than HB and RG stars so generally the production of 
new particles is smaller and hence the bounds weaker. This is not an exception for the case of the Sun, 
despite the fact that we can measure precisely its mass, luminosity, radius, metal contents, etc. Additional 
information from the Sun's interior is obtained by Helioseismology and the study of the neutrino flux. 
This detailed knowledge has provided a deep understanding of the solar dynamics and the possibility of 
very detailed models of the present inner structure of the Sun, namely the solar models of Bahcall et al. 
. The solar bounds for a weakly interacting ALP are reviewed and restated in Chapter [7] so there is no 
need for a further development here. 

The Sun shows two further and unique characteristics that we can use to obtain information about 
novel particles coupled to light. First, the energy transfer from the Hydrogen burning core to the surface 
is thought to be radiative for almost the whole Sun. The existence of a new "not so weakly interacting" 
particle that avoids the energy loss argument being trapped in the solar interior can affect this transfer 
mechanism and hence can also be constrained. Second, we can try to detect on Earth any novel weakly 
interacting particle produced in the Sun the same we detect neutrinos. I discuss these two possibilities 
in the next sections. 

3.1.7 The energy transfer argument 

Up to now I have focused on the possibility that new particles are weakly interacting so that they leave the 
stellar interior after being produced. There might be a way of avoiding the resulting cooling mechanism 
(and the corresponding bounds) if the new particle is reabsorbed inside the star. However, even in this 
case, the particle contributes again to the stellar dynamics, leaving its imprint in this case, in the speed 
of the energy transfer from the hot core to the outside. 

There are three fundamental ways of transporting heat from the inner core to the stellar surface, 
namely radiation, conduction and convection. In the first case, thermal photons from hotter parts of the 
plasma carry the energy to the colder parts, in the second case, electrons are the carriers of heat. Finally, 
in the third one, the plasma itself moves from the place where it gets hot to the place in which it can 
release the energy. Convection is the most effective energy transfer mechanism of all three so, wherever it 
takes place, the impact of a novel energy transfer mechanism is likely to be smaller than in environments 
in which radiation or conduction dominates. 

In the Sun, the energy transfer is thought to be completely radiative except for a small region near 
the surface, so a novel energy transfer can be limited very well there. RG and HB stars are brighter and 
thus have to resort using convection somewhere. A RG star is formed by a small Helium core surrounded 
by a Hydrogen burning shell and a huge convective envelope and a HB star has a Helium burning 
core dominated by convection and further alternate layers where radiation and convection dominates. 
The complicated structure of these developed stars obscures the response to a novel energy transfer 
mechanism, although some qualitative behavior has been discussed for RGs [37| . I shall focus on the very 
much controlled case of the Sun. 

Consider a plasma in local thermodynamic equilibrium. The flux of energy due to radiation of 
particles of energy oj is proportional to their velocity f3(u>), their mean free path \(u>), and, in the 
diffusion approximation, to the gradient of their energy density £(co), i.e. 

F( W ) = -Mj^WM = -^1 X{W) '2£^ VT (3.24) 

Alternatively one commonly uses the opacity k(cj) to replace the mean free path, both being related 
by \(lu) = (k(uj)p)~ with p the matter density. If these particles are photons in local thermodynamic 
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equilibrium, which is a good approximation for the Sun's interior, we find /3(u}) ~ 1 and £(u>) = G/(] w ' — 
oo). Integrating over energies we find the total flux 

4oT 3 
F = — — VT (3.25) 

with a = 7r 2 /15. This last equation implicitly defines the effective or Rosseland mean opacity k, an 
energy average over photon energies. We can easily write down the contribution of a new particle to this 
quantity 

1 rfaftu,)^ (3-26) 



x 4aT*J mx " v ; dT 

so that K~ ta j = k" 1 + k" 1 . In the mean free path we should include all processes that can lead to 
interchange of energy between the novel particles and the Sun's thermal bath .In the calculation of the 
opacity contribution one should include also the effects of self-interactions of the novel particles that 
might slow down the energy transfen 10 !. 

Photon opacities are difficult to calculate and measure, mainly because of the atomic transitions 
of Iron and heavy elements that, despite of its small abundance, have large contributions at concrete 
frequencies. The OPAL opacity tables agree very accurately with the observed properties of the Sun [73] 
from Helioseismology so one is led to the conservative 11 ! constraint that any further exotic opacity has to 
be smaller than that of photons 

k" 1 < k' 1 ~ 1 g cm" 2 . (3.27) 

where I have used the typical opacity in the center of the Sun. 

3.1.8 Helioscopes 

Any weakly interacting novel particle produced in the Sun could be detected at Earth if it encounters 
an efficient enough detector. For instance, neutrinos from the inner Sun have been detected in several 
underground experiments using different techniques (See the dedicated review in [74]) with the common 
characteristic of requiring a large target volume to compensate the smallness of the interactions. 

In 1983, Pierre Sikivie proposed to use the coherent conversion of ALPs in a magnetic field described 
in Chapter [2] to detect a hypothetical flux of ALPs coming from the inner Sun [75] (See Fig. I3.3[) . This 
experiment would benefit from a coherent enhancement if the particle is light enough or by filling the 
magnetic field region with buffer gas to tune the photon index of refraction and restore coherence [76] , as 
was explained in Chapter [2j Three of such "Helioscopes" have been built [77J [7H1 [79] which have found 
no signal of solar ALPs. The CAST collaboration has recently published the latest and more powerful 
results [79] using a 9.3 m long decommissioned LHC test magnet reaching B = 9 T. 

To derive bounds on the ^77 coupling one has first to calculate the expected flux of <^>'s from the Sun. 
This is usually done by assuming that the emission of ALPs does not alter the solar structure given by 
the precise solar models available, but just amounts to a small perturbation. I have performed such a 
calculation in [3] (see Fig J7.4H and my results agree with previous estimates. With an essentially similar 
flux, the CAST collaboration has derived [79] the following bound 

g < 1.16 1(T 10 GeV" 1 . (3.28) 



In [7T] and [72] we find two models where self-interactions are intended to evade the astrophysical limits of the ALP 
and MCP interpretations of the PVLAS particle. 

11 It is worth noting that a detailed study of the consequences of the violation of this constraint lacks in the literature. 
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Figure 3.3: Scheme of a Helioscope. The ALP <f> (or any other particle capable of mixing with the photon inside 
a magnetic field) is produced in the solar core and leaves the Sun without further interactions. After a flight of 
roughly 8 minutes it enters a shielded long cavity permeated with a magnetic field in which it can oscillate into 
photons by inverse Primakoff effect, which finally can be detected. 



Of course, these results assume a standard ALP history: one produced in the Sun by Primakoff effect 
with an energy of some keV, the ALP goes through the Sun's outer shells, several km of "not so empty" 
space and finally the CAST shielding without any change to end up converting into an X-ray detectable 
by any of the 3 CAST detectors. The failure of any of these three hypotheses could easily evade the 
bounds_3 

Notably, these Helioscopes are also sensitive to the existence of paraphotons even in the case in which 
the regenerating magnetic field is turned off. In [63] the results of the experiment [78| are analyzed 
and some prospects are made for further experiments. The resulting bounds worsen as the mass of the 
paraphoton decreases because then the Compton production in the Sun is increasingly suppressed by the 
photon plasma mass (See section [3.1.5p or, in the oscillation formalism used in [63], the oscillation 7 — -f s 
is suppressed in presence of a large value of ujp. They get 

pV 

e < 1(T 8 ^ (3.29) 

m 

valid for 10 _4 eV < fh < 100 eV. In the range 100 eV< fh < 3 keV the bound saturates at e < 10 -10 and 
disappears abruptly for higher masses since in this cases paraphotons cannot be created in the Sun (with 
a maximum temperature ~ 1.3 keV at the core). Such a study has not been performed for the CAST 
telescope, but since the improvement on the ALP coupling with respect to this experiment is roughly a 
factor of 10, we can naiveln 13 ] expect an increase in sensibility of the bound (J3.29P of the same ordeil 14 !. 

Interestingly enough, the Helioscope bounds are much more robust than stellar energy loss bounds 
if we accept a suppression of the ALP emission from the Sun due to new physics. This is due to the 
fact that an Helioscope signal is proportional to <? 4 (e 4 ) (it involves production and detection) while the 
stellar luminosity on exotica to g 2 (e 2 ). Therefore, in a suppressing scenario the CAST bound is the more 
demanding constraint for A LPs and paraphotons. 



12 Remarkably, I have found and discussed examples where these three hypothesis could fail. If the ALP strongly interacts 
with the solar medium it can be reabsorbed or redshifted by interactions and simply keV ALPs do not scape the Sun [T1I71]. 
The ALP can have also an invisible decay channel so it might not survive the travel from the Sun to the Earth [Q. Finally 
if the ALPs have a matter-density dependent mass m^ — m t / > (p) they could be not only produced below expectations, but 
also they might be deflected by the Helioscope shielding |54l l3], 

13 This bound will depend strongly on the CAST detection techniques, of course, so a more detailed study is worth. 

14 A revision of these bounds is in progress. 
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3.2 Bounds from cosmology 

The early universe also provides a sensitive laboratory for testing the existence of novel particles coupled 
to light. In fact novel "weakly interacting" particles are possibly required by the current standard 
cosmological model which requires a rough 23% of the critical density of dark matter and a 67% of dark 
energy [73]. 

The time of Big Bang Nucleosynthesis is commonly considered the frontier between the precision 
cosmology and "terra incognita" . At temperatures of 0.1 — 1 MeV the synthesis of light nuclei takes place 
and we are provided with the oldest check of our knowledge of nuclear physics. The accurate predictions 
constrain new thermal degrees of freedom at this epoch so we must ensure that novel particles have not 
been produced at early times. 

Much later on, at a temperature around 1 eV, nuclei and electrons recombine and the universe becomes 
transparent to light. The satellites COBE and WMAP as well as a plethora of dedicated balloon and 
Earth-ground experiments have been able to measure very precisely that ancient light, called the cosmic 
microwave background. Any novel particle interacting feebly with light could distort the thermal shape 
of the spectrum and would be strongly constrained. 

3.2.1 Big bang nucleosynthesis and active degrees of freedom 

The standard theory of big bang nucleosynthesis (BBN) predicts the relative abundances of light nuclei 
(Deuterium, He 4 , Li 3 , etc.) formed in the early universe from protons and neutrons at a temperature 
T ~ MeV. These light elements are able to combine despite the opposition of an expanding universe 
which is controlled by the radiation energy density. Hence, the number of degrees of freedom of radiative 
type (whose mass is smaller than the temperature) driving the expansion can be cross-checked with 
the observed abundances to limit possible exotic degrees of freedom. The deviation from the standard 
number of effective degrees of freedom gf" = 10.75 at that epoch is constrained to bqlj [H] 



0.-rf* = £* (|) 4 + lY,9i (I)' <2.8 (3.30) 

B ^ ' F ^ ' 

where I have included a sum over possible new bosons (B) and fermions (F) and gi is the number of 
internal degrees of freedom of the particle i. Here T{ is the temperature of the species i, which can be 
lower than T if it has decoupled before the reheating caused by annihilation of heavier standard particles 
(like QCD resonances, /i ± ,r ± ,etc). 

We see that a spinless boson (g = 1) or a Weyl fermion {g = 2) are compatible with the bound even 
with Tj = T, but this is not the case for a Dirac fermion (g = 4) or a massive paraphoton (g = 3). In this 
thesis I try to show that in order to solve the PVLAS- Astrophysics inconsistency we need new physics at 
a low energy scale, typically new particles and interactions. In fact, my best model includes at least two 
paraphotons and a Dirac millicharged fermion (although it could be a boson as well) with typical masses 
smaller than 1 eV. 

There are two arguments that can be used to circumvent these bounds. On the one hand we can 
assume that such particles are not produced in the primordial soup (coming form the decays of an inflaton, 
for instance) so they don't exist at the BBN epoch. Of course, we have to worry about its possible thermal 
production after inflation by the interactions with normal matter that I will introduce. The requirement 
of the freezing of these interactions all through the early universe (up to the BBN epoch) provides bounds 
for these new interactions that my models have to satisfy. On the other hand we can assume that these 



15 Stronger constraints can be derived by using populations out of equilibrium as in [80] 
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particles were produced indeed in the early universe but their interactions with standard model particles 
freezed out before BBN. Therefore, their temperature can be lower than that of photons and they can 
satisfy the bounc0 (|330l) . 

A particle species decouples when the rate of the interactions that keeps it in equilibrium with 
photons, r, is smaller than the rate of the expansion of the universe, H. In the early universe the rate of 
expansion is proportional to the squared temperature, H oc T 2 . However, the typical interaction rates of 
ALPs on the one hand, and MCPs and paraphotons on the other, behave very differently with T. ALP 
interactions proceed through the irrelevant couplings in (|2.4p and thus will have a strong dependence with 
temperature. From dimensional grounds and neglecting the ALP mass Talp °c 9 2 T s , while MCPs and 
paraphotons interact through marginal operators what will typically provide Tmcp~j °c T. Therefore, the 
ratio T/H behaves very differently for both cases: it increases with temperature for ALPs and decreases 
for MCPs and paraphotons. Consequently, ALPs in thermal equilibrium in the early universe could 
decouple before BBN but this is not the case for MCPs and paraphotons. 

Therefore in what follows I will focus on the first possibility for MCPs and paraphotons and in the 
second for ALPs to derive the required constraints from fj3.30f) . Notably most of my models in [JJ [2] 
are constructed to suppress the couplings of these particles to normal matter in the hot and dense 
stellar plasmas. In the more extreme primordial plasma of the early universe at the epoch of BBN these 
interactions will be even more suppressed, relaxing very much these bounds. 

ALP.- Although the current bound (13.30P allows for the existence of a spin zero boson it is interesting 
to know what would be the limit on the coupling g if the limit strengthens to g* — gl td < 1. Indeed this 
was the case some time ago when the general study [61] was performed. There we can find that ALPs 
with a coupling 

g > 2 1(T 7 GeV" 1 (3.31) 

will decouple before the reheating of photons due to pion and muon annihilation (T < 200 MeV) and 
therefore will have a smaller temperature, evading the BBN bound. For this bound, the value g* — gf 1 < 
0.5 was used [53]. See [61] for details. 

Lonely Paraphoton Model- Consider an extension of the standard model based just on a new gauge 
symmetry U(l)±, with its corresponding paraphoton and some sector responsible for its mass on the basis 
of Chapter [TJ 

A background of paraphotons will form if they couple to electrically charged particles with a small 
paracharge e, for instance, by Compton production je — ► j'e. To avoid such a production, the rate of 
this interaction T ~ e 2 a 2 T has to be smaller than the expansion rate of the Universe H ~ T 2 /mpianck , 
but note that in such a plasma photons have an effective mass m 2 = Wp ~ 47raT 2 /9 and thus, by means 
of (jl.26p the coupling is suppressed (e = eT=om 2 /u)p) with respect to its vacuum value. Finally one gets 



, / T T ^ s MeV 

< a-\ < HT 8 3.32 

V rapi anck m m 

where for the last equality I have used a minimum temperature for BBN of 0.1 MeV and mpi anc k = 10 19 
GeV. This plasma-mass suppression makes (|3,32p completely harmless for small m as it happens with 
the energy loss in stars [63], see equations (|3.2ip and (|3.22p . 

Paraphoton + Millicharged fermion Model.- Consider the extension of the later situation with a new 
Dirac fermion charged under a new symmetry £7(1) i. This particle will acquire an e-sized electric charge 
given by (|1.25|) . 



In this case they can contribute to the dark matter of the universe and eventually its energy density could exceed the 
accepted value Q ~ 1 and overdose the universe. This argument has been used to obtain bounds on ALPs [61] and MCPs 
82 , 83] but they apply typically for ALP and MCP masses above 10 eV and therefore are not interesting for this work. 
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Not surprisingly, a background of millicharged fermions forms much easier than one of paraphotons 
because the small electric charge of millicharged particles (|1.25|) is not suppressed by a relatively large 
photon plasma mass. The rate of plasmon decay into a fermion-parafermion pair is ~ ^r = QG?T and then 
we have to demand 

e T =o<a-\[—^<10- 8 (3.33) 

V "^Planck 

to prevent the formation of a background of millicharged particles. In order to achieve some suppression 
for this and other channels in stellar plasmas I invoked a model with two paraphotons which is presented 
in [2]. There we again recover a wp-dependent electric charge of parafermions in a plasma, and the bound 
is essentially similar to (|3.32p (regardless of small factors of a, 7r,etc.) and again completely harmless. 

3.2.2 Cosmic microwave background 

Novel particles weakly coupled to light during and after the time of decoupling of matter and radiation 
can distort the CMB Planckian spectrum which has been precisely measured. At temperatures around 1 
eV, the electrons and protons recombine and the universe becomes transparent to light. It seems therefore 
that the dispersion relation of photons is at most due to the neutral atoms and the few still unpaired 
charged particles, which give a negligible contribution. The CMB photons satisfy therefore the vacuum 
dispersion relation q 2 = 0. 

In this situation ALPs would be produced by photon coalescence (77 — ► <f>) or in pairs (77 — ► <f>4>)\ 
MCPs would be produced in pairs (77 — ► //) and paraphotons could be produced easily from them by 
Compton scattering (7/ — ► 7/). 

Very recently these arguments have been used to rule out millicharged particles with mass m < eV 
and millicharge e < 10~ 7 |85| . although the bounds can be much stronger (e < 10~ 9 ) in models with a 
paraphoton, like the ones considered here. The same article mentions that the value of the ALP coupling 
candidate for explaining PVLAS, namely eq. (|2.6ip . is still allowed by current precision of the FIRAS 
experiment at the COBE satellite. 

There could be ways to circumvent these bounds. For instance, it seems that if the millicharged 
particles (or paraphotons) had thermalized after BBN, but before the CMB epoch, there would be no 
distortion of the Planckian spectrum at decoupling and therefore no bounds 17 !. 

3.3 Laboratory bounds 

The powerful astrophysical arguments, specially the energy loss one, restrict very much the parameter 
space of novel low mass particles coupled to light, but they are certainly model dependent. The searches 
for these particles under the controlled conditions of an Earth-based experiment have the disadvantage of 
much smaller statistics but they benefit from a better precision, a controlled environment and therefore 
a higher reliability of the results. In this section I very briefly review the most relevant laboratory 
experiments constraining the existence of ALPs coupled to light, paraphotons and millicharged particles 
present in the literature. I divide the experimental efforts into collider and low energy experiments. 

Experiments at accelerators typically benefit of a high amount of energy available in the center of 
mass so that very massive particles can be produced. At the same time, clean collisions are rare and the 
lack of statistics makes impossible to look for very weakly interacting particles like the ones considered 
here. A higher precision is typically achieved in rare decays or spectroscopy of bound systems. We 
will see that Positronium decays and measurements of the Lamb shift give very strong bounds for these 



7 The impact of such possibility seems not to have been addressed by the moment. 
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Figure 3.4: ALPs with an ^77 coupling are produced in e + e colliders. 

weak interactions. Optical experiments with lasers offer high statistics and the possibility of coherence 
enhancement if the novel particles mix with photons, but they are performed at low energies where 
backgrounds can be difficult to control. Finally, in the cases in which the novel particle can mediate long 
range forces between macroscopic bodies, experiments testing the inverse square law of the Coulomb and 
Newton forces offer the most impressive bounds due to the cumulative effect over a macroscopic amount 
of particles. Let me review briefly all of them. 

3.3.1 Bounds from colliders 

If novel weakly interacting particles take part in the products of a high energy collision, the detectors will 
identify missing energy and momentum when reconstructing the standard particle tracks. For instance 
in the SLAC storage ring PEP a dedicated detector (the ASP, for Anomalous Single Photon detector) 
was built with this purpose. From the absence of e + e~ — ► 7 + p events one can bound a ^77 interaction 
like (fLT6|) (See Fig. E3D getting [EE] 

g < 5.5 IO^GcV- 1 (3.34) 

A recent paper [86] has analyzed these results and derived new bounds from the different LEP and 
KEKB experiments. Constraints as strong as g < 1.6 10~ 6 GeV~ , which will be sensitive to the PVLAS 
ALP, are suggested for LEP and KEKB experiments. However the missing energy analysis has not been 
really performed with LEP and KEK data, so these bounds are derived under the assumption than no 
deviations from the standard model predictions are to be found. 

Beam dump experiments have been used to limit the existence of millicharged particles. In such 
experiments a beam of electrons is shone against a target and the products remaining after some meters 
of shielding are analyzed (See Fig. 13. 5p . In [83] the results of the SLAC Beam Dump experiment were 
used to constrain millicharged particles with masses smaller than ~ 0.2 GeV 

e < 3 10~ 4 . (3.35) 



Also, very recently it has been reported a search for millicharged particles in reactor neutrino detectors 
|87| . yielding a slightly improved bound 

e < 10~ 5 . (3.36) 

3.3.2 Accelerator cavities 

Particles coupled to light can be produced in accelerator cavities like the ones used in big colliders 
to accelerate charged particles. In a recent publication [88] the Schwinger production of millicharged 
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Figure 3.5: Schematics of the SLAC beam dump experiment. 

particles in such cavities has been investigated concluding that low mass millicharged particles with 
e > 10 -6 would leave a noticeable imprint on the quality factor of the studied cavities of the TESLA 
experiment [89]. The prospects for future cavities would lead to ameliorate this bound in an order of 
magnitude. 



3.3.3 Optical experiments 

The PVLAS experiment, subject of Chapter [21 is only an example of the several experiments that can 
be performed searching for low mass particles interacting with light with lasers in optical setups. The 
simplest version of the family is the so called light shining through walls (LSW) experiment [9CH I91j . 
The basic idea is depicted in Fig. 13.61 There we see how a photon from a laser oscillates into an ALP 
in the presence of an external magnetic field B, then this ALP crosses an opaque wall to be further 
reconverted into a photon which can be easily detected at the end of the regenerating magnet. In 
contrast to the PVLAS type of experiments the signal would be due to the appearance of photons, so a 
positive result would provide a strong evidence for new particles. Because the right side of the opaque 
wall can be strongly isolated, the backgrounds of such experiments can be very much reduced. This is a 
key feature because in the PVLAS-type of experiments the disappearance of photons is proportional to 
the 077 coupling squared, while in LSW experiments the signal goes with the fourth power of such small 
quantity. 

The rate of regenerated/incoming photons can be easily calculated from the formulae derived for the 
probability of 7 — (p oscillation (12. 24ft in a magnetic field derived in Chapter [2j 
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(3.37) 
(3.38) 



With £ the angle between the laser polarization and the magnetic field direction (both of them assumed 
to be transverse to the beam direction). The 7 — eft mixing angle, 8, and A osc are defined in section [2T3 
and Li ; 2 are the lengths of the first and second magnet, respectively. The outgoing photons will have 
polarization parallel to the magnetic field if the ^77 coupling is odd, and perpendicular if the coupling 



is eve 



18 If parity is not conserved by the ^77 interaction, £ is the angle between the laser polarization and the schizon angle 6 B 
(See Appendix [S} . The polarization of the outgoing photons will be along the 9 B direction. 
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Figure 3.6: Schematic view of a "light shining through a wall" experiment (LSW). ALP production 
through photon conversion in a magnetic field (left), subsequent travel through an opaque wall, and 
final detection through photon regeneration (right). The same scheme holds for a sterile paraphoton 7 S 
replacing <j) without the external B field insertions. 



The BRFT collaboration performed such an experiment ending with no signal of novel particles [3DJ . 
This leads to the bound plotted in Fig. 12.91 which for mC 1 meV is 

g < 3.6 1(T 7 GeV" 1 (3.39) 

Notably, LSW experiments are also sensitive to paraphoton oscillations even without magnetic field. 
From eq. (|1.36p we find that the probability of photon regeneration is 

Nc ° Unt ^ ~ 75 ~ 7) = P( 7 - 7-)P(7- - 7) = (3-40) 

-< 'incoming 

a ( m 2 \ 9 Am 2 Li ,Am 2 L 2 

= 16e 4 9 1 9 sin 2 — sin 2 — 3.41 

Vmg-mf/ 4a; Auj 

Where now L\ 2 are the lengths from the photon source to the wall and from the wall to the detector and 
Am 2 is the square mass difference of the two propagating states. One might wonder if 7 — j s oscillations 
can be enhanced in the presence of a magnetic field due to low mass millicharged particles in loops [92] . 
The BRFT collaboration performed such a LSW experiment with null results [40] . Using wiq <C vn\ 
their result can be quoted as 

e < 1.7 W 5 - 4 llr '' V V (3.42) 

V mi y 

for m\ < 5.4 10 -5 eV. From this mass up to mi < 1.9 10 -4 eV, the probability of regeneration starts 
to oscillate and the bound (not quoted in [3D]) saturates at e ~ 1.7 10 -4 . Beyond that mass the BRFT 
apparatus was not sensitive to paraphoton oscillations. 

Other optical experiments could be sensitive to the interactions of ALPs, MCPs and paraphotons with 
light. Recently photon splitting has been proposed to test the PVLAS signal [93] and some experiments 
for 37 — ► 7 coalescence have been proposed to measure the QED vacuum effects (See [M] and references 
therein). Also, a new experimental set up has been proposed to measure vacuum birefringence [95] and 
some technical improvements have been suggested to increase the reach of the current LSW experiments 
[96] 197]. It seems that optical experiments still have much to say about ALPs, paraphotons and MCPs. 



3.3.4 Ortopositronium decays 

Ortopositronium (oPs) is the 1 3 S bound state of an electron and a positron. Being a spin 1 state, it 
cannot decay into a pair of photons because of the Yang-Landau theorem [98] [99] and its main decay 
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Figure 3.7: Contribution of ALPs (a) and millichargcd particles (c) to the invisible and semiinvisible decays of 
ortoPositronium. 

mode is oPs^ 37. This is a slow process that favors the discrimination of additional decay channels into 
ALPd_| [100] and millicharged particles [10141102] (see the diagrams depicted in Fig. 13. 7|) . Experimentally 
one finds 



BR(oPs -► 7 + invisible) < 4.2 10" 7 [103] 
BR(oPs -> invisible) < 8.6 10~ 6 [102] 

and using r(oPs) = 7.05^s~ 1 jlOlj and the widths 
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(3.43) 
(3.44) 



(3.45) 
(3.46) 



(where <fi are the ALPs that couple to photons according to (|1.16[) and (|1.18j) ) we find the following 
constraints on the relevant couplings 



< 1.5 KT d GeV" 



< 1.9 10 -J GeV" 



e < 3.4 10" 



(MCP) 



(3.47) 



There are proposals for improving the limits [1UU|, I102J by an order of magnitude (See the references 
in [27]). Finally note that the new vertices are probed at a momentum transfer squared q 2 = (2m e ) 2 , 
much higher than the energies related to the PVLAS experiment. As our models [TJ [2] are designed to 
decrease the novel interactions at these high energies, these bounds would turn out to be negligible in 
such models. 



3.3.5 Radiative impact on QED precision measurements 

Apart from their real production (and possible detection) in accelerators, optical setups or oPs decays, 
novel particles can affect QED precision observables radiatively, i.e. through its virtual production. The 
most powerful measurements of QED are the electron anomalous magnetic momentl 20 ! and the Lamb shift, 
and both of them have been used to constrain novel particles coupled to light. 



19 This reference does not consider the coupling ^77 but a direct Yukawa coupling of the ALP to electrons. 
The muon g — 2 is much more sensitive to contributions of particles with masses beyond the muon mass because in that 
case gi — 2 ~ rn^/M 2 with M the heavy mass. For smaller masses one can benefit from the fact that g e — 2 is measured 
more precisely. 
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Figure 3.8: Contributions to the electron g e — 2 from ALPs through the ^77 coupling (left) and from vacuum 
polarization produced by a millicharged particle /. (There are two additional diagrams to the ALP contribution 
differing in the way the photon lines are linked with the electron) 

Anomalous magnetic moment.- Aa ALP with a 0^77 coupling will induce also a contribution analogous 

to the pion-pole light by light one. This contribution can be extracted from a number of papers studying 

the pion contribution (See for instance |104j ) but it turns out to be a hard task for a poor bound. 

Therefore, and only for the sake of completeness, I provide just an order of magnitude bound estimating 

that 

' g e - 2\ a ml 



7T 



2 M 2 



(3.48) 



Similar expressions are obtained for the muon g^ — 2. 

A millicharged fermion / will also induce a contribution (See Fig. 13.81 for the relevant diagram). 
Keeping only the leading logarithm from the general expression, that can be obtained from [105] I find 



* V =T Log (3 - 49) 

\ 2 / 3 \-kJ rrif 

Experimentally, one finds the impressive result [106] 

^-— = 0.0011596521859(38) (3.50) 

which allows to set the following constraints accepting that the experimental value corresponds to the 
SM expectation^ 

g > 10~ 3 GeV" 1 (ALP) ; e < 3 10~ 4 (MCP) (3.51) 

where I have used mj = 1 meV for the last bound. Constraints from the anomalous magnetic moment 
of the muon are weaker. The contribution from paraphotons turns out to be negligible. 

Lamb Shift- Novel particles coupled to light can contribute also to the energy levels of atoms. Al- 
though it is not the most precisely measured shift! 22 !, the Lamb shift AE^s = E>2s 1/2 — E2P 1/2 of Hydrogen 
has been used to look for massive |80j {m,f > m e ) and light millicharged particles, because disentangling 
their contribution is easier. The calculation of the contribution from a (^77 coupling has not been ad- 
dressed by the moment. After a mild calculation the contribution for light millicharged fermions gives 

AE f = -^e* . (3.52) 

; 18vr v ; 



1 This argument is a bit too naive because in practice (g — 2) e it is used to determine the fine-structure constant a which 
is required to compute the SM prediction. 
22 The 1S-2S shift is better determined [T07l. 
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Figure 3.9: Contribution to the Lamb Shift from vacuum polarization from a novel millicharged fermion. 



Assuming that this additional splitting lies between the experimental uncertainty of the most precise 
measurement of the 2S 1/2 -2P 1/2 Lamb shift in Hydrogen pJH EDI EI] AE cxp /(2Trh) = (1075.85 ± 
0.01)MHz we find 

' 2.7 10" 4 (MCP - 0) [108] 
e < \ (3.53) 

1.1 10~ 4 (MCP- |) [TD5] 



3.3.6 Limits on a 5 force 

The exchange of light bosons between matter fields leads to long range forces between macroscopic bodies. 
The electrostatic Coulomb force and Newton's law of gravity are the only known forces of such type and 
they are well understood in terms of photon and graviton exchange. Both forces are extremely well tested 
experimentally and these tests can be used to constrain novel bosons coupled to electrons and protons. 
Goldstone bosons, however, only mediate spin-dependent forces at first order [112] which are ineffective 
for tests using unpolarized bodies. The reason is that Goldstone bosons couple derivatively to fermions 
and then a parity-even structure gives d^cfii/j^/^ip = by means of the Dirac equation of motion of the 
fermion 23 !. On the other hand, the parity-odd structure d^cpip^^j^ip = 2m^V , 75'0) which in the non- 
relativistic limit ~ 2m,/, a ■ q , where a are the Pauli matrices representing the spin and q is the 
momentum transfer carried by <fi. Although some work has been carried out to test these spin-dependent 
forces, the most accurate measurements are still those testing unpolarized bodiea 24 !. 

Regardless of its nature, I am concerned with ALPs that couple to light with both parity structures 
(jl.lfip and (|1.18|) . Chapter [5] is devoted to examine the existing bounds on these particles. 

Massive paraphotons will also couple to charged particles (recall Chapter [1]) as electrons and protons 
providing an electromagnetic type 5 force. The "paracharge" of a standard fermion will be in this case 
proportional only to its electric charge, and hence neutral bodies will also have zero paracharge. Then, 
paraphotons are to be looked for by providing a net electric charge to the test bodies. The reported 
negative results again place bounds on the e — m plane that can be seen in Fig. 13.101 



23 It is evident that this argument does not hold for flavor changing currents, but then we will not expect a "force" in the 
classical sense because the probes will change its flavor. On the other hand we can not forget that non-derivative interactions 
can be generated by terms that break explicitly the symmetry associated to the Goldstone. 

24 It is interesting to note that exchange of two goldstones [113] or two fermions [751 1114] leads to a spin-independent force 
again, and thus tests of spin independent forces lead to (weaker) bounds on spin dependent ones. 
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Figure 3.10: Constraints on an electromagnetic 5 th force. j3 corresponds to our e 2 while A 
from [115] . Labels in the figure can be found in the references cited in [115] . 



Figure taken 



3.4 Summary plots 

In order to summarize this Chapter I display the bounds on ALPs, paraphotons and MCPs in Figures 
I3.11l and l3.12i In Fig. I3.11l we see clearly that the hypothetical PVLAS ALP is strongly excluded by the 
HB lifetime argument, Helioscope experiments and 5 force searches (this last one holds only for a + 
ALP). As we see in Fig. 13.121 the MCP possibility (rriMCP < 1 eV and e > 10~ 6 ) is again very disfavored 
by the HB lifetime argument but also by distortions in the CMB spectrum and the number of radiation 
species at the epoch of BBN. With the possible exception of the 5 force searches, laboratory searches are 
consistent with both hypothesis, being the BRFT experiment the only one having competitive sensitivity 
with the PVLAS apparatus. 
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Figure 3.11: Summary plot for the constraints on ALPs. The 5 th bounds only apply for a parity-even ALP (0 -1 
The PVLAS signal is shown in green. 
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Figure 3.12: Excluded regions for the mass and coupling of paraphotons (up) and for spin-1/2 millicharged 
particles (down). For spin-0 MCPs the BBN bound disappears and the rest remain very close to the spin-1/2 case. 
Sec the text for details. 



Chapter 4 

Can the PVLAS particle be compatible 
with astrophysical bounds? 



In Chapter [3] I have shown that bare ALP and MCP interpretations of the PVLAS signal, involving 
each of them just a new particle and its corresponding interaction to light are completely excluded by 
astrophysical arguments (a parity-even ALP is also excluded by 5 th force experiments, as it will become 
clear in Chapter [8]) . In this thesis I have focused on the ALP interpretation and the general question I 
want to answer is: 

are there modifications of the ALP bare physics that avoid the astrophysical bounds? 

The answer is a modest "yes" which I am intended to explain. However, I can advance that my 
general strategy for invalidating the astrophysical bounds is to accept that these interactions are not the 
end of the history. Let me state the number of hypothesis in the PVLAS-astrophysics puzzle: 

• 0- </>'s are elementary. 

• 1- (f>'s are produced by Primakoff effect in stars by the couplings (|2.4[) 

• 2- m and g have the same values in PVLAS than in stellar conditions. 

• 3- We use unperturbed standard stellar models to calculate ALP emission. 

• 4~ 4 > ' s escape from the solar interior. 

• 5- We trust astrophysical calculations of evolutionary stellar timescales. 

• 6- </>'s from the Sun enter into Helioscopes and convert into photons by inverse Primakoff conversion. 

Note that there are two kinds of hypothesis, the ones that affect our knowledge about the ALP physics 
(0,2), completely unknown to us, (except for the existing bounds and maybe the PVLAS signal) and the 
ones regarding aspects of established fields of physics (1,3,4,5,6). As a natural strategy we will first have 
to try to adapt the new ALP to already understood physics and not the inverse, namely to adapt well 
established phenonrena or theories to englobe the bare ALP existence. This is well justified if one realizes 
that the ALP ^77 coupling required by PVLAS (I2.61J) violates the astrophysical bounds by five orders 
of magnitude!. It is not very likely that a soft modification of our orthodox ideas would leave room for 
the PVLAS ALP. 

In this thesis I have worked and discussed two general ways of evading the astrophysical constraints 
providing some new physics involving the PVLAS ALP, namely 

• The ALPs produced in the Sun are trapped inside, 
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• The Primakoff production in stellar conditions is suppressed. 

Trapping ALPs .- The idea consists in providing new (strong) interactions of the ALP with some of the 
constituents of a stellar plasmalj. Instead of escaping, and contributing largely to the stellar energy loss, 
the ALPs thermalize and slowly diffuse out to contribute typically as much as photons, not compromising 
very seriously the evolutionary time scales of stars. The diffusion of novel particles is, however, also 
constrained by the energy transfer argument and this imposes a new constraint on the feebleness of the 
novel strong interactions, as I review in Chapter [3j This possibility demands reconsidering hypotheses 3 
and 6 so it is likely to generate some problems. 

Suppression of stellar ALP production .- This idea supposes that the ^77 interaction is suppressed 
in stellar environments with respect to the PVLAS experimental setup. This can be achieved in several 
ways depending on the parameter that makes the difference. If we accept that the ALP interaction is 
always of the type (|2.4p we find that the coupling g in the Sun has to be roughly five orders of magnitude 
smaller than in PVLAS to evade the HB bound (|3.18p and about 10 orders of magnitude smaller to evade 
the CAST detection limit ([525]) . 

I this thesis I have considered both possibilities and indeed I have built models for both of them. 
However providing a strong coupling of <f> to the stellar plasma constituents is not an easy task if we 
examine its consequences out of the stellar environments. The second idea has showed up to be more 
successful, and to my knowledge the model [2] provides the only known way to reconcile PVLAS with 
particle physics and astrophysics. Next I dedicate a few words to contextualize the work exposed in 
Chapters 5 to 8, corresponding to the four articles in which I have published the main results of my PhD 
research. 

The key point in all my work arises when we realize that the discrepancy between PVLAS and 
astrophysics is really huge. In order to reconcile them I believe we need to require that the new physics 
involving the PVLAS ALP behaves very differently in the PVLAS and the stellar environments. 

If we consider eft and its coupling to two photons from an effective field theory point of view we can 
interpret the scale g~ l = M ~ 10 6 GeV of the (^77 coupling as the typical scale of the new physics. 
This is a very large scale, and the difference between the PVLAS and stellar conditions is very small 
when compared to it. If this were the case the ^77 interactions would be described fairly by (|2.4p up to 
energies of 10 6 GeV and any chance of reconciling PVLAS and astrophysics would be lost. 

From this viewpoint a reconciliation of the PVLAS ALP interpretatioro with the astrophysical con- 
straints requires that: 

the energy scale, Atvp, at which new physics appears giving the interactions measured in PVLAS, 
(|1.16p or p,18|) . should lie between the energies involved in the PVLAS experiment and those involved in 
the ALP emission from the stars, i.e. Ap vlas < Anp < A s t ar s- With A 8 t ars of the order of keV it is 
clear that the new physics should be at a very low energy scale. 

But, how to reconcile a small Anp with the value from the PVLAS measurements giving M ~ 10 6 
GeV ? Up to now we have only found one solution that all models in this thesis share. It requires the 
addition of a paraphoton 7' that mixes kinetically with the conventional photon 7 with a strength e. 
Then I made the hypothesis that the effective interactions (^77 come from a more fundamental ^yy'j' plus 
kinetic mixing of the two gauge bosons. If the 4>j'j' interactions come from the equivalents to (|2.4j) with 
paraphoton field strength F^ u replacing the photon field strength F^ u and g' is the coupling of the <fyy ,r )' 
vertex then we find that 

g' = e ~ 2 g ~ e" 2 10" 6 GeV" 1 (4.1) 



A further idea is that rf> interacts with dark matter accumulated in the Sun and other stars and not with the plasma. I 
regarded some models without success. Another idea was presented in [71] . 
2 Given in the terms of Chapter [2] 
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Figure 4.1: The structure of the 4>jj vertex responsible for the PVLAS signal in our models [TJ[2]. 

Using Anp ~ g'^ 1 ~ 1 eV we need to require a typical e ~ 10~ 7 . 

The nature of the (jr/^f interaction is what makes the difference in all the models proposed in [TJ [2J [3] 
(Chapters 5 to 7). While it is not discussed in the first model of [1], in the second we already introduce 
the possibility that it is produced by a triangle loop in which a novel low mass fermion, /, runqj (See Fig. 
14. ip . We couple 4> directly (and only) to this fermion, which is paracharged. This triangle loop is also at 
the heart of the model presented in [2]. Finally, in [3] we perform a general analysis without specifying 
the origin of the ^77 interaction. 

All this reasoning holds if the important parameter is the energy scale at which the vertex ^77 
is probed. But it might be as well that other environmental parameters would be the responsible for 
the suppression. For instance in [116] a model is proposed in which at the high temperatures of the 
stellar interiors the new physics overcome a phase transition that turns the ^77 coupling to zero. The 
temperature plays here the role of the energy in the above discussion and we should require Tpvlas < 
Tq < T s tars- Further examples are provided for instance in Chameleon models J117J . There, the mass of 
a scalar field, the chameleon, turns out to be dependent of the energy density in this point, which in a 
non-relativistic universe is nothing but the matter density. Let our ^77 coupling depend on such a field 
and we will get a matter density dependent coupling. 

Moreover, the same model proposed in [2] (Chapter [6]), albeit with some subtleties, produces a coupling 
which depends on the plasma frequency ujp. 

Needless to say, all the values of these parameters, temperature, matter density, photon plasma mass, 
and others, are very different in the PVLAS and stellar conditions (See Table I7.ip so hopefully models 
can be built in which such a difference is the responsible for the suppression of the production of ALPs 
in stars. These very interesting possibilities could be studied altogether in a model independent way and 
I have performed such an analysis in [3] focusing on the case of the Sun. The reason has already been 
mentioned, given a nonstandard suppression of the Primakoff ALP flux, the CAST bound scales with the 
fourth root while energy loss arguments only with the square of the suppression. This work is presented 
in Chapter [71 

Finally, in Chapter [8] we noticed that if the PVLAS ALP is parity-even the ^77 vertex will induce 
radiatively long range forces between macroscopic unpolarized bodies. In this case, our work gives the 
most demanding bound for the ALP coupling, around 5 orders of magnitude more stringent than the 
best astrophysical bounds for the mass suggested by the PVLAS experiment m ~ 1 meV. In the same 
paper we show that our model [2] avoids this constraint. 



We use a fermion for simplicity and because low mass fermions are more natural from the effective field theory viewpoint. 



Chapter 5 

Evading astrophysical constraints on 
axion-like particles 



(In collaboration with E. Masso. Published in JCAP 0509, 015 (2005) [T]) 



5.1 Introduction 



An pseudoscalar axion-like particle (j) coupled to photons 

^77 = Jm^^ ( 5 -!) 

(F^ u is the electromagnetic field tensor, F^ u its dual, and 4> the axion-like field) would be able to transform 
into photons when electromagnetic fields are present 

C^ = ^E-B4> (5.2) 

This is completely analogous to the well known Primakoff effect that involves the 7r°77 coupling. In this 
paper we will be interested in the case that (f> is very light since then a number of interesting effects may 
happen. 

When the <fi mass itij, < T c , with T c ~ 1 — 10 keV the typical temperature of stellar cores as those 
of the Sun or horizontal-branch stars, <j) particles are produced by the Primakoff-like effect due to the 
interaction (J5.ll5.2p . If one further assumes that the produced (ft flux escapes freely from the star and 
thus constitutes a non-standard channel of energy-loss, the strength of the interaction (|5.1|5.2p can be 
bounded using observational data on stellar evolution time scales [118] . For the Sun, one has the limit 

M > 4 x 10 8 GeV (5.3) 

There is general agreement that these arguments applied to globular clusters lead to an even stronger 
bound 

M > 10 10 GeV (5.4) 

(Here and thereafter, we understand that the bounds are on the absolute value of M.) 

In the allowed range for M, axion-like particles are still produced in the Sun and a calculable flux 
reaches the Earth |76^ I119J . A proposal to detect this flux was given in a pioneer paper by Sikivie [75j . 
The idea is that the interaction (|5.1|5.2j) allows solar (/>'s to transform back into X-ray photons in a cavity 
with an external magnetic field. Such helioscopes have been built and limits on M have been obtained 
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from the non-observation of this inverse Primakoff process in the cavity [78] [77] 11201179] . The last result, 
from the CERN Axion Solar Telescope (CAST) collaboration [7J5], is 

M > 0.9 x 10 10 GeV (5.5) 

(95%CL) comparable to the bound based on the stellar energy loss arguments, eq. (15.4p . The strong limit 
(|5.5p can be established when a coherent (f> — ► 7 signal is expected, which happens for tua, < 0.02 eV. The 
CAST prospects [79] are to further improve (15.51) and to extend the results to masses m^ up to 1 eV. 

The research we present in this letter has been motivated by the observation of a rotation of the 
polarization plane of light propagating through a transverse, static, magnetic field by the PVLAS col- 
laboration [43] . A possible interpretation of this result is the existence of a light axion-like particle <j) 
coupled to two photons [33]. However, if interpreted this way the scale appearing in (|5.1|5.2p must be 

M ~ 4 x 10 5 GeV (5.6) 

It results in such a strong coupling that it is in contradiction with the bounds (|5.4p and (|5.5p . Yet, it is 
consistent with the bounds coming from particle physics experiments |62j. 

Let us stress that if there exists a particle with the coupling M _1 ~ 2.5 x 10 -6 GeV~ as given in 
(15. 6h . it definitely cannot be the standard QCD axion. The naming "axion-like" we use in the paper refers 
to the particle being very light and to its pseudoscalar nature, reflected in the form of the interaction 

dSU). 

Since at present there is no an alternative explanation of the PVLAS data, we are faced to the 
challenge of finding a consistent model that could explain the constraints (|5.4p . (|5.5|) . and (|5.6p in terms 
of a light particle coupled to two photons. The route we have followed has been investigating ways to 
evade the astrophysical bounds. We have worked out two possibilities that could solve the problem. 
The first is that (^-particles are indeed produced in the Sun but that they interact so strongly that are 
trapped by the solar medium. Then, the energy of the emitted (^-particles is much lower than in the usual 
free-streaming regime and thus the CAST telescope is not able to detect them. In Section T5.2I we propose 
a simple model with paraphotons that provides a way (^-particles are trapped. However, we will see that 
it leads to too strong photon-paraphoton interactions that are not consistent with other observations. 
Even having this problem, we present the model because, first of all, it remains to be seen whether a 
sophistication of these ideas may lead to a consistent model. Second, some of the issues we are faced 
are helpful in Section 15.31 The second possibility we examine is that the Primakoff process is suppressed 
when occurring in a stellar medium. Then, there would be far less ^-particles emitted than expected. 
We discuss in Section 15.31 how a composite <fi and the corresponding form factor of the (^77 vertex could 
be responsible for such a suppression of the c/>-flux. Finally, we present our conclusions and additional 
comments in Section [5.41 

5.2 Trapping regime 

Let us start analyzing the possibility that (^-particles are produced in the solar core but that interact so 
strongly with the medium that their fate is analogous to what happens to the stellar photons, namely, 
they abandon the Sun after a lot of interactions, having followed a random walk path. In this trapping 
regime, local thermodynamic equilibrium applies and (ft would contribute to the radiative energy transfer. 
The total opacity, including the exotic contribution, 

Kotal = k -y + hi 
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should not be much different from the standard solar opacity fcT 1 ~ 1 g/cm 2 , if we do not want to ruin 
the standard solar model. Specifically, one imposes [27J 11^1] 

V 1 < i A ( 5 - 7 ) 

The key point of course is to try to implement this possibility within a particle physics model. The 
scenario we shall examine assumes paraphotons provide the trapping interaction. These vector particles 
were proposed by Okun in [26] (see also [122] ) and further developed by Holdom in [12311124] . The basic 
idea is a modification of QED that consists in adding an extra U{\) abelian gauge symmetry. If j„ is the 
electromagnetic current involving charged particles ju, ~ e7^e + ... we start with the lagrangian 

Co = ~F^Fi ia , + eij lt A^ (5.8) 

This lagrangian has a U{\) gauge symmetry group, and would be the photon part of the QED lagrangian. 
The paraphoton model assumes two groups U{l)i x U{\)2 as the gauge symmetry, so that one has two 
gauge fields A\ and A 2 . 

In the line of [1231 1124] we will assume that there are very massive particles carrying charges under 
both U\ and U 2 groups. At low energies, these massive particles running in loops can be integrated out 
leaving the lagrangian 

C = -\{1 + 2 ei i)FfiV - |(1 + 2e 22 )F 2 ^iV 

+\e 12 F^F 2 ^ + ei j„Al (5.9) 

The parameter ei 2 is the induced mixing in the kinetic terms, and en and £22 are also modifications 
due to these loops. At first order in the small e-parameters, we define new fields that diagonalize and 
normalize the kinetic terms, 

A lt = (l + e n )A lll (5.10) 

A'^ = (1 + e 22 )A 2fl - e 12 A lfl 

We end up with the photon A^ coupled to charged particles, with e = ei(l— en), and with the paraphoton 
A' 

C = -\F^F, V - \f'^FI v + ejpA» (5.11) 

Different authors have added some physics to (|5.9p and (|5.1ip so that phenomenological consequences 
arise. In [26] and [122] the effects of a paraphoton mass were discussed. In [12311124"] . it was shown that 
the existence of light particles having U 2 charge leads to these particles having an electric charge of size 
ee. In [125] the paraphoton was identified with a mirror photon and some implications were analyzed. 
The most recent work [126] considers higher-order operators to describe the interaction of the paraphoton 
with matter. 

What we propose is to add to C in (15. 9ft the interaction 



£07272 = jjy-F^F^cf) (5.12) 

with M 2 a low energy scale. The axion-like particle is therefore strongly coupled to the U 2 gauge boson. 
After diagonalizing (I5.10p we get a strong coupling of 4> to paraphotons 



1 

4M 2 ~ 



£ ^V = TW F '^> ( 5 - 13 ) 
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a weaker coupling with a mixed term 



= ^r| FA 



and finally we get a term that couples <j> to photons, i.e., an interaction as in (|5.ip with the identification 

h = % (5 ' 15 > 

with M giving the strength (|5.6p . 

We now have the necessary ingredients to have a large opacity of (ft in the solar medium. The dominant 
contribution to the opacity comes from the process 

07 -► <fr/' (5.16) 

where a virtual 7' is exchanged. The secondary paraphotons are further scattered 

y 7 _^ yy (5.17) 

where now a is exchanged. In both reactions, (|5.16p and (|5.17p . the initial 7 is of course from the stellar 
plasma. 

Having exposed the main idea, we proceed to the calculation of the opacity, where we shall content 
ourselves with order of magnitude estimates. The head-on collision in (|5.16p has a total cross-section 

\ 2 

s (5.18) 



384vr \M$ 

with s = {jptf, +Pj) 2 - The total cross-section for (|5.17p is 

2 



768vr \M$ 



(5.19) 



with s = (p 7 + py) 2 - I n A5.18D and (I5.19J) all particle masses are neglected in front of s. To estimate the 
opacity we set s ~ 4(E" 2 ) and (E?) ~ 10.3 T 2 , where T is the temperature of the medium. We get 



<**> 



5xur '^(s!) (k^) c ™ (5 ' 20 > 



and 

(v> =^='« ">-'<* {B) A (iS)" ™ < 5 - 2i > 

Requiring a large enough opacity, eq. (|5.7|) . for the conditions of the Sun core, T ~ IkeV, p ~ 
100 g cm" 3 , we are lead to 

Ml ~ k^ (5 - 22) 

This condition comes from the reaction f)5. 1T|) : the process (I5.16|) gives a weaker condition. 

In our model, the Sun is a copious emitter of low energy axion-like particles and paraphotons. How- 
ever, there could be no axion-like particles reaching the Earth, because of the decay <j) — ► W- The 
lifetime of cj> with energy E^ ~ 3T cscape is 
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cf) would decay before reaching the Earth when the parameters of our model are such that r^ < 500 s. 
In this case, only paraphotons, from emission or decay, would survive the journey from the Sun to the 
Earth. 

Using (|5.15p and the experimental value (|5.6p and then imposing condition (|5.22p we find the allowed 
values for M 2 and £12- There is a maximum value for the mixing 612 < 5 x 10" 7 , and also a maximum 
value for the scale M2 < 25 eV. 

Let us now discuss the cosmological constraints on the new interactions. In the early universe, 
production of paraphotons proceeds trough the reaction 

77 -> 7V (5.24) 

The interaction rate T has to be compared to the expansion rate H of the universe to see whether the 
process (15.24p is effective. The calculation is similar to the one leading to (|5.20p and (I5.2ip . Assuming a 
matter-dominated universe, we have 

r .. 1 / en ,< (My fry" 



H 200 VIO- 7 / \M 2 J \eVj 

Clearly, for the values of the parameters M2 and ei2 discussed before and for T > 1 eV, T/H > 1 
and thus a cosmic background of paraphotons will be born (when it is radiation that dominates, (J5.25P 
has to be modified, but we reach the same conclusion). Once there is a 7' population, the situation 
is catastrophic since the interaction (|5.16p is only ef 2 -suppressed while (|5.24p is ef 2 -suppressed. As a 
consequence photons and paraphotons would be in equilibrium for T < 1 eV, in contradiction with the 
observation of having a transparent universe for these low temperatures. 

There might be other constraints on 7 — 7' interactions at high energies coming from example from 
photon-photon interactions in accelerators. However, here we should consider the issue that the vertex 
could be subject to form factor effects. We will discuss about this topic in the next Section. 

5.3 Suppression of the solar production 

Let us investigate now a framework where the production in stellar cores is considerably diminished. A 
first thing to notice is that we should look at (15. ip as an effective lagrangian and consequently we should 
not expect it to be valid at arbitrarily high energies. The well studied 7r°77 vertex is similar to (|5.ip and 
it is useful as a guideline. The crucial point is that when one of the photons (or both) is off mass-shell 
the effects of the 7r°-photon transition form factor become manifest. 

There are indeed a variety of measurements where the transition form factor of pseudoscalar mesons 
can be observed, from moderate q 2 up to large momentum transfer [12T^ I128| 1129^ [130] . Let us emphasize 
that the appearance of a form factor is expected on general grounds. From the theoretical point of view, 
apart from the phenomenological VMD parameterization, one gets a form factor when using a quark- 
triangle model [131J. I132J . when calculating in perturbative QCD and when using some other methods 
|133l I134J, I135J, 1136} I137J . All these approaches are consistent among themselves and are able to fit the 
data. For example, when the 7r°77 vertex is described by a quark triangle loop with off-shell photons, 
the explicit calculation of the diagram leads to a form factor that can be identified with VMD provided 
one assigns constituent masses to the internal up and down quarks [131} I132| . Then, for high q 2 one has 
a suppression M 2 /q 2 ~ M 2 Jq 2 . 

These facts have encouraged us to postulate that the axion-like particle ^ is a confined bound-state 
of quark-like particles, that we will call preons in accordance with tradition. If for simplicity we consider 
one fermion / as the only preon, (j) would be the J p = 0~ // bound state and the coupling to two photons 
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would proceed through a triangle loop with / circulating in it. This would result in the appearance of a 
form factor effect at high energies. When both photons are on-shell there is no suppression; these are the 
conditions in the PVLAS experiment and in the detection setup in CAST. However, in the solar medium 
there would be a suppression of the <j> emission rate. 

Let us calculate which is the required suppression F in the Primakoff amplitude for having a consistent 
scenario. If we call M pv i as the value in (15. 6D and M cas t the lower bound in (15. 5h . we should have 

(5.26) 

where in square brackets there is the relevant factor referred to production in the Sun and outside the 
brackets the factor corresponding to detection in CAST. In the lhs we assume there is suppression, while 
in the rhs we assume none because the CAST limit is obtained assuming no form factor suppression in 
the solar production. Introducing numbers we obtain 

\F\ < 2 x 1(T 9 (5.27) 

We now turn our attention to the theoretical prediction for F, that we obtain from the calculation of 
the preon-triangle diagram amplitude. For invariant masses s± and S2 of the photons, and values of the 
masses of <p, m^, and the internal fermion /, Mf, the amplitude F(si, S2,m ( j 1 ; Mf) can be put in terms 
of dilogarithms [13H 1132] . Let us comment that F is in general a complex quantity and also that, as a 
form factor, we normalize it as F(0, 0,m^,; Mf) = 1. 

The values for s\ and S2 in the solar core will be in the keV range. Indeed, in the interior of the 
Sun the Primakoff production is started by a photon of the thermal bath with approximately up ~ 
(0.4 keV) 2 ~ si, with cup the plasma frequency. The virtual photon connecting the vertex to a proton 
(or to any charged particle) is subject to screening effects, as discussed in |138| . These effects amount to 
cut the momenta contributing to the Primakoff effect with the Debye-Hiickle scale k^u , that in the solar 
core is k 2 BH ~ (9keF) 2 ~ si- 

Provided the mass Mf is much less than s\ and s<i, we obtain a strong suppression compatible with 
(|5.27p . With the values of si 2 mentioned above and for ms < 10~ 3 eV (these are the values for which a 
coherent effect in vacuum is expected in the PVLAS setup) we obtain numerically that F satisfies ([5.27(1 
for 

M f <2x 10~ 2 eV (5.28) 

To see a bit more clearly how the suppression arises, we have verified that the exact value for F, in the 
limit I S2 1 S> |si| ^> Mf S> m^ has the behavior 

\F\ ~ 10 2 - , V (5-29) 

\ s 2\ 

Thus, Mf plays the role of the cut-off energy scale of the ^77 vertex form factor. The scale of new 
physics is again a low energy scale. 

Let us comment that, before, we have identified k 2 DH with S2 and that it is an approximation since 
the i-channel carries other momenta. However, |s2 1 ^ ^nw always, so that, at the view of (|5.29p . the 
approximation is conservative. 

There is a parameter that we have not discussed, namely the electric charge qfe of the preon /. With 
the coupling 

^fl^fd^ (5.30) 
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the result of the calculation of the triangle for on-shell photons is 

qla 1 l 

irv M~4xl0 5 GeV K } 

where we have already identified the result with the coupling 1/M in (|5.ip and with the experimental 
value in ()5.6|) . Also, we have introduce a factor N corresponding to having a confining SU(N) gauge 
group in the preon sector. 

Now we have to take into account the bound on light millicharged particles [831 [80] coming from BBN 
constraints, 

q f <2x 1(T 9 (5.32) 

In a paraphoton model such a small electric charge could arise naturally [1231 1 124] . Together with (|5.3ip 
it implies again low energy scales 

v<W~ 5 eV (5.33) 

(We have used N = 3). 

Let us point out that the SU(N) gauge group should be totally independent of the color SU(3) C 
standard gauge group. Otherwise, among other undesired consequences, we would have in nature hadrons 
with charges near ±(2/3)e, ±(l/3)e, and ±(4/3)e that would form when binding a preon or antipreon 
with quarks or antiquarks, of the kind uf, uuf, etc. Also, at tree level we should take / as a singlet 
under the standard model, with the small electric charge arising as a higher order effect. With all these 
assumptions, we think the new force and particles could have been not noticed in other experiments. 
Yet, there are consequences, like the existence of bound states with higher spin, as for example a state 
with J = 1 that would be unstable since it would decay into (fry. Of course, from a phenomenological 
perspective, it would be interesting to look for signals of the preon model. From a more theoretical point 
of view, here we do not attempt to build a full model, for we think it would be premature. Rather, we 
have shown a possible way to evade the astrophysical limits on axion-like particles. 

5.4 Conclusions 

A recent review by Raffelt |59] emphasizes that the PVLAS result (|5.6p . interpreted in terms of a new 
light axion-particle coupled to photons, would lead to the Sun burning much faster that what is actually 
observed. There is the pressing issue of explaining the results in a consistent model. Of course an 
independent check of the results would be most welcomed; in fact there are interesting proposals for such 
type of laboratory experiment where a high sensitivity would be reached [139] . 

In this paper, we report the work we have done trying to evade the astrophysical bounds and thus 
accommodating a light particle coupled to two photons. The astrophysical limits assume 1) a flux 
calculated with the interaction (|5.1|5.2p in the stellar core, and 2) that the produced particles escape the 
star without further interaction. 

Our first attempt has been trying to find a model where 2) is not true. In our paraphoton model with 
(I5.22J) . axion-like particles 4> are trapped in the stellar interior and so are the paraphotons 7' produced 
by cftj scattering. The large opacity makes 4> production not a problem. It follows that the astrophysical 
limit (]5.4p . that assumes <fi freely escapes, is no longer valid. These arguments have nothing to say about 
an axion-like interpretation of the PVLAS result, because it is an Earthbound laboratory experiment, 
with cj) produced and detected in the laboratory. 

The model, however, leads to photons to interact with paraphotons so strongly that it is excluded, 
at least in the simple framework we have exposed where a cosmological background density of parapho- 
tons emerges. Perhaps a more elaborated model with paraphotons, or another model with a different 
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strong interaction can do the job of trapping particles in the Sun without entering in conflict with other 
experiments. Let us point out that there are also astrophysical constraints on the ^77 coupling coming 
from red giants and from SN observations. As far as the SN is concerned, the low value (|5.6p makes 
(p to be trapped in the SN core in such a way that one does not need extra interactions [62]. In any 
case, the exercise we have presented in Section 15.31 shows that it is not trivial to evade the astrophysical 
constraints. 

The astrophysical bound could also be evaded if the ^77 vertex, while fully operating at PVLAS 
energies, is suppressed in the conditions of stellar interiors. In this case the condition 1) above does not 
hold. We have explored the possibility that <p is a composite particle and has a form factor leading to a 
suppression of the production. We have shown that this scenario is able to explain the puzzle. Our ideas 
are highly speculative since they involve preons with a new confining force and probably a minuscule 
electric charge, but notice that we have been inspired by the pion and the 7r°77 vertex, that after all 
have the nice property of being real. In any case, it would be crucial to look for other phenomenological 
consequences of the preon model. We have not found any that rules out it obviously. 

Either in the case of a strong interaction leading to trapping or of a suppression of the production, 
the astrophysical bounds on axion-like particles could be evaded. If indeed they are evaded, there are 
drastic consequences for CAST, since then the non-observation of X-rays does not imply a limit such as 
(I53|). 

Our main conclusion is that the explanation of the PVLAS result in terms of a light particle coupled 
to photons is not necessarily in contradiction with other experiments and observations. Let us emphasize 
that taking alone the PVLAS data, if interpreted in terms of new light particle coupled to photons, it 
would already mean an interesting piece of new physics. But there is even more. The result, taken 
together with the astrophysical limits and the CAST data, means that we have to go beyond the "mere" 
existence of a new pseudoscalar (j) coupled to photons and even more exotic physics has to be invoked. 
In the scenarios of the sort we have proposed the new physics scale is at very low energies. If confirmed, 
it would be an exciting discovery. Otherwise, if finally the models that try to evade the astrophysical 
constraints are shown not to be valid, the situation will be no less exciting since an alternative explanation 
for the PVLAS laser experiment result will be needed. 
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results 
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6.1 Introduction 



Very recently, the PVLAS collaboration has announced the observation of a rotation of the plane of 
polarization of laser light propagating in a magnetic field [33]. This dichroism of vacuum in magnetic 
fields may be explained as the oscillation of photons into very light particles <j>. If true, this would be of 
course a revolutionary finding |140| . 

The lagrangian that would describe the necessary (^77 coupling is 

C <hv = W e ' u "" F '" FPff<t> (61) 

when cj) is a pseudoscalar, and when it is a scalar is 

c ^ = m F " VF ^^ (6 - 2) 

with F^ v the electromagnetic field tensor. We shall refer to eft in both cases as an axion-like particle 
(ALP). Let us remark that a transition to a spin-two particle contributes to the polarization rotation 
negligibly [53] , 

Either (|6.ip or (|6.2p lead to 7 — (f> mixing in a magnetic field and, if (f> is light enough, to coherent 
transitions that enhance the signal [33]. Interpreted in these terms, the PVLAS observation [13] leads to 
a mass for the ALP 

1 meV < rritj, < 1.5 meV (6-3) 

and to a coupling strength corresponding to 

2 x 10 5 GeV < M < 6 x 10 5 GeV . (6.4) 

Of course we would like to have an independent test of such an interpretation. There are ongoing projects 
that will in the near future probe 7 — <j> transitions [141j . In the meanwhile we should face the problem 
of the apparent inconsistency between the value (|6.4p and other independent results, namely, the CAST 
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observations [79] on the one hand, and the astrophysical bounds on the coupling of ALPs to photons on 
the other hand [55] . 

The CAST collaboration has recently published [79] a limit on the strength of (|6.ip or (|6.2p , A light 
particle coupled to two photons would be produced by Primakoff-like processes in the solar core. CAST 
is an helioscope |75] that tries to detect the <j> flux coining from the Sun, by way of the coherent transition 
of (j)'s to X-rays in a magnetic field. As no signal is observed they set the bound 

M > 0.87 x 10 10 GeV , (6.5) 



which is in strong disagreement with (|6.4j) . 

Also, the production of <^'s in stars is constrained because too much energy loss in exotic channels 
would lead to drastic changes in the timescales of stellar evolution. Empirical observations of globular 
clusters place a bound [65], again in contradiction with (|6.4p . 



M > 1.7 x 10 10 GeV . (6.6) 

As it has been stressed in [1], once we are able to relax (|6.6p we could also evade (|6.5p . Indeed, the 
CAST bound assumes standard solar emission. From the moment we alter the standard scenario we 
should revise (|6.5p . In [HEHES] two ideas on how to evade the astrophysical bound (|6.6p are presented. 
One possibility is that the produced ALPs diffuse in the stellar medium, so that they are emitted with 
much less energy than originally produced pQ. A second possibility is that the production of ALPs is 
much less than expected because there is a mechanism of suppression that acts in the stellar conditions. 
We will present in this letter a paraphoton model with a low energy scale where the particle production 
in stars is suppressed enough to accommodate both the CAST and the PVLAS results. 

6.2 Triangle diagram and epsilon-charged particles 

The physical idea beyond this letter is that to understand PVLAS and CAST in an ALP framework we 
have to add some new physics structure to the vertices (|6.ip . (|6.2p . The scale of the new physics should 
be much less than O(keV), the typical temperature in astrophysical environments. 

We will assume that this structure is a simple loop where a new fermion / circulates; see Fig. (l7.3p . 
The amplitude of the cp'yy diagram can be easily calculated and identified with the coefficient in (|6.1|) or 
(1531) 

M IT V 

Here a = e 2 /An, and the charge of the fermion / is eqj. The value of the mass-scale v depends on 
the (j)ff vertex. If ^ is a pseudoscalar vps = m f/9PS while if is scalar vs = ffafj-m^) not far from 
vs ~ mt ~ m^ if rrif ~ m^. Finally if eft is & Goldstone boson vqb is related to the scale of breaking of 
the related global symmetry. 

From (16. 7p we see that M, the high energy scale (16. 4p . is connected to v. As we need v to be a low 
energy scale, qj should be quite small. 

Paraphoton models [1231 1124] naturally incorporate small charges. These models are QED extensions 
with extra U(l) gauge bosons. A small mixing among the kinetic terms of the gauge bosons leads to 
the exciting possibility that paracharged exotic particles end up with a small induced electric charge 
[123] H24]. 

Getting a small charge for / is not enough for our purpose since we need also production suppres- 
sion of exotic particles in stellar plasmas. With this objective, we will present a model containing two 
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paraphotons; if we allow for one of the paraphotons to have a mass, we will see we can evade the astro- 
physical constraints and consequently the model will be able to accommodate all experimental results. 
We describe it in what follows. 



6.3 A model with two paraphotons 

Let us start with the photon part of the QED lagrangian 
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where Jq^ is the electromagnetic current involving electrons, etc., Jq^ ~ e^/^e + .... From the Uq(1) gauge 
symmetry group, we give the step of assuming Uq(1) x U\(l) x [^(l) as the gauge symmetry group, with 
the corresponding gauge fields Aq, Ai, and A%. With all generality there will be off-diagonal kinetic 
terms in the lagrangian, like eoi FqF\ and €02-^0-^2 (Lorenz index contraction is understood). We expect 
these mixings to be small if we follow the idea in [12311124] that ultramassive particles with 0,1,2 charges 
running in loops are the responsibles for them. We will assume that these heavy particles are degenerate 
in mass and have identical 1 and 2 charges so that they induce identical mixings eoi = £02 



f. 



To write the complete lagrangian we use the matrix notation A = (Aq, Ai, A2) t and F = (Fq,F\,F2 



C = -±F T M F F + lA T M A A + eY,3iAi 



(6.9) 



We call Aq, Ai, and A2 interaction fields because the interaction term in (|6.9p is diagonal, i.e. the 
interaction photon is defined to couple directly only to standard model particles. Here the kinetic matrix 
contains the mixings, 
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(6.10) 



In general the diagonal terms are renormalized, 1 — > 1 + 5, and there are terms Mfi2- However, they do 
not play any relevant role here and we omit them. 

As said, we need one of the paraphotons to be massive but it will prove convenient to work with a 
general Ma — D««5 { m oi m i > m i}- Also, in the last term of (16. 9p we see the currents j\ and j'2 containing 
the paracharged exotic particles. To reduce the number of parameters we have set the unit paracharge 
equal to the unit of electric charge, so that there is a common factor e. 
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Figure 6.1: Triangle diagram for the (pjj vertex. 



12 



Compatibility of CAST search with axion-like interpretation of PVLAS results 



Diagonalization involves first a non-unitary reabsorption of the e terms in (|6.10[) to have the kinetic 
part in the lagrangian in the canonical form, (— 1/A)F T F. After this, we diagonalize the mass matrix 
with a unitary transformation that maintains the kinetic part canonical ending up with the propagating 
field basis A. We have A = UA, with 
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(6.11) 



We see that the interacting and the propagating photon differ by little admixtures of 0(e). (We work at 
first order in e). 

We have developed a quite general two paraphoton model. The specific model we adopt has the 
following characteristics. First, only one paraphoton has a mass, say mi = \x ^ 0, and TO2 = 0. Second, 
in order to get the effects we desire we have to assign opposite 1 and 2 paracharges to /, so that the 
interaction for / appearing in the last term of (16. 9p is 



ef^f(Al-A%) (6.12) 

Let us show why we choose these properties. The coupling of / to photons in the interaction basis 
is shown in Fig. (|7.2p . It proceeds through both paraphotons, with a relative minus sign among the two 
diagrams due to the assignment (j6.12|) . The induced electric / charge is thus 

q f = U10-U20 (6.13) 

We see from (|6.1ip that TO2 = implies U20 = — e. However, the value for U\q has to be discussed 
separately in the vacuum and the plasma cases. In vacuum, we have too = 0, so that U\q = and thus 
qf = e. In this case A$ and Ai are degenerate and we can make arbitrary rotations in their sector. 
This corresponds to different charge assignments that of course leave the physics unchanged. Due to our 
method of handling the diagonalizations, eq. (|6.1ip is bad behaved for too = 77*2, except for the case of 
our interest, too = TO2 = 0, in which the order in which we take the limits too — ► and TO2 — ► gives 
different charge assignments according to the rotational freedom. Here we have made TO2 — > before 
too — > to provide f with a millielectric charge as in [1231 1124| . Changing the order of the limits would 
end with a paracharge to electrons. 

In the classical and non-degenerated plasmas we consider the dispersion relation can be taken as 
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Figure 6.2: Diagrams of the interaction of f with photons. 
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Provided we have a low energy scale \i <C top ~ keV, we reach our objective of having a strong 
decrease of the / charge in the plasma, i.e., qf(uj p ) <C 9/(0) = e. 

The cancelation of the two diagrams of Fig l7. 21 requires that the equality e\ = e2 holds up to terms of 
order 0(fj, 2 /up). Note that even if e\ = e2 at some high energy scale because of a symmetry, a difference 
in the beta functions could also spoil our mechanism at low energy. The parafermion / contributes 
equally to both beta functions so the problem are the contributions from the sector that gives mass only 
to A\. However these contributions can be made arbitrarily small by sending the Higgs mass to infinity 
in the spirit of the non-linear realizations of symmetry breaking, by considering Higgsless models like 
breaking the symmetry geometrically, or by considering gauge coupling unification e\ = e2 at an energy 
not far from the typical solar temperature. A further possibility is to consider e\ = e2 <C e which would 
suppress the loop-induced effects at the prize of making the model less natural. 

6.4 The role of the low-energy scale 

We now discuss the consequences of our model. The PVLAS experiment is in vacuum, so / has an 
effective electric charge (?/(0) = e, which from (16. 7p has to be 

e 2 ~10 _12 — (6.15) 

eV 

Concerning the astrophysical constraints, we notice that the amplitude for the Primakoff effect ^Z —* 
(j)Z is of order qj = e 2 and that there are production processes with amplitudes of order e which will be 
more effective. One is plasmon decay 7* — > ff. Energy loss arguments in horizontal-branch (HB) stars 
[831 [80] limits qj to be below 2 x 10~ 14 , which translates in our model into the bound 

„2 

e-^<4x 1(T 8 (6.16) 

(we have used wp ~ 2 keV in a typical HB core). Other processes like bremsstrahlung of paraphotons 
give weaker constraints. 

Equations (|6.15j) and (|6.16p do not fully determine the parameters of our model. Together they imply 
the constraint 

vfi 4 < {OAeVf (6.17) 

We can now make explicit one of our main results. In the reasonable case that v and \i are not too 
different, we wee that the new physics scale is in the sub eV range. 

Let us consider now the CAST limit. The CAST helioscope looks for 0's with energies within a window 
of 1-15 keV. In our model, /'s and paraphotons are emitted from the Sun, but we should watch out <f> 
production. This depends on the specific characteristics of (p. We consider three possibilities. A) cp is 
a fundamental particle. As we said the Primakoff production is very much suppressed, so production 
takes place mainly through plasmon decay 7* — > ff(fi- The c^-flux is suppressed, but, most importantly, 
the average (ft energy is much less than top ~ .3 keV, the solar plasmon mass. The spectrum then will 
be below the present CAST energy window. B) cj) is a composite // particle confined by new strong 
confining forces. The final products of plasmon decay would be a cascade of 0's and other resonances 
which again would not have enough energy to be detected by CAST. C) cj) is a positronium-like bound 
state of //, with paraphotons providing the necessary binding force. As the binding energy is necessarily 
small, ALPs are not produced in the solar plasma. 

Let us now turn our attention to other constraints. Laboratory bounds on epsilon-charged particles 
are much milder than the astrophysical limits, as shown in [83|. 180] . In our model, however, even though 
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Figure 6.3: Constraints on the parameters of our model. The black area is excluded by Cavendish-type 
experiments, and the grey area by the astrophysical constraint A6.15]) . The dashed line corresponds to 
v = /i, and the dot to v = ji ~ 1 meV. 



paraphotons do no couple to bulk electrically neutral matter, a massive paraphoton A\ couples to electrons 
with a strength e and a range // . This potential effect is limited by Cavendish type experiments 

mams]. 

In Fig. (|6.3p we show these limits, as well as the astrophysical bound (|6.16j) . In the ordinates we can 
see both e and v, since we assume they are related by fl6. 15f) . At the view of the figure, we find out that 
there is wide room for the parameters of our model. However we would like that v and \x do not differ 
too much among them. We display the line v = fj,, the region where this kind of naturality condition is 
fulfilled. The most economical version of the model would be obtained when the new scales are, on the 
order of magnitude, about the scale of the ALP mass of 0(1 meV), (|6.3j) . We have also indicated this 
privileged point in the parameter space. 

Also, we should discuss cosmological constraints, i.e. production of paraphotons and /'s in the early 
universe. Taking into account that the vertices have suppression factors in the high temperatures of such 
environment, we find that there is not a relic density of any of them. 

Finally, let us come back to the physics responsible for the A\ mass. If this comes from an abelian 
Higgs mechanism then the Higgs boson acquires a millicharge ee\ and could be produced in the Sun and 
in the early universe, particularly in the period of primordial nucleosynthesis. However, this is not a 
problem if the mass of the Higgs is large enough, a constraint that we required at the end of section 16.2 
when discussing charge running. 



6.5 Conclusions 



We have presented a model of new physics containing a paracharged particle / and two paraphotons, 
one of which has a mass \i that sets the low energy scale of the model. With convenient assignments of 
the / paracharges and mixings, we get an induced epsilon-charge for / that moreover decreases sharply 
in a plasma with wp 3> fi. Our model accommodates an axion-like particle with the properties (16.3f) and 
(I6.4p . able to explain the PVLAS results, while at the same time consistent with the astrophysical and 
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the laboratory constraints, including the limit obtained by CAST. 

We have some freedom in the parameter space of our model; however if we wish that the energy 
scales appearing in it are not too different, we are led to scales in the sub eV range. A preferred scale is 
O(meV), because then it is on the same order than the axion-like particle mass. 

If the interpretation of the PVLAS experiment is confirmed, which means the exciting discovery of an 
axion-like particle, then to make it compatible with the CAST results and with the astrophysical bounds 
requires further new physics. In our model, the scale of this new physics is below the eV. 

Note added: Recently, a paper has appeared [143] that justifies our model in the context of string 
theory. 

6.6 acknowledgments 

We acknowledge support by the projects FPA2005-05904 (CICYT) and 2005SGR00916 (DURSI). 



Chapter 7 
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Published in Phys. Rev. D 75, 013004 (2007) [3]) 



7.1 Introduction 



Recently the PVLAS collaboration has reported the observation of a rotation of the polarization plane 
of a laser propagating through a transverse magnetic field [33]. This signal could be explained by the 
existence of a new light neutral spin zero boson (j>, with a coupling to two photons |33|. [50] 



4^ = ^^^^ or 4 +) = ^ (+) ^^ (7-1) 

depending on the parity of <f>, related to the sign of the rotation which up to now has not been reported_|. 
Such an Axion-Like Particle (ALP) would oscillate into photons and vice versa in the presence of an 
electromagnetic field in a similar fashion as the different neutrino flavors oscillate between themselves 
while propagating in vacuum. 

The PVLAS signal, combined with the previous bounds from the absence of a signal in the BFRT 
collaboration experiment |40j . implies [43] 



1 meV < m < 1.5 meV, 2 x 10 5 GeV < M < 6 x 10 5 GeV, (7.2) 

with m the mass of the new scalar. 

It has been widely noticed that the interaction (|7.ip with the strength (J7.2I) is in serious conflict 
with astrophysical constraints [591 1144] , while it is allowed by current laboratory and accelerator data 
|61|. 186] . This has motivated recent work on building models that evade the astrophysical constraints 
HI EU [Ml El [Ul], as well as alternative explanations to the ALP hypothesis [1351 |27]|l33"]. 

At the same time, many purely laboratory-based experiments have been proposed or are already on the 
way to check the particle interpretation of the PVLAS signal [1391 11351 H371 11351 [1391 [1501 1431 11511 1571 139] . 
It is important to notice, for the purpose of our paper, that these experiments are optical, and not high- 
energy, accelerator experiments. 

Quite generally, these experiments will have enough sensitivity to check values of M equal or greater 
than 10 6 GeV, but, apart from Ref. [139] . they do not have the impressive reach of the astrophysical 



1 The PVLAS collaboration has also found hints for an ellipticity signal. The sign of the phase shift suggests an even 
particle <f> [56] . 
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Figure 7.1: Primakoff processes in which a photon turns into an ALP in the electric field of a charged 
particle like a proton or electron. 

considerations, implying M > 10 10 GeV. Thus, if the PVLAS signal is due to effects other than <fi — 7 
oscillations and the astrophysical bounds are applicable, these experiments can not detect any interesting 
signal. 

However, the astrophysical bounds rely on the assumption that the vertex (|7.1h applies under typical 
laboratory conditions as well as in the stellar plasmas that concern the astrophysical bounds. It is clear 
that, if one of the future dedicated laboratory experiments eventually sees a positive signal, this can not 
be the case. 

In this work we investigate the simplest modification to the standard picture able to accommodate 
a positive signal in any of the forthcoming laboratory experiments looking for ALPs, namely that the 
structure of the interaction (|7.ip remains the same in both environments, while the values of M and 
m can be different. Interestingly enough, the environmental conditions of stellar plasmas and of typical 
laboratory experiments are very different and thus one could expect a very big impact on M and m. 

We consider qualitatively the situation in which the dependence of M and m on the environmental 
parameters produces a suppression of ALP production in stellar plasmas. The main work of the paper 
is devoted to compute this suppression using a realistic solar model and to investigate how it relaxes the 
astrophysical bounds on the coupling (|7.ip . This leaves room for the proposed laboratory experiments 
to potentially discover such an axion-like particle. 

In section [?T2l we revisit the astrophysical bounds and discuss general mechanisms to evade them. In 
the following section [731 we present our scenario of environmental suppression and calculate the modified 
bounds. We present our conclusions and comment on the reach of proposed future laboratory experiments 
in section 17.41 

7.2 Astrophysical bounds and general mechanisms to evade them 

Presuming the 1^77 vertex (|7.ip , photons of stellar plasmas can convert into ALPs in the electromagnetic 
field of electrons, protons and heavy ions by the Primakoff effect, depicted schematically in Fig. 17.11 
If M is large enough, these particles escape from the star without further interactions constituting a 
non-standard energy-loss channel. This energy-loss channel accelerates the consumption of nuclear fuel 
and thus shortens the duration of the different stages of stellar evolution with respect to the standard 
evolution in which ALPs do not exist. 

In general, the astrophysical observations do agree with the theoretical predictions without additional 
energy-loss channels so one is able to put bounds on the interaction scale M [65] . The most important 
for our work are those coming from the lifetime of the Sun |152j . the duration of the red giant phase, and 
the population of Helium Burning (HB) stars in globular clusters |138[I153] . The last of them turns out 
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Env. param. 


Solar Core 


HB Core 


PVLAS 


T [keV] 


1.3 


8.6 


~0 


q 2 [keV 2 ] 


~ 1 


~ 1 


~ lO" 18 


up [keV] 


0.3 


2 





p [g cm -3 ] 


1.5 x 10 2 


10 4 


< 10~ 5 



Table 7.1: Comparison between the values of environmental parameters, such as the temperature T, 
typical momentum transfer q, plasma frequency Up, and matter energy density p, in the stellar plasma 
and in the PVLAS experiment. Other parameters to consider could be the Debye screening scale k s , or, 
to name something more exotic, the neutrino flux, or the average electromagnetic field. 



to be the most stringent, implying 



,10 



M > 1.7 x 10 1U GeV = M HB , 



(7.3) 



for m < 0(1 kel^). Moreover, if ALPs are emitted from the Sun one may try to reconvert them to photons 
at Earth by the inverse Primakoff effect exploiting a strong magnetic field. This is the helioscope idea 
|75j that it is already in its third generation of experiments. Recently, the CERN Axion Solar Telescope 
(CAST) collaboration has published their exclusion limits [79] from the absence of a positive signal, 



M > 8.6 x 10 9 GeV = M C ast, 



(7.4) 



for m < 0.02 eV. 

One should be aware that these astrophysical bounds rely on many assumptions to calculate the flux 
of ALPs produced in the plasma. In particular, it has been assumed widely in the literature that the same 
value of the coupling constant that describes <f> — 7 oscillations in a magnetic field in vacuum describes 
the Primakoff production in stellar plasmas, and the mass has been also assumed to be the same. We 
want to remark that this has been mainly an argument of pure simplicity. In fact, there are models in 
which M depends on the momentum transfer q at which the vertex is probed pQ or on the effective mass 
ujp of the plasma photons involved [2]. These models have been built with the motivation of evading the 
astrophysical bounds on ALPs, by decreasing the effective value of the coupling 1/M in stellar plasmas in 
order to solve the inconsistency between the ALP interpretation of PVLAS and the astrophysical bounds. 
This has proven to be a very difficult task because of the extreme difference between the PVLAS value 
(I7.2p and the HB (1T.3() or CAST (J7.4D exclusion limits. These models require very specific and somehow 
unattractive features like the presence of new confining forces or tuned cancelations (note, however, 
|143| ). Anyway, they serve as examples of how M (and eventually m) can depend on "environmental" 
parameters n = q,ujp, etc... (for other suitable parameters, see Table I7TT1) . 



M — > M(rj), m^m(ri), 



(7.5) 



such that the production of ALPs is suppressed in the stellar environment. 

In the following, we will not try to construct micro-physical explanations for this dependence but 
rather write down simple effective models and fix their parameters in order to be consistent with the 
solar bounds and PVLAS or any of the proposed laboratory experiments. 

A suppression of the production in a stellar plasma could be realized in two simple ways: 

(i) either the coupling 1/M decreases (dynamical suppression) or 



so 
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Figure 7.2: Environmental parameters as a function of the distance to the solar center. Temperature 
(solid, red), matter density (dashed, blue), Debye screening scale (double dashed, green) and plasma 
frequency (triple dashed, black), normalized to their values in the solar center, To = 1.35 keV, po = 
1.5 x 10 2 gcm~ 3 , k s0 = 9 keV, uj po = 0.3 keV for the solar model BS05(OP) of Bahcall et al. [153] . 



(ii) m increases to a value higher than the temperature such that the production is Boltzmann sup- 
pressed (kinematical suppression). 

All the environmental parameters considered in this paper are much higher in the Sun than in labora- 
tory conditions (see Table [7TT1 and Fig. 17. 2p . so we shall consider M{rj) and m(rj) as monotonic increasing 
functions of rj with the values of M(~ 0) and m(~ 0) fixed by the laboratory experiments. 

Clearly, both mechanisms are efficient at suppressing the production of ALPs in the Sun, but there 
is a crucial difference that results in some prejudice against mechanism (ii). Mechanism (i) works by 
making the already weak interaction between ALPs and the photons even weaker. The second mechanism, 
however, is in fact a strong interaction between the ALPs and ordinary matter, thereby making it difficult 
to implement without producing unwanted side effects. We will nevertheless include mechanism (ii) in 
our study, but one should always keep this caveat in mind. 

As we said, rj in the stellar plasma is generally much higher than in laboratory-based experiments. 
It is then possible that new ALP physics produces also a big difference between the values of the ALP 
parameters, m and M, in such different environments. 

Let us remark on the a priori unknown shape of M(rj) and m{rf). In our calculations we use a 
simple step function (cf. Fig. IT.3H . which has only one free parameter: the value for the environmental 
parameter where the production is switched off, ?7crit- In most situations this will give the strongest 
possible suppression. The scale ?7 cr it can be associated with the scale of new physics responsible for the 
suppression. In what follows, we will consider only the effects of one environmental parameter at once 
although it is trivial to implement this framework for a set of parameters. 

For simplicity, we restrict the study of the environmental suppression of ALPs to our Sun because 
we know it quantitatively much better than any other stellar environment. The group of Bahcall has 
specialized in the computation of detailed solar models which provide all the necessary ingredients to 
compute accurately the Primakoff emission. We have used the newest model, BS05(OP) [154] . for all 
the calculations of this work (our accuracy goal is roughly 10%). The variation of some environmental 
parameters is displayed in Fig. 17.21 as a function of the distance from the solar center. 
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Figure 7.3: Coupling as a function of an environmental parameter n: The simple form used in our 
calculations (solid line) and a generic, more realistic, dependence (dashed). 

7.3 Numerical results 

Let us first state how a suppression S of the flux of ALPs affects the bounds arising from energy loss 
considerations and helioscope experiments. If the flux of ALPs from a stellar plasma is suppressed by a 
factor S, the energy loss bounds on M are relaxed by a factor of \/S while the CAST bound relaxes with 

Vs, 

M\ oss - > yS M\ oss , energy loss bound, (7.6) 

Mcast -> v^Mcast, CAST bound, (7.7) 

since the former depends only on the Primakoff production, ~ 1/M 2 , and the latter gets an additional 
factor ~ l/M 2 for the reconversion at Earth resulting in a total counting rate ~ 1/M 4 . 
According to (J7.2J) and (J7.4J) to reconcile the CAST and PVLAS results we would need 



/Mpvlas\ 



l 



SCAST < ( J * ) ~ 10 . ( 7 - 8 ) 

while to reconcile the PVLAS ALP with the Sun energy loss bound we need a much more moderate p 

7.3.1 Dynamical suppression 

We consider first a possible variation of the coupling that we have enumerated as mechanism (i). Treating 
the emission of ALPs as a small perturbation of the standard solar model, we can compute the emission of 
these particles from the unperturbed solar data. The 7 — 4> Primakoff transition rate for both interactions 
in eq. (|7.ip can be written as (neglecting the plasma mass lop for the moment jf| 

t,/ N Tk S f 1 , n 1 + COS0 1 , n , 

r wL_,f, = — - — dcosu—z -Z -TT-, xo 1 (7.10) 

v n 64vr y_! k 2 + 1 - cos 6 M(n) 2 ' v ; 

where u is the energy of the incoming photon, and 

^ 2 = ^K + £^)> (7.11) 



2 We have recalculated the Solar energy loss bound using the latest Solar model BS05(OP). See next section for details. 
3 We are using natural units H = c = 1 with the Boltzmann constant, fcs = 1. 



Need for purely laboratory-based axion-like particle searches 



' 


CO 




ON 


o 


a 


o 


u 


^H 


> 




CD 




^ 










a; [keV] 

Figure 7.4: Our spectrum of ALPs at Earth (black solid) agrees reasonably well with that of the CAST 
collaboration [79] (dashed orange) for M = 10 10 GeV. 



is the Debye screening scale, rii, Z^ are the number densities and charges of the different charged species of 
the plasma, a ~ 1/137, n e is the electron number density, cos# is the relative angle between the incoming 
photon and the outgoing ALP in the target frame (we take the mass of the target to be infinite since the 
masses of protons and electrons far exceed the typical momentum transfer of order < keV; this implies 
that the incoming photon and the outgoing axion have the same energy) and k 2 = k 2 /2ui 2 . Integration 
over the whole Sun with the appropriate Bose-Einstein factors for the number density of photons gives 
the spectrum of ALPs (number of emitted ALPs per unit time per energy interval), 



d 2 N{u) 
dujdt 



Air 



R 



- 2 u^rv)^ 

K aii—^ — -^ 



(7.12) 



(Remember that T, k 2 , etc. depend implicitly on the distance R from the solar center.) 

As a check of our numerical computation we have computed the flux of standard ALPs at Earth which 

is shown in Fig. 17.41 and does agree with the CAST calculations |79j . 
It is very important to differentiate two possibilities: 

1. r] is a macroscopic (averaged) environmental parameter given by the solar model and depending 
only on the distance R from the solar center. Then the suppression acts as a step function in the 
R integration (|7.12p for the flux. 

2. r] depends on the microscopic aspects of the production like the momentum transfer q 2 . Then the 
step function acts inside the integral in eq. (|7.1U|) . 

We now start with the first possibility and let the second, which requires a different treatment, for 
subsubsection 17.3.11 



Dynamical suppression from macroscopic environmental parameters 

If 1/M{rj) is a step function, ALP production is switched off wherever rj > rj cr i t . Let us call R cr u the 
radius at which the coupling turns off, i.e. rj(R cr i t ) = r] CT i t . Since the functions r](r) shown in Fig. 17.21 are 
monotonous, we can calculate the suppression as a function of R cr u and then determine r\ cr it = r\[R cr it)- 
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Figure 7.5: Suppression of the flux of ALPs S(ujq = lkeV,i? cr it) as a function of R cr it- 



We define the suppression efficiency, S(to, R cr it), as the ratio of the flux of ALPs with energy uj with 
suppression, divided by the one without suppression, 

d*N(u;;R crit ) f d^N^ y 1 
S(u>; Rcrit) = 1— 7- —j— 77- • (7.13) 



dujdt 



doj dt 



The CAST experiment is only sensitive to ALPs in the range of (1 — 14) keV. Hence, we must suppress 
the production of ALPs only in this energy range. In order to provide a simple yet conservative bound 
we use the factor S(ujq, R C rit) evaluated at the energy IkeV < ujq < 14keV which maximizes S. We have 
checked that, in all cases of practical interest, ujq is the CAST lower threshold, 1 keV. In Fig. 17.51 we plot 
S (IkeV, Rcrit). In Tab. 17.21 we give some values for 5 together with the corresponding values of 7/ cr i t . 

Looking at Table 17.21 we find that it is possible to achieve the suppression required in Eq. (|T.8j) and 
reconcile PVLAS and CAST, but the critical environmental parameters are quite small; for example, the 
critical plasma frequency is in the eV range. Moreover, the results are sensitive to the region close to the 
surface of the Sun where log(S) changes very fast and our calculation becomes somewhat less reliable. 

We now take a look at the solar energy loss bound (|7.6p . The age of the Sun is known to be 
around 5.6 billion years from radiological studies of radioactive crystals in the solar system (see the 
dedicated Appendix in |155j ). Solar models are indeed built to reproduce this quantity (among others, 
like today's solar luminosity, solar radius, etc.), so one might think that a model with ALP emission 
can be constructed as well to reproduce this lifetime. However, this seems not to be the case for large 
ALP luminosity [153] and it is concluded that the exotic contribution cannot exceed the standard solar 
luminosity in photons. For our purposes this means 



L ALP <L Q = 3.846 x 10 ZD W ~ 1.60 x l(f u eV 



26 



-.30 



with 



f°° , d?N 
Lalp = dioto- 
Jo 



dtodt 



We have computed the ALP emission in BS05(OP), 



L ALP = 1.8 x 10" 



10 iU GeV \ 
M J 



L, 



0- 



(7.14) 



(7.15) 



(7.16) 
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Rcrit/R® 


T crit [keV] 


p cr it [g cm -3 ] 


wp, CT it [keF] 


s 





1.35 


150 


0.3 
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0.2 


0.81 


35 


0.16 


0.67 


0.5 


0.34 
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0.03 


0.08 
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0.2 


0.2 


0.01 


2 x 10~ 3 


0.8 


0.12 


0.09 


0.008 


2 x 10~ 5 


0.85 


0.08 


0.05 


0.006 


2 x 10~ 7 


0.9 


0.05 


0.03 


0.004 


4 x 10~ n 


0.95 


0.025 


0.009 


0.0025 


~ io- 20 



Table 7.2: Several values of S(u>o = 1 keV, -R C rit) with their respective values of the suppression scales 
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Figure 7.6: Suppression S of the total flux of ALPs as a function of the critical radius -R C rit- 

This value is slightly bigger than that of Ref. [7S], which relies on an older solar model |156j . probably 
as a consequence of the different data. 

For the total flux, we find a suppression 



S(R 



crit j 



-^ALp(^crit) 

£alp 



(7.17) 



which we plot in Fig. 17.61 

Remember that in order to avoid a conflict between the PVLAS result and the energy loss argument 
we required Eq. (|7.9p . Si ss < 10~ 10 . Looking at Fig. 17.61 we find that this bound alone requires values 
for the critical environmental parameters that are larger (and therefore less restrictive) than those from 
the CAST bound. 
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Dynamical suppression from microscopic parameters: q 2 

In the previous subsection, we have considered macroscopic environmental parameters like, e.g., the 
temperature T. However, suppression could also result from a dependence on microscopic parameters 
like, e.g., the momentum transfer q 2 in a scattering event (not averaged). 

In this section we discuss the well motivated (cf. [T]) example of a possible dependence M = M(q 2 ) 
on the momentum transfer involved in the Primakoff production (Fig. 17.1] ). Again, we use a step function 
to model the dependence on q 2 , 

©(<£it " \<1 2 \) = 77777 ^crit " <L - 2W1 - cos0 ))> ( 7 - 18 ) 



M{q 2 ) M(0) v * mt n " M(0) 

where k^, fey are the moduli of the momenta of the ALP and the photon. q m = \k</> — fey| is the smallest 
possible momentum transfer. Here, we will use the approximation m = 0, but it will be crucial to take 
into account that photons have an effective mass 

2 2 ^an e 

m y = <jj P = , (7-19) 

' m e 

so q m (u>) = u> — \/u 2 — (dp. Note that the plasma mass is crucial because it ensures that q m > 0, i.e. it 
removes ALP production processes with very small momentum transfer which would be unsuppressed. 
With this modification, Eq. (17.10P reads 



Tk 2 f +1 sm 2 9 1_ 

64^/„ 1 C ° S (x - cos 6) (y - cos 6) M 2 (q 2 



r 7-<^) =^Z I dcose 7r—^K7^—^T727^ ( 7 - 2 °) 



with x = (fe? + k 2 )/2k < pk 1 and y = x + k 2 /2k ( j,k 1 . The step function implies that only values of cos# 
satisfying 

2 2 

cos 6 > 1 ^ ~ qm (7.21) 



contribute to the integral. Hence, we find that the effect of the step function (J7.18P is to restrict the 
integration limits of Eq. (17.201) , 



Tk 2 f +1 , „ sin 2 ' 



T y _Ju>) = - ' . / dcos9- — — , (7.22) 

7 ^ V ' 647tM 2 (0)7 5H (x-cos6)(y-cos6) , V ; 



with 



S(u 



( 1 for felt < q m {u) 

1 ~ n To. T tor 9cnt > 1m\W), 1 — 7=5 = > —1 

\ - i for *»** >*»m. ^S^H- 1 ; 



(7.23) 



When <5(w) = 1, the integral is zero and Primakoff conversion is completely suppressed. This happens 
for values of the plasma frequency wp and the energy to for which the minimum momentum transfer is 
already larger than the cut-off scale q C tit- We point out that this is an energy dependent statement. For 
u S> top large enough, q m is small enough to satisfy q C rit 3> qm- When this is the case we have only 
partial suppression. The integral goes only over the small interval [S(u), 1] where 5(lo) «1- q 2 lit /2uj 2 , 



so 
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Figure 7.7: Suppression factor S for CAST, and S for the energy loss arguments as a function of (fait- 



x ~ 1 and y ~ 1 + k 2 /2u> 2 . Then the integral can be easily estimated by the value of the integrand at 
cos# = 1, 

.2 a, ,2 „2 



r 7 _ (w) 



64vrAf 2 fc 2 2u; 2 ' 



for w » wp, k s , q crit . 



(7.24) 



Notice that although we have used the strongest possible suppression, a step function, at the end of the 
day, at high energies, the transition rate is only suppressed by a factor q 2 rit /k 2 . This means that the 
7* — <f> transition is suppressed at most quadratically. 

This holds even for a generic suppressing factor F(q 2 ) = M(q 2 )/M(0). The limitation comes from the 
part of the integral which is close to cos 9 = 1. There the integrand is a constant, 1+cos 9/(y—l) ~ Au; 2 /k 2 . 
By continuity, the suppression factor F(q 2 ), whatever it is, must be close to unity because q 2 is very close 
to zero and normalization requires F(q 2 = 0) = 1. This holds for values of q 2 up to a certain range, limited 
by the shape of F(q 2 ). Defining q 2 Tit as the size of the interval where -F(|(? 2 | < Qcrit) ~ ■*-> then q 2 vit = \q 2 \ 



gives a minimum value for cos 9 for which the integrand is nearly constant (cos9 r , 
leading to 
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(7.25) 



Note that here we have ALP emission from every place of the Sun which is in contrast to the macroscopic 
environmental suppression scenario. Indeed, the ALP production rate typically increases towards the 
solar center. 

Proceeding along the lines of the previous section we can calculate the suppression factors for the 
CAST experiment S and the corresponding S that appears in the energy loss considerations. The results 
are plotted in Fig. 17.71 

Using the required suppression (17.8J) . S ~ 10~ 20 , for CAST and (17. 9|) . S ~ 10~ 10 , for the energy loss 
arguments, we infer that sufficient suppression requires 



Qcrit < 10" 2 eV. 



(7.26) 



Although this seems rather small it is nevertheless quite big compared to the typical momentum transfer 
in the PVLAS experiment, 

,2 



■m 



QPVLAS 



-i~6x 10~ 7 eV. 
2a; 



(7.27) 
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Figure 7.8: Schematic view of a "light shining through a wall" experiment. (Pseudo-)scalar production 
through photon conversion in a magnetic field (left), subsequent travel through an (opaque) wall, and 
final detection through photon regeneration (right). 

7.3.2 Kinematical suppression 

So far, we have suppressed the production of ALPs by reducing their coupling to photons. Now, we 
consider the possibility that the suppression originates from an increase of the ALP's effective mass. 
Clearly, if the latter is larger than the temperature, only the Boltzmann tail of photons with energies 
higher than the mass can contribute to ALP production. 

If we consider macroscopic environmental parameters r]{R) and, again, assume the simplest depen- 
dence on these parameters, 

m{rj < r/crit) = m (~ meV), 771(77 > r/ crit ) = 00, (7.28) 

the suppression is identical to the one computed in Sect. 17.3.11 since the Boltzmann tail vanishes for infi- 
nite mass. Accordingly, Figs. 1731 and 17.61 give the correct suppression also for the case of an environment 
dependent mass. 

Before we continue let us point out that a strong dependence of the mass on environmental parameters 
such as in Eq. (J7.28P is problematic because it requires a strong coupling between the ALP and its 
environment. This still holds even if we require only m(j] > rj cr it) > lOkeV. The strong coupling is 
likely to lead to unwanted side effects, as we commented in Sec. II, but let us however discuss some 
phenomenological aspects which could distinguish kinematical suppression from a dynamical suppression 
via the coupling. As an explicit example, we discuss a dependence on the density p. The wave equation 
for the ALP will be 

□0 + m 2 ( / o(x))0 = O. (7.29) 

The effective mass, m(p(x)), acts as a potential for <ft. This can actually lead to a new way to avoid 
the CAST bound. For example consider a situation where ALPs are emitted with energy uj. When 
they encounter a macroscopic "wall" with m{p wa \\) > uj on their way to the CAST detector, they will 
be reflected due to energy conservation (tunneling through a macroscopic barrier is negligible). In other 
words, they will not be able to reach the CAST detector and can not be observed. In this case only the 
energy loss arguments require a suppression of the production (17. 9p whereas the stronger constraint (|7.8p 
from CAST is circumvented by the reflection. 

This effect will also play a central role in the interpretation of the PVLAS result in terms of an ALP. 
Note that the interaction region (length L) of the PVLAS set up is located inside a Fabry-Perot cavity 
which enlarges the optical path of the light inside the magnetic field by a factor N r ~ 10 5 accounting for 
the number of reflections inside the cavity. In the standard ALP scenario, the ALPs created along one 
path cross the mirror and escape from the cavity. Coherent production takes place only over the length 
L. The net result produces a rotation non-linear in L but only linear in N r |50| . 

™ = »>{£*)' *>(£)■ (7 - 30) 
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However, if m = m{p) the ALPs have a potential barrier in this mirror and they will be reflected in the 
same way as the photons. In fact, the whole setup now acts like one pass through an interaction region of 
length N r L. The ALP field in the cavity will increase now non-linearly in N r L modifying the predicted 
rotation in the following way 



(*,\ 



NrLm 2 



\A9 modlfied \ = ^— 2 ) sin [-^^j , (7-31) 

where u is the frequency of the laser. For small enough m < few x 10 -6 eV this grows as 

jy2£2£j2 

\AOmodified\ ~ J" 6M 2 • (7.32) 

Under these conditions the PVLAS experiment cannot fix m using the exclusion bounds from BFRT. 
Using Eq. (|7.3ip the rotation measurement suggests, however, a much more interesting value 

Modified ~ 10 8 GeV, for m < few x KT 6 eV, (7.33) 

where we have usedL ~ 1 m, iV r ~ 10 5 and u> ~ 1 eV for the PVLAS setup. That could be reconciled 
more easily with astrophysical bounds within our framework. 

Such an effective mass will also play a role in "light shining through a wall" experiments (cf . Fig. 17. 8p . 
Typically, the wall in such an experiment will be denser than the critical density p cr it required from the 
energy loss argument. Consequently, an ALP produced on the production side of such an experiment will 
be reflected on the wall and cannot be reconverted in the detection region. Hence, such an experiment 
would observe nothing if a density dependent kinematical suppression is realized in nature. 

7.4 Summary and conclusions 

The PVLAS collaboration has reported a non-vanishing rotation of the polarization of a laser beam 
propagating through a magnetic field. The most common explanation for such a signal would be the 
existence of a light (pseudo-) scalar axion-like particle (ALP) coupled to two photons. However, the 
coupling strength required by PVLAS exceeds astrophysical constraints by many orders of magnitude. 
In this paper, we have quantitatively discussed ways to evade the astrophysical bounds by suppressing 
the production of ALPs in astrophysical environments, in particular in the Sun. 

The simplest way to suppress ALP production is to make the coupling 1/M of ALPs to photons small 
in the stellar environment. Motivated by microphysical models |116j . we considered a dependence of M 
on macroscopic environmental parameters, such as temperature, plasma mass wp, or density p whose 
values typically depend only on the distance from the center of the Sun. One of our main results is that 
it is not sufficient to suppress production in the center of the Sun only. One has to achieve efficient 
suppression also over a significant part of the more outer layers of the Sun. As apparent from Tables [TTTT 
17.21 and Eq. (|7.8|) . it is possible to reconcile the PVLAS result with the bound from the CERN Axion 
Solar Telescope (CAST) if strong suppression sets in at sufficiently low critical values of the environmental 
parameters, e.g. p ~ 1CU 3 g/cm 3 , or up ~ eV. The bounds arising from solar energy loss considerations 
are less restrictive (cf. Eq. (J7.9I) and Figs. 17.21 17.61) . As an alternative suppression mechanism, we have 
also exploited an effective mass that grows large in the solar environment. This case, too, requires that 
the effect sets in already for low critical values of the environmental parameters (cf. Figs. 17.51 17. 6|) . 

A somewhat different possibility is that the coupling 1/M depends on a microscopic parameter of 
each production event like, e.g., the momentum transfer q 2 [HE]. The typical values of such a microscopic 
parameter in the Sun may be different from those in laboratory conditions. For example, the typical q 2 
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in the Sun is ~ keV whereas it is ~ 10~ 12 eV~ in a laser experiment like PVLAS. Suppression is then 
achieved by making the coupling small for the typical values in the solar environment. However, untypical 
events occur from time to time and the required value, q CT \t < 10~ 2 eV is smaller than the q cvlt < 0.4 eV 
estimated in pQ (it is a bit more difficult to provide a similar estimate for [2] since there are energy 
loss channels other than Primakoff production). This holds even though we have used a step-function 
suppression factor which gives a stronger suppression than the form factors obtained in the microphysical 
models of [D [2]. 

Most proposed near-future experiments to test the PVLAS ALP interpretation are of the "light 
shining through a wall" type (cf. Fig. 17.80 . In these experiments, the environment, i.e. the conditions 
in the production and regeneration regions, may be modified. The above mentioned critical values are 
small enough that they may be probed in such modifications. For example, a density dependence may 
be tested by filling in buffer gas. 

In conclusion, the PVLAS signal has renewed the interest in light bosons coupled to photons. The 
astrophysical bounds, although robust, are model-dependent and may be relaxed by many orders of 
magnitude. Therefore, the upcoming laboratory experiments are very welcome and may well lead to 
exciting discoveries in a range which was thought to be excluded. 



Chapter 8 

Light scalars coupled to photons and 
non-newtonian forces 



(In collaboration with A. Dupays, E. Masso and C. Rizzo. 
Published in Phys. Rev. Lett. 98, 131802 (2007) @|) 

Spinless light particles are a common prediction of many theories that go beyond the standard 
SU(3) C x SU(2)l x U(1)y gauge theory. Probably the most famous of them is the axion [T71 [18], 
an unavoidable consequence of the introduction of a new global U(l) symmetry designed to solve the 
QCD CP-problem [151 116] . There are other examples of light particles; some are pseudoscalar like the 
axion itself and some are scalar particles. We call them axionlike particles (ALPs). 

In general, a spinless particle couples to two photons, and of course in principle it has also couplings 
to matter. But the 77 coupling is particularly interesting because many experimental searches for ALPs 
are based on it. One of these searches is based on the so-called haloscope [75], where axions or some 
other similar particles forming part of the dark matter in the galactic halo can convert into photons 
in a cavity with a strong magnetic field. Another search is to look for ALPs produced in the Sun 
converting into photons in a detector having a strong magnetic field; this is called helioscope [75]. Still 
another experiment searching for ALPs and that uses the 77 coupling looks for optical dichroism and 
birefringence in a laser that propagates in a magnetic field [33] . 

While the most recent haloscope [157] and helioscope [79] experiments have not found any signal and 
thus have put limits on model parameters, the third type of experiment described above has reported a 
positive signal. Indeed, the PVLAS collaboration finds a rotation of the polarization plane of the laser 
as well as an induced ellipticity [431 1158] . There is an exciting interpretation of these results in terms 
of a new scalar particle 0, that should have a mass m^ ~ 10 -3 eV and a ^77 coupling scale of about 
M ~ 10 5 GeV -see equation (I8.ip below for the definition of M. 

The purpose of this paper is to show that in the case that the particle <f> is indeed a light scalar, the 
1^77 coupling leads to the existence of long-range spin-independent non-newtonian forces. Our calculation 
will allow us to find a very stringent limit on the coupling, using the experimental null results coming 
from searches for new forces. 

A scalar particle couples to two photons with the lagrangian 

A = ^-iV (8.1) 

The key point is that (|8.ip induces radiatively a coupling to charged particles, for example to protons 
(see Fig 1). The induced coupling will have the form 

£ 2 = y 0f* (8.2) 
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Figure 8.1: Loop diagram. 

where \I/ is the proton field and y the Yukawa coupling. The loop diagram is logarithmically divergent. 

In order to treat physically this divergence, we notice that (|8.ip corresponds to a non-renormalizable 
term in the lagrangian and as such is expected to be valid only up to a high-energy scale A, where new 
physics has to appear. We integrate momenta in the loop of Fig.l only until A, so that we cut the 
divergence with the scale A. 

The logarithmically divergent part of the diagram of Fig.l is well-defined and is the leading radiative 
contribution to y. Approximating y by this term we obtain 

3 a m v A 

2 tt M m p 

In principle there is also a tree level Yukawa term (|8.2p in the theory, and there could be some 
cancelation between the tree level and the radiatively induced coupling (|8.3p . We regard this possibility 
as very unnatural, and will use (|8.3|) for our estimates. 

The Yukawa coupling (|8.2p leads to non-newtonian forces between two test masses m\ and m,2, due 
to (j) exchange. The total potential between 1 and 2 is 

V(r) = G^^ + y l^le-^ (8.4) 

Here rii is the total number of protons in the test body i. In (|8.4p we have neglected the new electron- 
proton and electron-electron force since the corresponding value for y in the case of the electron is smaller 
than (|8.3p by a factor m e /m p . When the two test bodies are constituted each by only one element, with 
atomic numbers Z\ and Z2, and mass numbers A\ and A2, we can approximate (|8.4p by 



v{r) ^ G H^1 



Gm\ 4tt \Aj 1 \A ' 



(8.5) 



where now in the second term inside the square brackets i.e., the term containing the correction to the 
newtonian potential, we have approximated toj ~ Aim p . 

The non-newtonian part of the potential we have obtained has two clear properties: it has a finite 
range mT and depends on the composition of the bodies, i.e. on their Z/A values. We can use the 
abundant experimental bounds on fifth- type forces to limit our parameter M as a function of m<h. 

To see how we proceed we find now the limit in the interesting ranges mT ~ (meV) -1 ~ 0.2 mm, with 
the PVLAS results in mind. Bounds on new forces have been obtained by experiments designed to measure 
very small forces. In the submillimeter range, in 1997 authors of ref. [159j using a micromechanical 
resonator designed to measure Casimir force between parallel plates gave limits on new forces, and they 
also estimated the corresponding limits on the mass and inverse coupling constant of scalar particles. 
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Figure 8.2: Constraints in the log M — m^ plane. Lines labeled Earth/moon and Astr. 1998 show 
constraints from astrophysical observations, Refs. [16311164] and [165] respectively. Lines labeled Be/Cu, 
Irvine, Eot-wash, Stanford 2, Stanford 2, and Lamoreaux show experimental constraints, Refs. |166| . 
[HZKEH], [M1[I62], [EH], [170], and HZUHS] respectively. The shaded region is excluded. 



Also, strong experimental limits on new forces have been published in |160j . where they use a torsion 
pendulum and a rotating attractor in the framework of tests of the gravitational inverse-square law. The 
most strict bounds have been obtained very recently by using torsion-balance experiments [16141162] . 
The limit from the experiment presented in ref. [1611 1162| for m^ = 1CP 3 eV is 



4-7T Gmj, 



< 1.3 x 10" 



which leads, introducing the conservative values for (Z/A) of 0.4, 



y < 7.8 x 10 



-20 



(8.6) 



(8.7) 



To obtain now a limit on M, we will simply put the value of the log in (|8.3p equal to 1; this will lead 
to a conservative limit since A 3> m p . With this, we obtain 

M > 4.2 x 10 16 GeV (8.8) 

Such a high lower bound implies that no signal of a + particle of mass m^ ~ meV should be seen 
in experiments like PVLAS [431 1158] or CAST [79]. In order to reach our conclusion we have to assume 
that the lagrangian (18. 1|) is valid up to high energy scales A 3> m p . We can proceed in analogous way 
and find bounds for other values of ml . They are shown in Fig. 2, and as we see, the limits on M are 
very tight. 

Our final discussion is about modified 1^77 vertices. Examples of such models have been developed 
in [U [2] with the motivation of making compatible the PVLAS particle interpretation with the bounds 
coming from stellar energy loss and a fortiori with the CAST results. In these models the ALP does 
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couple to new paraphotons and has not a direct coupling to photons. The (^77 vertex arises because 
there is kinetic photon-paraphoton mixing. Also, a paraphoton mass \x induces an effective photon form 
factor such that the coupling is reduced for \q\ S> (X. 

When there is a form factor with a low scale /i, the analysis shown in this paper should be modified 
accordingly. In such a case the lagrangian in (|8.ip is valid only for photons with momentum q such that 
\q\ >C fj,. The induced Yukawa will be suppressed with respect to the value (|8.3j) . and thus the lower 
bound (|8.8p can be very much relaxed if indeed /i is a low energy scale. 

In order to quantify our last assertion, we have calculated the induced Yukawa coupling when the 
photons have a modified propagator 



1 1 fi 2 



q 2 q 2 /z 2 — q 2 



(8.9) 



Now the diagram of Fig.l is finite. The calculation of the coupling to protons y' , at leading order in 
n/in p gives 

j/ = SLJL (8.10) 

y 4 M y ' 

To find the potential between bodies we have to take into account the coupling to protons as well as 
to electrons, because (I8.10P is independent of the mass of the fermion. The potential when having a form 
factor (|8.9p is given by 



V{r ) * G^™ 



1+ Gml^ UJiUv 



(8.11) 



We see that the new non-standard potential (|8.1ip has a Yukawa coupling y' that compared to (|8.3p 
is suppressed by a factor of order fi/m p . The parameter [i introduced in [TJ [2] is not fully specified by 
the theory. In order to solve the PVLAS-CAST puzzle and if we do not wish too different scales in the 
model, fj, should be in the subeV range with a preferred value fi ~ meV. For this value, the bound (|8.8j) 
would relax by 12 orders of magnitude, bringing it close to the PVLAS value M ~ 10 5 GeV. Remarkably 
enough, experiments searching for new forces and testing Casimir forces at the submm lengths may be 
sensitive to the potential (|8.11|) that corresponds to a vertex with the form factor (|8.9|) . 

Note added : The fact that a scalar-photon-photon coupling gives rise to new forces and leads to a 
bound on M has been independently realized by Shmuel Nussinov |173j . 



Chapter 9 

Final Discussion and Conclusions 



The PVLAS signal has attracted a lot of theoretical work since the publication of their signal. As 
a consequence, many arguments have arose that are relevant for my work. On the other hand it is 
mandatory to include some comments on the actual status of the field together with future prospects. In 
next section I address both tasks. Finally I present a short summary and my conclusions. 

9.1 Final discussion 

The original work of this thesis starts in Chapter EJ corresponding to reference pQ . In this paper we 
quantified the inconsistency between the PVLAS ALP interpretation and the astrophysical bounds (HB 
lifetime and Helioscope experiments) and studied two mechanisms for avoiding these bounds. Both of 
them require the existence of new particles together with the PVLAS ALP. The energy scale related to 
this new physics is closer to the suggested ALP mass m ~ meV than to the inverse coupling g" 1 ~ 10 6 
GeV. This can be achieved considering that the ALP coupling to two photons proceeds through photon- 
paraphoton mixing and a mediator MCP which couples the paraphoton field with the ALP. 

The first mechanism consists in providing a strong interaction of ALPs with the stellar medium in 
such a way that ALPs produced in stellar media can be reabsorbed inside the stellar plasma. The second 
proposes that the Primakoff production of ALPs at the relatively high temperatures of stellar interiors 
is suppressed with respect to the PVLAS experiment, performed in vacuum. Although these general 
mechanisms remain to be valid, (and have probably inspired other works trying to solve the PVLAS- 
astrophysics puzzle) the concrete models we presented are not reliable anymore. The trapping mechanism 
has to pay the price of introducing a new relatively strong interaction between low mass novel particles 
and photons. In the same article pQ we see that a background of paraphotons and ALPs will be formed in 
the early universe and this would contradict the actual observations of the cosmic microwave background, 
pointing to a transparent universe at temperatures around T ~ eV. 

Nevertheless, we decided to include this model for two reasons. First, because it allowed us to introduce 
paraphoton models as our theoretical tool to interpret the PVLAS signal as a smoking gun of low energy 
physics and second because it illustrates one of the known mechanisms to avoid the astrophysical bounds, 
so it might be inspiring for further investigations along that line. Indeed, the trapping mechanism has 
originated at least a couple of works trying to evade the stellar bounds on ALPs [Til [72] and paraphoton 
models have been further used for model building around the PVLAS signal. 

The second mechanism has proven to be more fruitful because it does not require a new strong 
interaction, but contrarily, it suggests that the new physics can come into the game to suppress the 
Primakoff production in stars, i.e. to make the ALP-photon interactions feebler. The concrete model 
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presented uses a new confining force, in the spirit of QCD, that makes the ALP a composite particle of 
millicharged preons, providing a production form factor that suppresses the stellar emission. This models 
has nowadays become obsolete. The reason is that, in this case, the main production channel of ALPs 
in stars would not be the Primakoff production but the plasmon decay into a pair of preons that would 
further hadronize into many ALPs and other preonic bound stateqj. Although the final products of the 
decay will be very different, such a process occurs at the same rate than the plasmon decay into a pair 
of MCPs I used in Chapter [3] so that the stellar energy loss bounds apply also in this case. 

In our second paper [2], presented in Chapter El we recovered the idea of the suppression of stellar 
ALP production within the framework of paraphoton models to provide a new model. There, the coupling 
between the MCP and photons proceeds through two possible intermediate states, two paraphotons, in 
such a way that they interfere negatively canceling the amplitude of any physical process when the square 
of the momentum carried by the paraphotons tends to infinity. The suppression of the stellar production 
achieved with this model can be very strong but has the disadvantage that the cancelation between the 
contributions from the two paraphotons has to be carefully tuned. Beyond this prejudice, the model is 
able not only to reconcile the PVLAS ALP with astrophysics (HB lifetime and Helioscopes) but also with 
5 th force searches if the PVLAS ALP turns out to be finally a + particle [5] . 

At the moment of writing this paper we were already aware of the above criticism to the suppression 
model presented in [JJ so this model focused on suppressing the plasmon decay into MCP pairs and not 
the Primakoff production. 

An alternative interpretation of the PVLAS signal can be made in terms of the pair production of 
millicharged particles (MCPs) with mass < 0.1 eV and charge ~ 10 -6 [27]. Basically this means that 
our model [2] is a candidate for explaining the PVLAS signal even in the absence of an ALP. Of course, 
the MCP-PVLAS interpretation based on [2] is more economical than the one originally presented in [2] 
because it requires one particle less. 

In order to know if the ALP is superfluous or not we will have to wait until more experimental data 
is available. The PVLAS experiment is sensitive to the parity of the ALP or the bosonic/fermionic 
character of the MCP. Also, a positive signal in a light-shining-through-walls (LSW) experiment could 
be the smoking gun of the existence of the ALP, since MCPs can not contribute to it. The program 
of discriminating between these two possibilities has been started in [57], and I say has been started 
because an exact study of the PVLAS signal in our model [2] has not yet been performed after we know 
that just MCP production can do the job. A combined analysis of coherent ALP production and MCP 
production is in current development. Yet a new source of rotation has been recently discovered, the 
coherent production of paraphotons in magnetic fields in models with light paracharged particles [92] . 
This adds up a new difficulty to this combined analysis because now the model [2] without an ALP will 
give a positive signal in future LSW experiments. 

Remarkably, the advent of the MCP interpretation has motivated some works proposing new bounds 
for the existence of MCPs [881 11081185] . In the last reference, the impact of MCPs in the CMB radiation 
is studied and it is nowadays the most firm constraint for the model [2]. Some work is in progress to 
ensure if the bounds presented are really model independent or if they can be circumvented in any region 
of the parameter space or by means of soft modifications of the model [2]. Notably, the CMB constraints 
do not rule out the bare ALP interpretation. 

Going back to the features of the model [2] one might feel uncomfortable not only with so many 
parameters and fields introduced but also with the apparently arbitrary choice of charge assignments. 
At this respect it is interesting to note that our model has been justified from string theory [143] . were 
the existence of additional U(l) symmetries is generally predicted as well as particles in the required 



I acknowledge interesting conversations at this respect with Joerg Jaeckel, Andreas Ringwald, Alex Pomarol, Georg 
Raffelt. 
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representations . 

In [3] I present a general study of the suppression scenario suggested in [I] . We parameterize a general 
suppression mechanism allowing the ALP mass and coupling constant to depend on any environmental 
physical quantity that can make a difference between the clean PVLAS setup and the hot and dense 
stellar plasmas. Because in the suppression scenario the Helioscope bounds are more robust than the 
energy loss bounds we focus exclusively on the case of the Sun. 

In the models worked out in [H [2] we performed a rough estimate of the stellar lifetime bounds in 
the presence of suppression of the ALP production: we computed the whole ALP emission from the 
values of the relevant parameters at the stellar center. This is indeed a quite good approximation in the 
absence of suppression but it is dangerous within our suppression models for the following reason: our 
models are designed to be more effective in suppressing ALP production wherever the conditions of the 
production are more extreme. For instance, in [1] the suppression of the Primakoff production depends on 
the screening scale, k s , which is maximum in the stellar center because the density of screening charges 
is much higher there than in outer shells. Therefore, the suppression computed is overestimated in the 
external layers of the stars. In order to provide better estimates we computed more accurate bounds by 
integrating the ALP emission over the most recent solar model available. 

In [3] I differentiate two possibilities for a general suppression scenario: (a) The sensitive parameter 
is an average environmental parameter like temperature or matter density that, given a solar model, 
depends only on the position of the Sun and (b) the parameter is the momentum transfer of the Primakoff 
conversion. Our conclusions are very interesting, in the first case the conciliation of PVLAS and CAST 
demands that the suppression of the Primakoff ALP production is complete in almost the whole volume 
of the Sun (up to a 95% of the solar radius). In the second, the scale at which the photon propagator 
must be cut-off is ~ 10~ 2 eV. It is interesting to note that this is the scale suggested in our model pQ by 
means of a much simpler estimation. 

While at the time of publication of pQ our model was the unique solution for the PVLAS-astrophysics 
puzzle some alternatives begun to show up soon after. Interestingly enough, at least two of them fit 
exactly our scheme of suppression driven by an environmental parameter in [3]. Thus our paper can be 
used to further constrain these new proposals. 

In [116] the authors propose an intricate mechanism for getting a ALP coupling g(T) that vanishes 
abruptly at a temperature between PVLAS and the stars due to the phase transition of another scalar 
field. This model fits exactly our model independent scheme [3]. Our study clearly predicts that the 
value of the critical temperature must be below 25 eV (See table I7T2J) in order to satisfy the CAST bounds 
because the ALP emission must be suppressed not only at T> keV but essentially at every place of the 
whole Sun. 

Other model has recently appear in [174] . There the authors propose that the PVLAS ALP can be 
of the chameleon type |117j . an exotic + field whose mass depends on the environmental matter density, 
m = m{p). This possibility had been indeed briefly examined in [3]. If the mass of the chameleon 
increases so much under the solar conditions such that its production is kinetically impossible then we 
have the suppression required to evade the CAST bounds. In this case the conclusions of our paper are 
much more dramatic than in the previous case. The authors use a typical matter density in the Sun 
PSun ~ 10 gem -3 while we find that the suppression should be fully active above p c ~ 9 10 -3 gem -3 . 

All the models presented in this thesis show that reconciling the PVLAS ALP with astrophysics 
(in particular with the Helioscope measurements) does require new physics at accessible scales. As I 
pointed out in Chapter H] the ALP physics must dramatically change in between the PVLAS and stellar 
environments. Notably, this suggests that the PVLAS particle (ALP or MCP) is not likely to exist alone 
and thus the new generation of low energy experiments can eventually discover several particles. 

A program to discriminate between the different possibilities is currently under way. As the conditions 
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in the solar plasma are not so extreme, they can be testable in laboratories and therefore new experimental 
efforts are crucial to test the PVLAS particle interpretation. 

Interestingly enough, many experiments have been recently proposed to test the PVLAS signal. At 
least three of them have currently began to work. The BMV and Q&A collaborations are building 
PVLAS-type experiments looking for vacuum dichroism and birefringence while the ALPS collaboration 
is currently building a LSW experiment. The same PVLAS collaboration will soon perform a LSW run 
after a soft modification of their setup. Other LSW experiments have been proposed at Jefferson Lab 
and CERN, the LIPPS and OSQAR collaborations. A detailed review of the status of the experimental 
tests of the PVLAS signal can be found in |175| . 

Finally it is worth mentioning that the PVLAS collaboration has been releasing new data in several 
conferences during the last years. Together with the published measurements of the rotation of the 
polarization plane of a lu ~ 1.2 eV laser (red) they have presented rotation data for a smaller wavelength 
light (u ~ 1.4 eV, green) and ellipticity measurements with both lasers. The absolute values of rotation 
and ellipticity data agree with an ALP interpretation but it is not the case when one takes into account 
the measured signs. Ellipticity data favors an + ALP while rotation data with the green laser favors a 0~ 
ALP. If these results are confirmed then the bare ALP interpretation would be ruled out. It remains to be 
investigated if the model in [2j could accommodate this results when the effects of MCP and paraphoton 
production are taken into account. 

The fact that PVLAS ellipticy measurements point to a + ALP motivated us to consider the 5 
bounds presented in [2j. As already noted, these bounds are the strongest constraints for a + ALP 
interpretation. However, they were unknown at the time of publishing [2] and [3J and therefore the 
requirements we derive there for a + ALP to satisfy the HB lifetime and Helioscope bounds are not 
sufficient to ensure a PVLAS ALP interpretation. As a remarkable fact, the model [2] turned out to 
evade the 5 force bounds and chameleon models as |174j are in principle able to satisfy them as well. 
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9.2 Conclusions 

In this thesis I have been concerned with some of the problems that arise when trying to interpret the 
measurements of the PVLAS collaboration [13] in terms of a novel, low mass, spin-zero particle coupled 
to light [33] we call ALP. 

My main concern has been to investigate whether the strong astrophysical bounds, excluding such a 
particle, could be circumvented in any particle physics model. I have discussed two mechanisms, trapping 
the new particles in the stars and suppressing their production in stellar environments with respect to 
the PVLAS experiment. The second has given more reliable results. 

In order to build models that either trap or suppress the number of stellar ALPs I have invoked new 
particles, namely paraphotons and millicharged particles. 

Notably the scale of the masses required for these new particles is very low < eV, a very suggestive 
value. I have also made a model-independent study of general suppression mechanisms. This work also 
points to new phenomena accessible in earth-ground laboratory experiments. Therefore, new laboratory 
experiments have been encouraged. 

In the near future, many experiments testing the PVLAS signal will be performed and this thesis 
suggests that the PVLAS particle could not be alone. Therefore, having different kinds of experiments 
will be crucial for confronting different models. 

Recently, new strong constraints on the PVLAS ALP as well as on paraphotons and millicharged 
particles have been derived. I have shown that 5 force searches set strong limits on the PVLAS particle 
if it is a + . On the other hand, the CMB has been proven to be sensitive to the existence of millicharged 
particles. Notably, my best model is able to avoid both constraints. 

If the PVLAS signal is confirmed and firmly established, I believe that few particle physics models 
could solve the strong conflict with astrophysics. In these thesis I have developed solutions that involve 
physics at very low energies. 



Appendix A 

Schizons 



If we let our ALP <p s to couple to light through both interactions in (|2.4|) parity is not conserved by 
the ALP interactions. Therefore, the ALP has not definite parity and because of that it can be called a 
"schizon" |176j . In this case, we find a complicated 3x3 mixing problem 
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Where in the relativistic limit ? = ±1 for photons traveling along the zLz-direction. In absence of Cotton- 
Mouton and faraday effects (vacuum) the problem simplifies notably. By means of a rotation Rl of angle 
91 in the Am — A± space it can be decoupled into 
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where A s ,A n are the photon states coupling to the schizon and neutral respectively, and g = 

is the effective coupling constant. We have found exactly equation (|2.16p with the following replacements 
{s, n} <-> {||, _L}, so equations (|2.23p . (|2.24p are also solutions for the schizon case. In order to give a result 
for the induced rotation and ellipticity in PVLAS we have to rotate back to Aim . In the weak mixing 
regime, I get 
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The schizon couplings add at least two new ingredients to the ALP interpretation. First, they provide a 
phase shift between the photon initial polarization and the final rotation. By sending gn j_ to zero we recover 
the two pure parity cases, because in this cases sin6 s = 0, 1, namely 9 S = 0,ir/2 and sin(2£o — 0,ir) = 
±sin2£o. The rotations shown again to be of opposite sign. Secondly, the photon state coupling to <fi s 
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depends on the direction of propagation (91 depends on the sign ?, namely 9 S = —9f) 

(j>s(— ►) couples with g A± + g~Au (A. 5) 

(j) s (<—) couples with — g + A±+g~A\\ . (A. 6) 

As a consequence the effect of multiple reflections inside a Fabry- Perot cavity can be considerably dimin- 
ished |177j . a fact that worsens the PVLAS discrepancy with astrophysics. This can be easily seen from 
(|A.3|) and (|A.4p since the sum of the two rotations(ellipticities) acquired after two consecutive passes in 
opposite directions is 

*C(<-) = 2cos20 s x ^!^srn 2 ^£ (A.7) 

mi 4w 

which is maximum for 9 S = 0, 7r/2..., namely, the pure parity cases. A similar expression holds for 
the ellipticity. A remarkably fact in this case is that the information of g and 9 S cannot be extracted 
separately and thus PVLAS is not able to discriminate is the ALP is an schizon or not. 
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